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was  the  slowest  step  in  the  process,  requiring  90-95  percent  of  the  drop  life¬ 
time,  even  in  regions  where  maximum  agglomerate  reaction  rates  were  observed. 

An  analysis  was  developed  to  provide  predictions  of  both  liquid  end  agglomerate 
heat-up  and  reaction.  Area-reactivity  factors  were  used  in  the  carbon-agglom¬ 
erate  reaction  analysis  to  account  for  the  effects  of  pores,  reaction  within 
the  particle,  intrinsic  variations  of  carbon-black  reactivity  and  catalysts. 
Transport-rate  enhancement  factors  were  also  employed  in  the  carbon-agglomerate 
reaction  analysis  to  account  for  the  increased  mass  and  energy  transport  due  to 
the  flow  percolating  through  the  open  porous  structure  of  the  carbon  agglomerate. 
Two  separate  models  were  examined.  The  first  considered  oxygen  and  hydroxal 
reacting  with  carbon,  and  the  second  involved  oxygen,  carbon  dioxide  and  water 
vapor  reacting  with  carbon.  The  analysis  yielded  good  predictions  of  both 
particle  size  and  temperature  variations  for  flame  equivalence  ratios  of  0.272- 
1.350.  The  use  of  a  catalyzed  slurry  was  found  to  increase  agglomerate  burning 
rates  in  the  lean  portions  of  the  flame,  extending  the  lean  limit  of  agglomerate 
reaction. 

A  1  minar-premixed  flame  apparatus  was  also  constructed,  in  order  to  examine 
the  variation  of  the  empirical  factors  over  a  broader  range  of  test  conditions. 
These  tests  emphasized  smaller  freely  moving  particles,  more  representative  of 
practical  slurry  sprays,  and  test  conditions  in  the  kinitically-controlled 
‘regime.  The  apparatus  consisted  of  a  particle  gene  :ator  which  produced  a  stream 
of  carbon  agglomerate  particles,  with  initial  initial  diameters  in  the  10-100  Vm 
range,  which  passed  through  the  post  flame  region  of  a  flat-flame  burner. 
Measurements  were  made  of  particle  size,  mass,  temperature  and  velocity  as  a 
function  of  position  above  the  burner. 

A  relationship  between  apparent  particle  density  and  reacted-mass  fraction 
was  obtained  which  was  relatively  independent  of  particle  size  and  local  particle 
environment.  The  variation  of  both  the  transport-enhancement  factors  and  the 
area-reactivity  multiplication  factors,  for  both  reaction  mechanisms,  were 
found  to  depend  on  particle  diameter  and  reacted-mass  fraction  to  some  degree. 
However,  constant  values  similar  to  those  found  for  the  large  supported 
particles,  also  yielded  a  reasonable  correlation  of  the  measurements.  The 
burning  rate  increased  with  increasing  flame  temperature  and  with  decreasing 
equivalence  ratio. 
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ABSTRACT 

A  theoretical  and  experimental  investigation  of  the  combustion 
process  of  carbon  slurry  fuels  is  described.  The  combustion  of 
individual  drops  (400-1000  pm  in  diameter)  supported  at  various 
positions  within  an  open  turbulent  diffusion  flame  was  observed;  When 
a  slurry  drop  was  exposed  to  the  flame,  the  liquid  fuel  evaporated  in 
the  first  stage  of  the  process  leaving  a  porous  carbon  agglomerate 
formed  from  the  carbon  particles  in  the  slurry.  The  second  stage 
involved  heat-up  and  reaction  or  quenching  of  the  agglomerate. 
Consumption  of  the  agglomerate  was  the  slowest  step  in  the  process, 
requiring  90-95  percent  of  the  drop  lifetime,  even  in  regions  where 
maximum  agglomerate  reaction  rates  were  observed.  An  analysis  was 
developed  to  provide  predictions  of  both  liquid  and  agglomerate  heat¬ 
up  and  reaction.  Area- reactivity  factors  were  used  in  the  carbon- 
agglomerate  reaction  analysis  to  account  for  the  effects  of  pores, 
reaction  within  the  particle,  intrinsic  variations  of  carbon-black 
reactivity,  and  catalysts.  Transport- rate  enhancement  factors  were 
also  employed  in  the  carbon-agglomerate  reaction  analysis  to  account 
for  the  increased  mass  and  energy  transport  due  to  the  flow 
percolating  through  the  open  porous  structure  of  the  carbon 
agglomerate.  Two  separate  models  were  examined.  The  first  considered 
oxyg'n  and  hydroxyl  reacting  with  carbon,  and  the  second  involved 
oxygen,  carbon  dioxide  and  water  vapor  reacting  with  carbon.  The 
analysis  yielded  good  predictions  of  both  particle  size  and 
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temperature  variations  for  flame  equivalence  ratios  of  0.272-1.350. 
The  use  of  a  catalyzed  slurry  was  ' found  to  increase  agglomerate 
burning  rates  in  the  lean  portions  of  the  flame,  extending  the  lean 
limit  of  agglomerate  reaction. 

A  laminar- premixed  flame  apparatus  was  also  constucted,  to 
examine  the  variation  of  the  empirical  factors  over  a  broader  range  of 
test  conditions.  These  tests  emphasized  smaller  freely  moving 
particles,  more  representative  of  practical  slurry  sprays,  and  test 
conditions  in  the  kinetically-controlled  regime.  The  apparatus 
consisted  of  a  particle  generator  which  produced  a  stream  of  carbon 
agglomerate  particles,  with  initial  initial  diameters  in  the  10-100  |io 
range,  which  passed  through  the  post-flame  region  of  a  flat-flame 
burner.  Measurements  were  made  of  particle  size,  mass,  temperature 
and  velocity  as  a  function  of  position  above  the  burner. 

A  relationship  between  apparent  particle  density  and  reacted-mass 
fraction  was  obtained  which  was  relatively  independent  of  particle 
size  and  local  particle  environment.  The  variation  of  both  the 
transport-enhancement  factors  and  the  area-reactivity  multiplication 
factors,  for  both  reaction  mechanisms,  were  found  to  depend  on 
particle  diameter  and  reacted-mass  fraction  to  some  degree.  However, 
constant  values  similar  to  those  found  for  the  large  supported 
particles,  also  yielded  a  reasonable  correlation  of  the  measurements. 
The  burning  rate  increased  with  increasing  flame  temperature  and  with 
decreasing  equivalence  ratio. 
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CHAPTER  I 


INTRODUCTION 

1 . 1  General  S tatement  of  the  Problem 

A  promising  technology  for  coal  utilization  is  to-  mix  finely 
crushed  carbon  or  coal  in  oil  and  then  directly  burn  the  resulting 
slurry  in  stationary  combustors  such  as  furnaces  and  boilers.  This- 
concept  is  attractive  since  it  directly  substitutes  coal  for  oil  as  an 
energy  source.  Also,  the  crushing  and  mixing  processes  involved  in 
the  slurry  preparation  are  far  more  economical  than  energy-expensive 
coal-liquefaction  and  gasification  processes.  Finally,  the  slurry  is 
pumpable  and  can  readily  be  used  in  conventional  combustors  with 
minimal  equipment  modifications. 

In  addition  to  economy,  the  high  density  of  the  combustible  solid 
carbon  gives  the  slurry  a  significantly  higher  volumetric  energy 
density  than  conventional  liquid  fuels.  Because  of  its  high  energy 
density,  a  carbon  slurry  is  a  very  promising  military  fuel  for  air- 
breathing  propulsion  systems,  such  as  cruise  missiles,  where 
compactness  Is  an  important  design  criterion  [1]. 

In  order  to  meet  the  high  energy  density  requirements  of  cruise- 
missile  applications  a  low-volatile,  low-ash  solid  carbon  fuel 
suspended  in  a  high  energy  density  liquid  fuel  was  desired.  The  fuel 
proposed  was  a  medium- thermal  carbon  black — liquid  JP-1Q  mixture. 
JP-LG  or  exo-tetranyarodi  (cyclopentadiene)  is  a  pure  hydrocarbon  with 
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the  chemical  formula  of  C^qH^,  and  its  properties  and  structure  are 
given  in  Appendix  A.  A  50%  by  mass  solid — 50%  by  mass  liquid  fuel 
mixture  is  believed  to  be  a  good  compromise  between  high  energy 
density,  greater  than  50,000  MJ/m ■*,  and  low  viscosity  to  maintain  the 
fuel's  liquid  characteristics  [2].  The  formulation  would  also  have  to 
contain  additives  to  deflocculate  the  thermal  black,  stabilize  the 
dispersion,  and  promote  combustion  [2]. 

The  general  combustion  properties  of  carbon  slurries  have  been 
examined  in  gas  turbine  combustors  [2,3]  and  well-stirred  reactors 
[4].  The  substitution  of  slurries  for  liquid  fuels  in  combustors, 
however,  is  not  without  problems-  The  liquid  fuel  tends  to  bum  off 
before  the  medium- thermal  black  can  ignite,  especially  in  low- 
temperature  combustors  [2].  Furthermore,  the  carbon  residue  remaining 
after  the  liquid  evaporates  is  difficult  to  burn  out.  Therefore,  new 
combustion-chamber  designs  optimizing  the  combustor  aerodynamics  are 
required  to  accommodate  the  dual  burning  nature  of  slurry  fuels  [2]. 

The  present  study  emphasized  fundamental  aspects  of  the  carbon- 
slurry  combustion  problem  by  observing  and  analysing  the  combustion 
properties  of  individual  slurry  drops.  The  droplet  combustion  process 
is  fundamental  to  successful  operation  of  a  combustor  since  it  affects 
the  residence-time  requirements  and  combustion  efficiency.  Thus,  the 
combustion  of  carbon-slurry  droplets  has  been  studied  under  conditions 


which  simulate  those  of  a  combustion  chamber. 


1.2  Description  of  Slurry  Drop  Combustion 

Slurries  have  not  been  widely  used  in  Che  pasc  and  Che  existing 
literature  on  their  combustion  characteristics  is  limited.  A  recent 
paper  by  Law  et  al.  [5]  reports  some  results  on  the  combustion 
characteristics  of  single,  supported  coal  slurry  drops.  Combustion 
was  accomplished  in  still  air  under  natural  convection  conditions. 
The  slurry  consisted  of  dried,  200  mesh,  coal  powder  (coal  type 
unidentified)  in  No.  2  heating  oil.  The  drops  were  suspended  from  a 
small  probe  (initial  drop  diameters  of  850-1250  Dm)  and  ignited.  The 
initial  combustion  process  was  fueled  by  the  liquid.  As  the  liquid 
continued  to  evaporate,  small  masses  of  material  (presumably 
containing  coal  particles)  were  intermittently  ejected  from  the  drop. 
Ai'-ev  all  the  oil  had  vaporized,  the  flame  extinguished,  leaving  an 
irregularly  shaped  agglomerate  of  coal  particles  on  the  support.  The 
mass  of  material  remaining  on  the  probe  increased  as  the  initial 
percentage  of  coal  in  the  slurry  was  increased.  For  initial  coal 
weight  percentages  greater  than  15-20%,  the  mass  of  the  coal 
agglomerate  remaining  after  extinction  reached  70-80%  of  the  initial 
mass  of  coal  in  the  drop. 

Based  on  existing  information  on  coal  combustion  (6-10],  it  is 
likely  that  persistent  combustion  of  a  coal-slurry  drop  in  air 
involves  three  major  processes:  (1)  evaporation  and  combustion  of  the 
liquid  fuel,  (2)  devolatilization  of  the  coal  with  subsequent 
oxidation,  and  (3)  reaction  of  the  remaining  solid  carbon-ash 
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as  che  solid  phase;  thus^ a  devolatilization  stage  is  not  expected. 
Judging  from  the  results  or  Law  et  al.  [5],  carbon  levels  in  the 
present  test  slurries  are  sufficiently  high  to  result  in  significant 
agglomeration. 

1.2.1  Preliminary  Test  Apparatus 

In  order  to  gain  further  insight  into  the  combustion  processes 
of  a  carbon-slurry  droplet,  a  single  fuel  drop  (50.42  dispersed  carbon 
by  weight  in  JP-10)  suspended  from  a  quartz  fiber  in  still  air  was 
ignited  with  a  match.  As  a  baseline  a  pure  JP-10  droplet  was  also 
burned  in  the  same  manner.  The  drop  combustion  process  was  observed 
with  a  16  mm  Photosonics,  model  16-B,  motion  picture  camera.  The 
camera  optics  provided  a  2:1  magnification.  The  camera  was  powered 
with  a  Kepco,  SM  36-5  AH  d.c.  power  supply.  The  framing  rate  was 
indicated  with  a  timing  light  on  the  camera,  activated  with  an  Adtrol 
Electronic  pulse  generator,  model  501.  Kodak  Plus~X  reversal  film  was 
used  for  all  tests. 

Backlighting  for  shadowgraph  measurements  was  provided  by  a  Pek, 
model  401A,  75W  mercury  arc  lamp.  The  light  from  the  arc  was 
collimated  and  directed  toward  the  drop  location.  A  diffuser  screen 
was  employed  behind  the  drop  to  equalize  the  light  intensity  of  the 
background.  The  background  intensity  was  adjusted  so  that  envelope 
flames  around  the  drop  could  also  be  observed. 

The  film  records  were  analyzed  on  a  frame-by-frame  basis,  using  a 
Vanguard  motion  picture  analyzer.  Photographs  of  objects  of  known 
size  at  the  drop  location  provided  a  calilraf.ion  of  distances  on  the 
film. 
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Carbon  residue  lefc  on  che  probe  after  combustion  was  examined 
with  a  scanning  electron  microscope  (SEM),  International  Scientific 
Instruments,  model  M-7.  Permanent  records  of  the  appearance  of  the 
residue  were  obtained  with  Polaroid  film. 

1.2.2  Preliminary  Test  Observations 

The  pure  JP-10  drop  burned  in  a  conventional  manner  for  a 
liquid  fuel  [11].  An  envelope  flame  was  observed  around  the  droplet 
until  all  the  fuel  was  consumed.  After  an  Initial  heat-up  period,  the 
square  of  the  droplet  diameter  decreased  in  a  linear  fashion,  which  is 
typical  behavior  for  supported  drops  under  natural  convection 
conditions. 

The  initial  phase  of  combustion  of  the  slurry  drops  was  similar 
to  that  of  pure  JP-10.  An  envelope  flame  was  observed  around  the  drop 
and  the  drop  diameter  decreased.  This  phase  appeared  to  involve  the 
combustion  of  the  JP-10  in  the  slurry.  However,  after  a  time,  the 
envelope  flame  was  extinguished  and  a  black  residue  remained  on  the 
probe.  Although  gravimetric  measurements  were  not  made,  it  appeared 
that  most  of  the  original  solid  carbon  in  the  drop  remained  unburned. 

A  slurry  with  a  catalyst  to  initiate  carbon  combustion  was  also 
burned.  The  catalyzed  slurry  drops  exhibited  the  most  interesting 
combustion  behavior  in  air.  A  typical  dark  field  motion  picture  for 
this  case  appears  in  Figure  l.  The  early  frames  clearly  show  the 
envelope  flame  around  the  drop  persisting  for  nearly  a  second.  After 
this  flame  was  extinguished,  a  small  amount  of  luminosity  was  present 
for  a  period  of  2.8  s.  Then  the  residue  itself  began  to  glow. 


indicating  reaction  of  the  carbon.  After  a  period  of  glowing  (12.6  s, 
of  which  only  the  last  2.5  s  are  shown  in  the  figure)  the  particle 
extinguished  again,  which  left  a  carbon  residue  suspended  from  the 
probe. 

The  dark  period  prior  to  the  glowing  period  probably  involved 
reaction  of  the  carbon,  causing  it  to  become  heated.  -Glowing  was 
observed  where  heating  has  proceeded  to  a  temperature  high  enough  for 
the  particle  to  emit  significant  radiation  in  the  wavelength  range  for 
which  the  photographic  system  is  sensitive.  Thus  the  absence  of 
glowing  is  not  necessarily  indicative  of  no  reaction,  and  even  the 
uncatalyzed  drop  may  have  undergone  significant  reaction  during  the 
dark,  period. 

These  results  indicate  that  the  catalyst  has  an  important  effect 
on  the  carbon  reaction  rate.  The  catalyst  increased  the  reactivity  of 
the  carbon,  and  thus  combustion  was  sustained  at  higher  rates  in  spite 
of  heat  losses  te  the  surroundings,  at  least  for  a  time.  The  slow 
step  in  the  reaction  process  involved  combustion  of  the  carbon,  which 
will  be  the  controlling  step  during  spray  combustion  in  a  practical 
combustion  chamber. 

These  observations  indicate  that  the  overall  behavior  of  the 
present  carbon-slurry  drops  was  qualitatively  similar  to  the  behavior 
obser/ed  for  coal  slurries  [5J.  In  the  present  case,  however,  there 
was  less  disruption  of  the  particles  during  combustion,  probably  due 
to  the  absence  of  the  volatilization  step  for  the  solid. 
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Further  insight  concerning  the  slurry  combustion  -"echanism  was 
obtained  by  observing  the  carbon  residue  with  che  SEM.  As  a  baseline. 
Figure  2a  is  a  photograph  of  the  residue  from  a  catalyzed  slurry  drop, 
which  was  evaporated  in  air  without  combustion.  In  this  case,  the 
suriace  of  the  agglomerate  is  very  smooth.  The  smoothness  suggests 
that  the  particles  in  the  slurry  tend  to  pack,  together-  in  a  close 
array,  with  perhaps  the  smallest  sized  particles  more  prevalent  at  the 
surface . 

A  SEM  photograph  of  the  uncatalyzed  slurry  residue,  after 
combustion  in  air,  is  shown  in  Figure  2b.  In  this  case,  the  residue 
appears  smooth  in  some  regions  and  rough  in  others.  The  residue  also 
appears  to  be  partly  hollow.  It  was  difficult  to  prevent  some 
disturbance  of  the  sample  while  it  was  transported  to  and  mounted  in 
the  SEM;  therefore,  the  fractured  zones  may  be  due  to  handling 
problems.  Certainly  some  areas  are  very  similar  to  the  appearance  of 
the  surface  when  no  reaction  occurred,  e.g.,  as  shown  in  Figure  2a. 
The  rough  areas  could  be  typical  of  the  subsurface  structure  in  cases 
of  little  reaction  or  could  result  from  consumption  of  smaller 
particles. 

Figure  3a  is  an  illustration  of  SEM  photographs  of  the  residue  of 
the  catalyzed  slurry  drop  following  combustion  in  air.  This  structure 
is  very  different  from  the  preceding  cases.  The  residue  has  a  more 
open  structure  and  no  smooth  regions  are  observed.  The  reaction  of 
the  solid  appears  to  extend  a  significant  distance  below  the  outermost 
surface  of  the  particle.  The  smallest  elements, 


observed  at  a 
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magnification  of  10,000,  appear  to  have  a  characteristic  diameter  of 
approximately  0.3  pm,  which  corresponds  to  the  ultimate  particle  size 
of  the  slurry.  These  sub-elements  are  arranged  in  groups  having  a 
grape-cluster  appearance.  The  protuberance  visible  in  Figure  3b  has  a 
diameter  of  roughly  25  pm,  which  would  be  more  typical  of  the  residue 
from  a  drop  produced  under  combustion-chamber  conditions.  It  is 
evident  that  particles  of  this  size  are  irregular  in  shape  and 
relatively  porous. 

1.2.3  Burning  Rates  in  Air 

Figure  4  is  an  illustration  of  the  variation  of  particle 
diameter  as  a  function  of  time  after  ignition  iu  air  for  three  test 
fuels:  pure  JP-10,  noncatalyzed  slurry  and  catalyzed  slurry. 
Throughout  the  suspended-particle  test  of  this  investigation,  the 
reported  particle  diameter  is  the  equivalent  diameter  of  an  ellipsoid 
having  the  same  volume,  i.e.. 


d  -  (d.i„2  Wl/3  u-i) 

where  d  is  the  reported  diameter,  and  d__„  and  d_,-„  are  the  major  and 

uldX  IwxU 

minor  diameters  of  the  particle. 

After  ignition,  the  diameter  of  the  pure  JP-10  droplet  decreased 
at  an  accelerated  rate  until  all  of  the  fuel  was  consumed.  An 
envelope  flame  was  observed  to  surround  the  liquid  droplet  as  the 
diameter  was  decreasing.  The  slurry  drops  also  exhibited  a  reduction 
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Droplet  diameter  as  a  function  of  time  for  pure  , 
noncatalyzed  and  catalyzed  slurry  droplets  burnei 


in  diameter  in  the  region  just  after  ignition,  when  an  envelope  flame 
was  observed.  This  initial  diameter  reduction  was  due  to  the 
combustion  of  liquid  JP-1Q,  and  the  rate  of  diameter  decrease  was 
approximately  equal  to  that  for  the  pure  JP-10  droplet. 

Once  the  envelope  flame  of  the  noncatalyzed  slurry  drop 
extinguished,  no  further  luminosity  was  observed  and  the  particre 
diameter  remained  constant.  With  a  catalyst  present,  however,  the 
particle  glowed  after  a  dark  period,  and  then  extinguished  again 
leaving  a  residue.  During  the  glowing  period,  however,  there  was 
relatively  little  variation  in  the  diameter  of  the  particle.  This 
agrees  with  the  SEM  observations  where  the  structure  was  rather  porous 
following  combustion  of  a  catalyzed  slurry  drop,  suggesting  reaction 
in  depth  rather  than  just  at  the  surface.  The  bulk  of  the  reaction  in 
the  glowing  period  is  probably  confined  to  the  subsurface  region  where 
heat  losses  are  reduced.  As  the  reaction  proceeds,  the  structure 
becomes  more  porous  while  the  outer  portions  remain  in  place.  Final 
extinguishment  then  becomes  a  complex  process  dependent  upon  heat 
losses,  oxidant  penetration  into  the  pores  of  the  structure,  and 
perhaps  loss  of  the  catalyst. 

The  variation  of  droplet  diameter  observed  in  Figure  4  can  be 
explained  with  the  combustion  mechanism  of  a  slurry  droplet  presented 
in  Figure  5.  After  ignition  and  the  establishment  of  an  envelope 
flame  surrounding  the  drop,  the  droplet  begins  to  heat  up.  As  the 
liquid  temperature  increases,  the  rate  of  evaporation  increases  due  to 
the  increased  vapor  pressure  of  the  fuel  at  the  liquid  surface. 


Sketch  illustrating  the  details  of  the  slurry  combustion  process 
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Eventually  a  condition  is  reached  where  all  the  heat  transferred  from 
the  flame  to  the  drop  surface  is  utilized  for  the  heat  of  vaporization 
of  the  evaporating  fuel,  and  the  liquid  temperature  stabilizes  at  what 
is  termed  the  wet“bulb  temperature.  The  wet-bulb  temperature  is 
generally  somewhat  below  the  boiling  temperature  of  the  fuel  at  the 
ambient  pressure.  The  vaporization  rate  is  greatest  at  the  wet-bulb 
temperature.  Therefore,  while  the  drop  diameter  decreases  through  the 
period  of  liquid  combustion,  the  rate  of  decrease  is  greatest  at  the 
wet-bulb  condition. 

As  the  liquid  evaporates,  the  solid  carbon  particles  pack 
together  to  form  an  agglomerate.  Eventually,  the  remaining  liquid 
fuel  can  only  be  found  in  the  pores  of  the  agglomerate,  and  although 
liquid  is  evaporating  and  fueling  the  flame,  the  diameter  of  the 
agglomerate  remains  constant.  When  all  the  liquid  fuel  has  evaporated 
the  envelope  flame  is  extinguished.  At  this  point  the  particle  will 
either  extinguish  as  in  the  noncatalyzed  slurry  case;  or  the 
temperature  will  begin  to  rise  very  slowly  as  in  the  catalyzed  slurry 
case,  due  to  Ignition  of  the  carbon  by  the  hot  gases  produced  by  the 
envelope  flame.  As  the  carbon  particle  heats  up,  the  rate  of  reaction 
increases.  At  sufficiently  high  temperatures,  the  surface  layer  of 
the  carbon  agglomerate  vigorously  combusts  which  causes  the  diameter 
of  the  particle  to  decrease.  As  the  process  continues,  the  particle 
will  either  become  small  enough  to  be  extinguished  by  natural 
convection  effects,  or  will  burn  until  there  is  no  carbon  remaining  on 


the  support  probe. 
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1.3  Previous  Related  Studies 

The  preliminary  tests  suggested  that  the  life-history  of  a 
combusting  slurry  droplet  could  be  divided  into  two  separate  phases, 
the  liquid  heat-up  and  evaporation  phase,  followed  by  the  agglomerate 
heat-up  and  reaction  phase.  Therefore,  it  appears  logical  to  model 
the  slurry  droplet  combustion  process  as  two  independent  stages.  The 
review  of  previous  related  studies  is  separated  in  the  same  manner. 

1.3.1  Evaporation  of  Liquid  Droplets 

Drop-life-histories  must  be  calculated  in  order  to  determine 
interphase  heat,  mass  and  momentum  transfer.  In  determining  the  drcp- 
life  history,  the  drop  configuration,  convection  effects,  transient 
effects,  fluid  properties,  interface  effects,  radial  pressure 
gradient,  radiation  effects,  constitutive  relations  for  fluxes,  and 
whether  the  droplet  is  evaporating  or  combusting  all  play  a 
significant  role.  A  recent  review  of  drop-life-history  assumptions 
and  models  is  presented  by  Faeth  [11]. 

Most  studies  have  considered  a  drop  to  be  a  smooth  sphere 
although  deformation  of  the  drop  into  a  flattened  elliptical  shape 
occurs  at  higher  Reynolds  numbers  [12,13]. 

The  spherical  case  itself  presents  substantial  analytical 
difficulties  since  the  presence  of  convection  requires  numerical 
solutions,  except  at  the  limit  of  low  Reynolds  numbers.  Results  are 
available,  for  low  Reynolds  numbers  flow  over  a  rigid  sphere  [14, 15], 
and  some  studies  have  considered  higher  Reynolds  numbers  and 
circulating  liquid  drops  [16-13] . 


17 


Complete  solutions  for  drop  characteristics  are  of  value  to  check 
acre  approximate  methods,  but  they  are  too  complex  to  be  useful  for 
drop-life-history  analysis  at  tnis  time.  The  general  approach  has 
been  to  model  the  drops  as  a  spherically  symmetric  flow  field,  and 
then  employ  an  empirical  correction  for  convection  [19-26].  This 
procedure  frequently  involves  the  use  of  film  theory,-  where  the 
transport  process  is  analyzed  in  a  spherically  symmetric  layer  about 
the  droplet.  Only  the  radial  motion  due  to  mass  transfer  is 
considered:  convection  is  handled  by  selecting  the  film  thickness  to 
yield  the  correct  heat  or  mass  transfer  rate,  at  a  given  Reynolds 
number,  at  the  limit  of  low  mass-transfer  rates  [11]. 

During  its  lifetime  a  droplet  is  heating  up;  its  radius  is 
changing;  due  to  drag,  its  relative  velocity  is  changing;  and  the 
ambient  conditions  may  be  changing  as  the  droplet  travels  through  the 
flow  field.  Therefore,  both  gas-  and  liquid-phase  transient  effects 
must  be  considered  for  exact  analysis.  However,  just  as  in  the  case 
of  convection  for  a  sphere,  a  complete  transient  solution  is  too 
unwieldy  to  be  of  value  in  drop-life-history  analysis,  and 
simplifications  must  be  sought - 

At  high  Reynolds  number,  or  at  low  Reynolds  numbers  3nd  moderate 
pressures  where  the  thermal  diffusivity  of  the  gas  phase  is  two  to 
three  orders  of  magnitude  larger  than  the  thermal  diffusivity  of 
liquid  phase,  the  gas-phase  transient  time  is  much  smaller  than  the 
liquid-phase  transient  ti-w  [11].  Thus,  it  is  reasonable  to  assume  a 
quasisteady  gas  phase  where  the  boundary  layer  around  the  droplet  has 
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the  saae  properties  as  a  steady  boundary  layer  for  the  sane  conditions 
of  droplet  size,  velocity,  surface  and  ambient  temperatures,  etc.;  and 
the  effect  of  the  radial  velocity  of  the  liquid  surface  can  be 
neglected . 

Comparison  between  detailed  numerical  calculations  of  the 
transient  evaporation  process  and  the  quasisteady  gas-phase  model  have 
been  completed  by  Hubbard,  Denny,  and  Mills  [27}.  The  calculations 
demonstrate  that  the  quasisteady  gas-phase  approximation  yields  small 
errors,  at  least  for  pressures  below  10  bar.  Bellan  and  Summerfield 
[28}  present  some  simplified  rules  for  determining  the  region  where 
the  quaslstaady  approxiaa^ion  is  valid,  concluding  that  transient  gas- 
phase  effects  are  small  at  low  pressures.  Several  other  investigators 
also  present  discussions  on  various  aspects  of  the  quasisteady  gas- 
phase  assumptions  with  generally  similar  conclusions  [25,26,29]. 

For  volatile  fuels  at  low  pressures,  the  wet-bulb  temperature  is 
close  to  typical  Injection  temperatures,  and  liquid-phase  properties 
can  be  assumed  to  be  constant.  At  higher  pressures,  however,  reduced 
liquid  densities  cause  the  droplet  to  swell,  and  variable  properties 
must  be  considered  for  accurate  work  (30J. 

A  number  of  analyses  have  been  conducted  for  variable  gas-phase 
property  effects  on  the  evaporation  process  [27, >1,32].  While  early 
variable  property  studies  only  c>nsidered  variable  temperature 
effects,  recent  work  has  also  Created  variable  concentration  effects, 
since  the  effect  of  species  concentrations  can  be  appreciable  [27], 


Conditions  at  the  droplet  surface  must  also  be  specified;  the 
formulation  requires  an  estimation  of  the  radial  mass  flux  of  the 
ambient  gas  at  this  position,  as  well  as  a  relationship  between  the 
temperature  of  the  liquid  surface  and  the  concentration  of  the  fuel  in 
the  gas  phase  adjacent  to  the  surface.  At  moderate  pressures,  where 
the  quasisteady  gas-phase  assumption  is  valid,  liquid  fuels  have  low 
solubilities  for  the  nonfuel  portion  of  the  ambient  gas.  Therefore, 
Spalding's  [33]  insolubility  condition  is  most  often  employed  as  a 
boundary  condition,  and  the  mass  flux  of  the  ambient  gas  is  taken  to 
be  zero  at  the  liquid  surface. 

The  relationship  between  the  liquid  temperature  and  the  gas-phasa 
fuel  concentration  at  the  liquid  surface  is  usually  obtained  from  the 
equilibrium  vapor  pressure  of  the  pure  fuel  as  a  function  of 
temperature.  Vapor-pressure  data  can  be  correlated  with  the  Ciausius- 
Clapeyron  equation.  The  assumption  implied  by  the  use  of  pure-vapor- 
pressure  data  is  that  the  presence  of  other  gases  at  the  surface  does 
not  affect  the  equilibrium  vapor  pressure.  This  is  adequate  as  long 
as  the  liquid-phase  solubilities  of  the  ambient  gases  are  low  and 
pressure  levels  are  moderate,  i.e.,  in  the  range  where  an  ideal  gas 
model  is  valid  near  the  liquid  surface. 

Models  of  spherically  symmetric  evaporation  usually  ignore  the 
radial  momentum  equation  and  assume  that  the  pressure  is  uniform.  The 
bulk,  radial  flow,  however,  causes  a  radial  pressure  gradient  since  the 
mean  radial  velocity  decreases  with  increasing  radial  distance.  The 
worst  condition  occurs  during  the  steady  evaporation  period,  when 
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radial  velocities  are  largest.  However,  the  pressure  can  be  assumed 
to  be  nearly  constant,  except  for  drops  on  the  order  of  1  pa  or  less, 
at  pressures  above  atmosphere  pressure  [11].  Further  discussion  of 
radial  pressure  gradient  can  be  found  in  Reference  27. 

Previous  studies  of  burning  drops  have  indicated  that  radiative 
heating  is  relatively  unimportant  for  drop  sizes  representative  of 
most  sprays  [34-38].  Beriad  and  Hibbard  [36]  compared  the  measured 
absorption  spectrum  of  113  jjm  thick,  layers  of  typical  hydrocarbon 
fuels,  iso-octane,  cyclohexane,  toluene,  hexene,  and  JP-4,  with 
measured  emission  spectra  of  "luminous  and  nonluminous"  burner  flames. 
The  fuels  generally  absorb  radiation  near  3.5  ra,  where  flame  emission 
was  low  in  both  cases.  They  conclude  that  these  fuels  absorb  little 
radiative  energy  from  the  combustion  regions  of  gas  turbine 
combustors- 

Droplet  evaporation  involves  heat  and  mass  transfer  in  the 
presence  of  temperature  gradients  and  multi-component  mixtures.  The 
appropriate  expressions  for  mass  and  heat  transfer  flux  are  discussed 
by  Lees  [39],  Bird,  Stewart  and  Lightfoot  [40],  and  Williams  [25]. 
The  mass  transfer  flux  contains  concentration,  thermal,  pressure  and 
forced  diffusion  terms.  Forced  diffusion  occurs  only  for  charged 
species  in  the  presence  of  an  electrical  field  and  does  not  have  to  be 
considered  [40].  Pressure  diffusion  is  usually  small  for  typical  flow 
conditions  [40],  and  is  absent  in  the  constant-pressure  field  of  a 
spherically  symmetric  evaporating  drop.  Thermal  diffusion  is  more 
important;  however,  the  contribution  of  the  thermal-diffusion  term  is 
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still  generally  less  than  10%,  when  compared  with  concentration 
diffusion  [39] ,  This  leaves  concentration  diffusion  as  the  major 
driving  potential  of  the  mass  flux.  Rather  than  solve  the  multi¬ 
component  diffusion  equations,  most  analyses  assume  an  effective 
binary  diffusion  expression — Fick's  Law  [25].  This  approximation  is 
adequate  when  the  binary  diffusivity  of  all  species  is  the  the  same, 
or  when  all  species,  but  one,  are  moving  at  the  same  velocity  [40]. 
For  drop  evaporation  where  the  insolubility  condition  is  satisfied, 
only  the  fuel  has  a  non-zero  velocity  and  the  condition  is  met  even 
for  fuels  with  large  molecular  weights  in  comparison  with  the  ambient 
gas. 

The  energy  flux  contains  the  conductive  flux  as  well  as 
contributions  from  interdiffusing  species  and  the  Dufour  effect  [40], 
The  Dufour  effect  is  related  to  thermodiffusion  and  is  small  even  when 
thermodiffusion  must  be  considered  [25].  The  interdiffusion  term  is 
generally  considered,  although  the  equations  assume  a  particularly 
simple  form  when  the  Lewis  number  is  unity  and  this  approximation  is 
frequently  made  [25,39],  Even  when  the  Lewis  number  is  not  unity,  if 
only  the  fuel  has  a  net  mass  flux,  only  energy  transported  by  fuel 
transport  must  be  considered. 

When  a  drop  is  surrounded  by  a  hot  oxidizing  medium,  it  can 
ignite,  giving  rise  to  a  reaction  zone  in  its  immediate  vicinity.  For 
fast  rates  of  reaction,  the  flame  zone  is  very  thin  and  surrounds  the 
drop.  In  this  case,  the  process  is  diffusion  controlled  and  bears 
many  similarities  to  drop  evaporation  without  combustion.  Sew 


a  i  icu.i  cieb  isvoivi  uore  cxcrcme  property  variations,  and  the  need  to 
determine  the  point  of  ignition.  Criteria  are  also  needed  to 
determine  when  the  flame  zone  is  thin  enough  for  the  diffusion  flame 
approximation  to  be  valid  and  when  a  stable  flame  is  established 
arourd  the  drop. 

Many  investigators  [19-23,34,41-47]  have  been  able  to  correlate 
measurements  of  droplet  gasification  rates  with  diffusion  flame 
models.  This  suggests  that  either  the  gasification  rate  itself  is 
relatively  insensitive  to  the  details  of  the  flame  structure,  or  that 
property  uncertainties  are  masking  significant  effects.  (Predicted 
burning  rates  can  be  modified  within  an  order  of  magnitude  while  still 
making  reasonable  choices  for  "average"  properties;  even  variable 
property  models  provide  considerable  latitude  due  to  the  uncertainties 
of  concentration  effects  and  convection  correlations.) 

More  sophisticated  droplet  combustion  models  allowing  for 
variable  property  effects  and  employing  the  diffusion-flame 
approximation  may  be  found  in  the  literature.  Law  and  Law  [48] 
present  a  rather  complete  version  for  a  pure  fuel  at  its  wet-bulb 
state,  which  allows  for  both  concentration  and  temperature  effects. 
These  authors  also  review  a  number  of  the  common  variable-property 
models  which  only  allow  for  variations  with  temperature. 

1.3.2  Combustion  of  Carbon 

Carbon-slurry  combustion  is  complex  since  the  combustion  of 
gaseous  or  liquid  fuels  is  superimposed  on  the  heterogeneous  oxidation 
of  carbon  particles.  To  gain  some  insight  into  the  carbon-agglomerate 
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combustion  process  of  a  combusting-slurry  droplet,  a  review  will  be 
made  of  the  combustion  of  coal  and  carbon  particles,  carbon — gas-phase 
reaction  mechanisms,  and  carbon  structure. 

1.3. 2.1  Combustion  of  Particles 

While  in  comparison  with  carbon  combustion,  the  presence  of 
ash  and  volatile  material  in  coal  modifies  the  combustion  process, 
past  observations  of  coal-particle  combustion  afford  useful  background 
information  for  the  present  investigation.  Therefore,  earlier  studies 
of  both  carbon  and  coal  particles  combustion  are  discussed  in  this 
section. 


1.3. 2. 1.1  Suspended  Coal  Particles 

There  have  been  surprisingly  few  studies  of  single  coal- 
particle  combustion.  Important  among  the  earlier  studies  was  the  work 
of  Tu,  Davis  and  Hottel  [49],  where  carbon  spheres  approximately  25  mm 
in  diameter  were  burned  within  a  furnace.  Burning  rates  and  surface 
temperatures  were  measured  for  various  furnace  temperatures,  ambient 
gas  compositions  and  gas  velocities.  A  low-temperature  regime  where 
kinetics  dominates,  and  a  high-temperature  region  where  diffusion 
dominates  were  delineated  and  interpreted  using  a  film  theory  model. 
Later,  this  data  was  reinterpreted  in  terms  of  a  combined  absorption, 
desorption  and  diffusion  mechanism  [50]. 

Following  this  initial  work  there  have  been  a  number  of  studies 
of  the  combustion  characteristics  of  carbon,  coke  and  char  [7,51-60]; 
reviews  of  this  work  may  be  found  in  References  7,  8  and  10.  This 
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area  has  received  attention  since  carbon  and  chars  are  reasonably 
well-defined  materials,  without  the  complication  of  a  devolatilization 
process.  The  chars  most  often  have  been  formed  by  slow  heating  in 
order  to  preserve  the  structure  of  the  material.  For  overall  furnace 
analysis,  the  carbon-combustion  process  is  also  much  longer  than  major 
devolatilization  processes,  making  this  aspect  more  critical  in  the 
determination  of  performance  and  heat  loading.  Questions  remain, 
however,  concerning  whether  these  structures  are  at  all  typical  of 
materials  subjected  to  rapid  devolatilization  in  a  furnace. 

Many  of  these  studies  have  involved  determination  of  the  rate  of 
reaction  by  continuously  weighing  a  particle  supported  within  a 
furnace  [49,51,54-60].  The  structural  characteristics  of  the  solid 
have  also  received  some  attention.  Orenbakh  [55,56]  employed  mercury 
porousimetry  to  estimate  pore-size  distributions  during  burnout; 
others  have  determined  density  distributions  by  progressively  grinding 
away  the  surface  and  reweighing  the  sample  [37].  Kurylko  and 
Essenhigh  [58,60]  obtained  information  on  the  temperature  distribution 
within  the  particle  using  embedded  thermocouples.  By  combining 
photographic  size  determinations  with  continuous  weighing,  it  also  has 
been  possible  to  determine  the  average  density  variation  of  the 
particle  during  burnout  [57,58,60]. 

The  gas  phase  around  a  burning  carbon  particle  has  nor.  received 
much  attention.  Kurylko  [58]  made  qualitative  measurements  around  a 
sphere.  Quantitative  measurements  were  completed  in  the  boundary 
layer  of  a  carbon  cylinder  burning  in  oxygen  by  Wicke  et  al.  [52,53]. 
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In  general,  however,  the  role  of  carbon  monoxide  combustion  in  the  gas 
phase  around  a  burning  carbon  particle  is  not  well  understood, 
particularly  in  the  presence  of  hydrocarbon  compounds  and  their 
combustion  products  [7].  This  uncertainty  in  the  ratio  of  CO2  and  CO 
production  at  the  solid  surface  has  led  to  considerable  vagueness  in 
the  interpretation  of  experimental  results,  even  at  the  diffusional 
limit  [7 ] . 

Fewer  studies  have  been  made  of  captive  coal  particles  containing 
volatiles.  Essenhigh  and  Thring  [61]  studied  the  combustion  of 
individual  captive  coal  particles  having  sizes  greater  than  1000  pm. 
Flame  observations  indicated  evolution,  ignition  and  gas-phase 
combustion  of  volatiles  as  a  first  step  in  the  process.  The  last 
stage  involved  collapse  of  the  gas-phase  flame,  followed  by  combustion 
of  the  solid  residue.  The  times  of  these  two  stages  were  proportional 
to  the  square  of  the  initial  particle  diameter,  with  the  first  stage 
approximately  an  order  of  magnitude  shorter  than  the  second.  Later 
work  by  Essenhigh  and  coworkers  [62,63]  considered  the  constants  of 
proportionality  for  the  time  of  these  two  processes  for  various  coal 
ranks  and  oxygen  concentrations.  Measurements  were  also  made  of  the 
influence  of  swelling  and  shape  on  the  combustion  rate  of  single 
particles  [ 14 J . 

Anson  et  al,  [65]  considered  supported  coal  particles  heated 
radiatively  with  projector  lamps  in  a  cold-air  environment.  Particle 
size  was  on  the  order  of  300  pm,  with  both  low-  and  medium-rank  coals 
being  considered.  The  expulsion  of  volatiles  exposed  to  high 


radiation  fluxes  caused  swelling  and  the  formation  of  hollow  spheres, 
most  often  with  a  lacy  internal  structure.  The  authors  argue  that 
combustion  proceeds  externally  and  internally  at  constant  diameter 
until  at  an  advanced  stage  the  particles  fragment.  The  ash  was  not 
observed  to  participate  in  the  process,  or  to  wet  the  surface,  thereby 
reducing  oxidizer  access  to  the  combustible  material.  The' findings  of 
this  investigation  emphasize  the  need  for  considering  materials 
prepared  at  heating  rates  typical  of  operating  conditions  if  the 
results  are  to  be  relevant  to  actual  combustion  systems. 

1.3. 2. 1.2  Free  Coal  Particles 

Following  its  introduction  by  Orning  [66],  the  gas  swept 
drop  tube,  or  one-dimensional  laminar-suspension  flow  furnace,  has 
been  a  major  tool  in  research  on  the  combustion  of  small  coal 
particles.  With  the  exception  of  a  study  by  Saji  [67],  who  considered 
some  of  the  overall  characteristics  of  pulverized  coal  flames,  the 
bulk  of  this  work  has  dealt  with  devolatilized  materials,  or  materials 
of  low  volatility.  Early  work,  represented  by  References  67-71, 
employed  pure  carbons  to  eliminate  complications  due  to  the  presence 
of  both  volatiles  and  ash.  Reviews  of  much  of  this  work  may  be  found 
in  References  7,  8  and  10.  More  recent  use  of  this  technique  by 
workers  at  BCURA  and  CSIRO  [72-77]  has  primarily  considered  prepared 
chars  of  relatively  low  volatile  content. 

The  technique  has  been  to  introduce  the  particles  into  a 
preheated  stream  of  gas,  and  then  allow  them  to  flow  through  a  furnace 
at  low  loading  so  that  temperatures  (except  particle  temperature)  are 
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primarily  fixed  by  the  furnace  temperature.  Since  the  characteristics 
of  the  particle  are  influenced  by  its  heating  rate  during 
volatilization,  in  the  most  recent  work  [72-77],  chars  are  prepared  in 
a  similar  furnace  arrangement  to  achieve  initial  heating  rates 
characteristic  of  furnace  applications,  (10^  K/s). 

With  the  laminar-suspension  flow-furnace  technique,  burning  rate 
measurements  are  made  by  sampling  the  particles  after  various  times  of 
exposure  in  the  furnace.  Major  variables  considered  have  been 
particle  size,  feed  rate,  coal  type,  gas  velocity,  gas  composition 
(oxygen  concentration-  etc.)  and  furnace  temperature.  Burning  rates 
are  measured  both  gravimetrically  and  by  the  appearance  of  CO  and  CO2 
in  the  gas  stream.  Particle-density  variation  has  been  determined  by 
weighing  and  sizing  particles  at  various  sampling  stations. 

Using  a  similar  apparatus  in  tesis  of  bituminous  coal,  Anson  et 
al.  [65]  determined  surface  area  through  BET  techniques.  Guerin  et 
al.  [78]  employed  similar  methods  to  infer  the  characteristics  of  the 
pore  structure  of  a  pulverized  coal  char  contained  in  a  combustion 
boat  within  a  furnace  (non-flowing). 

The  one-dimensional  furnace  has  successfully  been  used  for 
studying  the  entire  combustion  process  of  coal  by  Howard  [79]  and 
Howard  and  Essenhigh  [30,81].  Samples  withdrawn  from  the  furnace  were 
analyzed  in  order  to  determine  the  variation  of  volatile  matter,  ash 
and  fixed  carbon  content  of  the  particles.  Temperature  and  gas- 
composition  measurements  were  also  made  in  the  furnace.  In  order  to 
simulate  combustion  at  elevated  pressures,  where  heating  rates  of  the 
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particles  are  even  higher  than  at  atmospheric  pressure,  a  few  studies 
have  been  conducted  in  shock  tubes  [82]. 

There  is  general  recognition  that  pulverized-coal  combustion  is 
so  complicated  that  there  is  little  likelihood  of  formulating  valid 
and  accurate  analytical  combustion  models  unless  they  are  based  on 
thorough  experimental  delineation  of  key  processes.  For  that  reason, 
current  coal-particle  combustion  models  are  oversimplified,  simply 
because  the  available  experimental  data  cannot  support  a  higher  degree 
of  sophistication.  Of  course,  the  degree  of  oversimplification 
depends  upon  the  user's  viewpoint;  if  the  burning  rates  are  of  primary 
interest,  the  present  models  give  credible  predictions  since  that  is 
the  principal  kind  of  experimental  data  upon  which  they  are  based. 

The  most  extensive  models  consider  devolatized  materials.  All 
employ  conventional  spherically  symmetric  film  theory  for  mass 
transfer  to  a  solid  sphere,  with  consideration  of-  some  type  of 
reaction  at  the  surface  [10,61-78];  although  just  what  is  diffusing  is 
not  clear  in  all  cases  due  to  competing  processes  of  CO  and  CO2 
formation.  Some  models  allow  for  diffusion  through  an  ash  layer 
[10,69];  in  some  cases,  diffusion  in  the  carbon  pore  structure  is 
considered  [8,57,58,60],  and  gas-phase  reaction  has  been  treated  for 
cylindrical  geometries  [52].  Most  often,  however,  gas-pbase  reaction 
is  neglected,  the  particle  is  assumed  to  burn  over  its  outer  surface, 
and  the  actual  products  of  combustion  are  assumed  [10,72-78]. 
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1.3. 2. 1.3  Theoretical  Models 

In  most  theoretical  studies  of  burning  carbon  particles, 
only  radial-diffusive  transport  in  the  gas  phase  and  the  heterogeneous 
reactions  are  considered  to  influence  the  overall  reaction  rate.  In 
early  work,  Spalding  [83]  and  Avedesian.  and  Davidson  [84]  included 
gas-phase  chemistry  through  a  flame  sheet  which  stands  off'  the  surface 
of  the  particle  where  CO  was  oxidized.  However,  these  analyses  are 
limited  by  the  assumption  of  negligible  radial  convection; 
furthermore,  the  heterogeneous  reactions  are  assumed  to  be  fast  enough 
so  that  the  carbon  consumption  is  diffusion  dominated.  The 
theoretical  work  of  Ubhayakar  and  Williams  [85],  which  included  the 
influence  of  radial  convection  in  the  gas  phase,  supports  their 
experimental  results-  In  this  case  the  authors  assumed  that  the  gas- 
phase  chemistry  is  frozen  and  that  only  the  carbon-oxygen  reaction  is 
operative. 

Studies  by  Amundson  and  coworkers  [86,87]  considered  burning  of 
spherical  carbon  particles  surrounded  by  finite,  quiescent  gas  films. 
Cnram  and  Amundson  [86]  also  provide  an  extensive  review  of  single¬ 
particle  carbon  combustion  studies,  including  information  on  the 
various  fluid-mechanical  and  chemical  models  used  in  each.  The 
emphasis  in  both  Caram  and  Amundson  [86]  and  Mon  and  Amundson  [87]  is 
on  the  gas  phase  alone. 

Libby  and  Blake  [88]  studied  the  particle-life  history  of  a 
single-carbon  particle  in  a  hot  oxidizing  ambient  environment.  Their 
theoretical  model  was  based  on  a  unitora-temperature ,  constant- 
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density,  spherical  particle  of  carbon  surrounded  by  a  quiescent  gas 
phase  having  constant  transport  and  thermophysical  properties.  The 
direct  oxidation  of  carbon  with  oxygen  and  the  indirect  oxidation  via 
carbon  dioxide  were  taken  into  account.  While  both  mechanises 
influenced  the  dynamic  behavior  of  carbon  particles,  mass  loss  by 
oxygen  dominated  the  process.  Theoretical  results  of  Libby  and  Blake 
[88J  indicate  that  particle  consumption  involves  complex  interactions 
of  convective,  diffusional  and  thermochemical  effects.  In  a  later 
paper,  Libby  and  Blake  [89]  extended  their  model  to  include  the 
effects  of  water  vapor  and  hydrogen.  Their  prediction  of  rates  of 
particle  mass  loss  for  specified  ambient  conditions  and  particle 
temperature  and  size  suggest  that  in  the  high  temperature  or  the 
diffusion  limit,  the  influence  of  water  vapor  and  oxygen  on  particle- 
mass  loss  are  approximately  equivalent.  However,  at  the  kinetic 
limit,  oxygen  yields  higher  rates  of  mass  loss. 

Combustion  of  carbon  particles  in  an  oxidizing  atmosphere  was 
also  studied  by  Matalon  [90].  In  this  investigation  the  author 
describes  the  departure  from  an  equilibrium  gas-phase  behavior  for  a 
pure  carbon  particle  immersed  in  an  oxidizing  atmosphere  over  a  wide 
range  of  particle  sizes  and  pressures.  Matalon  [90]  determined  that 
the  oxidation  of  the  carbon  particle  depended  on  both  direct  oxidation 
by  oxygen  and  indirect  oxidation  by  carbon  dioxide.  The  latter 
mechanism  was  modeled  as  a  CO  flame  sheet  surrounding  the  particle. 


1.3. 2. 2  Gas-Solid  Reaction  Mechanisms 

From  the  material  reviewed  in  the  previous  section,  it  is 
apparent  that  accurate  predictions  of  the  reaction  rate  of  carbon  in  a 
hot  oxidizing  environment  requires  a  thorough  understanding  of  the 
gas-solid  reaction  mechanisms.  Therefore,  this  section  provides  a 
review  of  the  mechanisms  which  are  pertinent  for  carbon  reaction. 

It  is  generally  accepted  that  there  are  three  principal  global 
reactions  of  carbon  which  are  of  major  interest  in  the  carbon-particle 
combustion  process: 

cs  +  °2  C02  ,  CO 

Cs  +  C02  *  2  CO  (1.2) 
Cs  +  H2°  ■*"  co  +  «2 

The  reactions  are,  of  course,  much  more  complex  than  is  indicated  by 
the  simple  schemes  shown;  a  common  feature  is  that  they  are  all 
accompanied  by  the  formation  of  surface  complexes.  Also,  the  behavior 
of  these  surface  complexes  is  determined,  in  part,  by  the  properties 
of  the  carbon  surface.  This  reaction  surface  cannot  be  defined 
adequately,  and  according  to  the  circumstances  in  which  the  reaction 
takes  place,  both  the  internal  and  external  surface  of  the  carbon 
particle  may  be  involved  in  the  reaction.  In  the  following,  each  of 
the  reactions  described  in  Equation  (1.2)  will  be  discussed. 


The  chemical  reaction  of  oxygen  with  solid  carbon  is  generally 
agreed  to  occur  by  the  processes  of  adsorption  and  desorption  [7]. 
Although  there  is  some  difference  of  opinion  concerning  the  detailed 
steps  by  which  these  processes  occur,  a  relatively  simple  adsorption- 
desorption  model  was  developed  by  Lewis  and  Simons  [91]  as  follows: 


Cs  +  °2  -  C(02) 

Cs  +  C(02)  +  2  CO 


(1.3) 


where  CCC^)  is  an  02  molecule  adsorbed  on  the  carbon  surface.  Other 
investigators  [8,9,72]  prefer  a  more  direct  oxidation  of  carbon  by 
oxygen  mechanism  as  expressed  in  Equation  (1.2).  However,  there  is  a 
great  deal  of  controversy  concerning  the  interpretation  of  existing 
data  which  has  led  to  a  wide  range  of  activation  energy  and  reaction 
orders  for  even  this  well-documented  reaction  [8,9,72].  For  example, 
data  on  pulverized  coal  and  carbon  combustion  for  temperatures  greater 
than  1200  K  given  by  Field  [72,73],  Mulcahy  and  Smith  [7],  Khitrin  and 
Golovina  [92],  and  Smith  [74,75]  indicate  the  reaction  to  be  first 
order  in  the  partial  pressure  of  oxygen  with  an  apparent  activation 
energy  of  20  to  40  kcal/mole.  Such  behavior  is  in  accord  with  the 
accepted  heuristic  consideration,  although  the  observed  activation 


energy  appears  to  be  large. 
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The  principal  mechanistic  studies  of  the  carbon-carbon  dioxide 
reaction  have  been  reported  by  Gadsby  et  al.  [93],  Ergun  et  al. 
[94-96],  Walker  et  al.  [97,98],  von  Fredersdorff  and  Elliott  [99], 
Lewis  et  al.  [100],  and  Reit  [101].  In  each  case,  the  authors  suggest 
the  following  oxygen-exchange  mechanism: 


Cs  +  C02  C  CO  +  CCO) 
C(O)  t  CO 


(1.4) 


This  mechanism  is  the  most  widely  accepted  explanation  of  the  course 
of  the  C02  reaction  below  or  near  atmospheric  pressure  [99];  however, 
the  other  investigators  [101,102]  support  a  mechanism  with 
unidirectional  oxygen  exchange  at  the  surface. 

Gadsby  et  al.  [102]  supported  the  second  carbon  dioxide  mechanism 
on  the  basis  of  charcoal  chemisorption  experiments  at  750°C  and  about 
300  to  400  mm  Hg  total  pressure,  which  indicated  that  variations  in 
CO 2  partial  pressures  corresponded  more  closely  to  variations  in  CO 
partial  pressure  than  to  changes  in  amount  of  oxygen  complex  absorbed. 
In  addition,  the  quantity  of  oxygen  complex  remained  relatively 
constant  with  variation  in  CO  partial  pressure  at  constant  C02  partial 
pressure.  This  was  taken  as  evidence  that  the  rate  of  inhibition  must 
have  been  set  by  CO  chemisorption,  not  by  the  reverse  oxygen  exchange 
reaction  of  Equation  (1.4). 
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In  addition  to  the  uncertainty  about  which  mechanise  correctly 
predicts  the  oxidation  of  carbon  by  carbon  dioxide,  once  a  mechanism 
is  chosen  there  is  a  wide  range  of  reported  activation  energies  and 
reaction  orders  [99].  As  in  the  case  of  oxygen  reaction,  the 
activation  energy  varies  by  a  factor  of  two,  from  30  to  60  kcal/mole 
[100-L02J. 

The  kinetics  of  the  carbon-steam  reaction  is  in  many  respects 
analogous  to  the  CO2  reaction.  A  unified  theory  of  the  course  of  gas 
reactions  of  carbon,  where  the  reactants  are  oxidizing,  holds  that 
their  mechanisms  or  elementary  steps  should  be  similar.  However,  the 
C-CO2  system  comprises  only  two  constituents — carbon  and 
oxygen — allowing  less  uncertainty  about  the  choice  of  its  limited 
number  of  possible  mechanisms.  The  carbon-steam  system,  on  the  other 
hand,  poses  the  added  complexity  contributed  through  possible  roles  of 
hydrogen  and/or  other  intermediate  dissociation  products  of  water 
vapor.  Furthermore,  the  carbon  dioxide  reaction  with  carbon  is 
occurring  simultaneously. 

The  principal  mechanistic  studies  of  the  carbon-steam  reaction 
have  been  reported  by  Gadsby  et  al.  [102],  Johnstone  et  «1.  [103], 
von  Fredersdorff  and  Elliott  [99],  Long  and  Sykes  [1Q4],  Strickland- 
Constable  [105],  Wicke  and  Rossberg  [106],  Binford  and  Eyring  [107] 
and  Eurgun  [95]-  In  each  case,  the  authors  recommend  the  following 
oxygen  exchange  mechanism: 


(1.5) 


Cs  +  H20  -  H2  tC(O) 

C(0)  ^  CO 


There  seems  to  be  general  agreement  that  the  primary  products  of  the 
steam  reaction  with  relatively  pure  carbons  at  atmospheric  or  lower 
pressures  are  carbon  monoxide  and  hydrogen,  whereas  C02  is  a  secondary 
product  arising  through  the  water-gas  shift  reaction  which  is  known  to 
be  catalyzed  by  carbon  surfaces  [93,108] . 

Similar  to  the  carbon-carbon  dioxide  mechanisms,  there  are 
several  other  proposed  carbon-steam  mechanisms  and  a  wide  range  of 
reported  activation  energies  and  reaction  orders  [99].  Activation 
energy  for  the  carbon-steam  reaction  has  been  reported  from  10  to  65 
kcal/mole  [99].  In  addition,  the  reaction  rate  of  water  vapor  and 
carbon  is  complicated  by  oxygen  and  carbon  dioxide  all  competing  for 
the  same  active  sites. 

In  addition  to  the  reaction  mechanism  which  investigators  have 
suggested  for  the  reaction  of  carbon  or  coal,  oth  investigators 
[109-121]  have  determined  separate  reaction  rat.  aechanisms  for 
graphite,  soot  and  carbon  black.  Nagle  and  Strickland-coi  stable  [109] 
formulated  a  semi-empirical  formula  based  on  reaction  of  oxygen  with 
pyrolytic  graphite.  Rosner  and  Allendorf  [110]  nave  determined  that 
for  pyrolytic  graphite  monatomic-oxygen  atoms  tend  to  be  more  reactive 
than  diatomic  oxygen  .  Wright  [111]  has  also  correlated  the  reaction 
rate  of  soot  and  pyrolytic  graphite  with  monatomic  oxygen  as  the 
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oxidant.  Ottertein  and  Bonnetain  [112]  and  Appleton  et  al.  [115,116] 
compared  the  rates  of  soot  oxidation  In  flames  [117-119],  pyrolytic 
graphite  oxidation  using  oxygen  [110,113]  and  carbon  black  oxidation 
in  a  shock  tube  [120],  with  the  rates  predicted  by  the  Nagle  and 
Str ickland-Constable  formula  [109].  They  concluded  that  the  magnitude 
of  most  of  the  soot  and  pyrolytic-graphite  oxidation  rates  can  be 
adequately  correlated  by  the  Nagle  and  Strickland— Constable  formula. 

Although  the  Nagle  and  Strickland-Constable  formula  seems 
adequate  for  predicting  oxidation  rates  of  soot  and  pyrolytic  graphite 
by  oxygen,  thare  is  a  lack  of  evidence  that  this  formula  will 
correctly  predict  soot  oxidation  rates  in  flames  where  there  are  more 
than  one  potential  oxidant.  In  fuel-rich  flames  the  02  concentration 
can  be  comparable  to  the  0  concentration  and  also  can  be  much  less 
than  the  OH  concentration.  The  reaction  probability  of  soot  with  0 
and  OH  can  be  expected  to  be  higher  than  with  02 .  Neoh,  Howard  and 
Sarofim  [121]  determined  that  at  high  temperature  and  fuel-rich 
conditions,  the  OH  radical  was  the  principal  oxidant,  with  oxygen 
dominating  only  in  the  lower-temperature  and  fuel— lean  conditions. 

1.3. 2. 3  Carbon  Structure 

The  chemistry  and  fluid  mechanics  of  the  combustion  and 
gasification  of  porous  carbon  char  have  been  studied  by  numerous 
researchers  and  summarized  in  two  extensive  review  articles,  one  by 
Walker  et  al.  [8j  one  by  and  Mulcahy  and  Smith  [7].  It  is  known  that 
the  gasification  process  can  depend  upon  the  gaseous  environment  and 
porous  structure,  and  is  either  diffusion  limited  or  kineticaliy 


37 


limited.  The  diffusion  process  may  be  controlled  by  Knudsen  diffusion 
within  the  pores,  or  by  continuum  diffusion  either  within  the  pores  or 
in  the  boundary  layer  surrounding  the  particle.  In  addition,  the 
heterogeneous  rate  constant  may  be  limited  by  either  the  adsorption  or 
the  desorption  process.  The  complications  introduced  by  the  presence 
of  many  limiting  processes  are  compounded  by  lack  of'  a  complete 
description  of  the  porous  structure.  Recognizing  these  difficulties, 
transport  processes  have  been  modeled  in  highly  idealized  situations. 
Thiele  [122]  described  the  diffusion  and  heterogeneous  chemistry 
occurring  within  a  single  isolated  cylindrical  pore  and  determined 
that  the  diffusion  limited  gasification  process  scales  as  the  square 
root  of  the  kinetic  rate.  Wheeler  [123]  applied  the  Thiele  solution 
to  a  pore  structure  which  consisted  of  a  large  number  of  identical 
parallel  pores.  He  demonstrated  that  the  gasification  process,  when 
limited  by  Knudsen  diffusion,  scaled  as  the  porosity  over  the  square 
root  of  the  kinetic  rate. 

While  the  works  of  Thiele  and  Wheeler  have  yielded  a  better 
understanding  of  the  gasification  process,  the  models  are  rigorously 
valid  only  in  extreme  limits,  since  they  do  not  include  a  realistic 
pore  size  distribution  or  a  description  of  pore  branching  and  pore 
combination.  These  effects  have  been  incorporated  by  Hashimoto  and 
Silverston  [124,125]  using  an  integral  technique  with  adjustable 
empirical  parameters  to  close  the  mathematical  system.  A  qualitative 
description  of  the  pore-distribution  function,  pore  branching  and  pore 
coabi.nati.on  was  not  developed,  and  the  adjustable  parameters  handled 
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this  part  of  the  model.  A  model  proposed  by  Kriegbaum  and  Laurenaeau 
[126]  considers  the  gasification  of  pulverized  or  crushed  coal  char 
particles.  This  model  characterizes  the  porous  solid  reactant  as  a 
system  of  individual  pores.  A  pore  is  assumed  to  consist  of  a  tubular 
capillary  bounded  radially  by  a  nonporous  carbon  sheath.  This  carbon 
sheath  is  assumed  to  be  bounded  radially  by  nonporous  ash.  The  ash 
serves  to  limit  the  radial  growth  of  the  pore.  The  model  combines 
Knudsen  capillary-flow  theory  with  standard  molecular-diffusion 
theory.  The  resulting  diffusion  coefficient  accounts  for  effects  of 
pore  growth  on  diffusion  since  it  is  a  function  of  pore  radius.  The 
model  also  uses  Langmuir-Einshelwood  kinetics  to  model  the  surface 
reaction.  Analysis  of  the  simultaneous  mass  transport  plus  reaction 
within  the  individual  pores  yields  burnoff  data  for  the  entire 
particle.  The  model,  however,  is  incapable  of  predicting  pore 
structure  development  with  burnoff,  since  pore  interactions  are  not 
considered. 

Beshty  [127]  developed  a  mathematical  model  to  analyze  the 
internal  combustion  of  a  porous  carbon  sphere  with  uniform  spherical 
pores  arranged  on  a  simple  cubic  lattice  extending  through  the 
particle  volume.  Beshty  assumed  that  as  the  reaction  proceeded  the 
pores  enlarged  and  eventually  merged,  causing  the  partial  collapse  of 
the  particle.  The  model  was  used  to  examine  the  effect  of  particle 
si2e,  temperature  and  extent  of  burnoff  on  the  rate  of  combustion  and 
the  change  in  internal  specific  surface  area.  Three  combustion 
regimes  were  predicted  by  the  model:  chemically  controlled  uniform 
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burning,  joint  chemically  and  pore-diffusion  controlled  combustion, 
and  a  pore-diffusion  region  leading  to  nonuniform  burning.  As  the 
temperature  of  the  particle  and  particle  size  increases,  the  effects 
of  pore  diffusion  becomes  more  significant.  Although  temperature  and 
particle  size  determine  the  initial  combustion  region,  as  burnoff 
progresses,  transition  to  the  chemically  controlled  region  gradually 
occurs.  3eshty  [127]  determined  that  this  is  a  result  of  the 
diminishing  effect  of  pore  diffusion  as  combustion  progresses  and 
pores  develop. 

Simons  et  al.  [92,128-130]  developed  a  model  to  describe  the  mass 
transport  and  heterogeneous  chemistry  which  occurs  during  the 
combustion  and/or  gasification  of  porous  coal  char.  The  authors 
develop  a  semi-empirical  model  to  describe  the  porous  coal  char 
structure.  The  pores  were  assumed  to  be  cylindrical  tubes  connnected 
to  other  pores  of  the  same  shape  to  form  pore-true  structures.  Data 
suggest  that  the  distribution  of  pore  number  density  is  inversely 
proportional  to  the  pore  radius  cubed.  Simons  and  Finson  [128] 
determined  that  the  smaller  pore  trees  were  kinetically  limited  while 
the  larger  trees  were  diffusion  limited.  This  is  the  opposite  to  that 
which  one  would  conclude  from  the  Thiele-Wheeler  [122,123]  approach, 
which  does  not  include  pore  branching  or  the  concept  of  a  pore  tree. 

Simons  [129]  developed  a  mechanistic  model  which  describes  pore 
combination  during  char  gasification  by  expressing  a  general 
population  balance  for  a  given  pore  size  which  includes  pore  growth 
and  pore  engulfment.  The  engulfment  of  the  smaller  pores  by  the 
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larger  ones  was  responsible  for  Che  relative  depletion  of  the  small 
pores  and  enhanced  growth  of  the  larger  pores.  Simons  also  concluded 
that  the  influence  of  the  parent  coal  and  effects  of  pyrolysis  and 
pore  growth  during  gasification  are  as  significant  as  the  effect  of 
pore  combination  and  must  be  considered  before  the  time  history  of  the 
pore  distribution  may  be  accurately  predicted. 

The  effects  of  Knudsen  diffusion,  continuum  diffusion,  and  both 
adsorption  and  desorption  kinetics  were  considered  by  Simons  [129]. 
He  concluded  that  each  physical  mechanism  may  limit  the  gasification 
rate  within  some  range  of  pore  size,  particle  size,  gas  temperature 
and  gas  pressure.  In  general,  however,  no  single  mechanism  controls 
the  entire  gasification  process. 

Although  there  are  several  models  which  can  predict  the  pore 
structure  of  coals  and  chars  reasonably  well,  no  models  considering 
agglomerate  pore  structure  were  found  in  the  literature. 

1.3.3  Combustion  of  Slurry  Droplets 

The  combustion  of  a  pure-component  fuel  droplet  can  be 
considered  to  be  reasonably  well  understood  [11],  and  its  burning  rate 
can  be  predicted  quite  accurately  by  applying  the  classical  d  -Law 
[131].  This  law  assumes  a  constant  droplet  temperature  when  the 
process  of  droplet  heating  is  taken  into  account  [131,132].  A  short 
period  spanning  about  10Z  of  the  initial  droplet  lifetime,  at  low 
pressures,  exists  during  which  the  droplet  is  actively  heated  up  and 
the  burning  rate  is  very  slow.  It  was  further  shown  that  the  bulk 
droplet  combustion  characteristics  depend  only  minimally  on  the 
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detailed  description  of  the  heat  transport  mechanism  within  the 
droplet  interior  [132]. 

However,  studies  [133-135]  on  miscible  multi-component  droplets 
show  that  their  combustion  behavior  depends  critically  on  the  liquid- 
phase  heat  and  mass  transport  processes.  In  particular,  two  limiting 
cases  have  been  identified.  In  the  diffusion  limit  [133],  the 
external  gas  motion,  which  exerts  shear  stress  at  the  droplet  surface, 
produces  essentially  no  motion  within  the  droplet  interior  because  of 
the  small  droplet  size  and  high  liquid  viscosity.  In  the  convective 
limit,  the  internal  states  of  the  droplet  are  assumed  to  be  maintained 
as  perpetually  uniform  since  the  rate  of  internal  circulation,  and 
therefore  the  rate  of  convective  transport,  is  fast. 

Recently,  the  above  concepts  have  been  extended  to  systems 
involving  more  than  one  liquid  phase.  The  use  of  emulsions  of  water 
in  oil  has  attracted  considerable  attention  as  a  means  of  reducing 


emissions 

of 

particulate  solids 

and 

saving  fuel 

because  of 

the 

reduction 

in 

combustion  air  that 

can 

be  achieved 

[136-143] . 

The 

addition 

of 

small  quantities  of 

water  significantly  enhances 

the 

potential  and  intensity  of  droplet  explosion  [144]. 

In  the  combustion  of  the  carbon-oil  mixtures,  the  oil  is  much 
more  volatile  than  the  carbon  particles.  Therefore,  it  is  reasonable 
to  expect  that  the  oil  will  be  vaporized  first  and  that  a  large  carbon 
particle  agglomerate  will  be  left  behind.  This  possible  agglomeration 
implies  that  the  sizes  of  the  carbon  particles  which  actually  undergo 
combustion  in  the  form  of  agglomerates  are  determined  by  the  initial 
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concentrations  of  the  carbon  particles  and  the  initial  droplet  sizes, 
but  not  by  the  sizes  of  the  individual  carbon  particles.  Therefore, 
it  is  probably  not  necessary  to  use  extremely  fine  carbon  powders  in 
the  preparation  of  slurries  [142],  at  least  from  the  combustion 
standpoint.  These  agglomerates  will  also  affect  the  rate  of  radiation 
transfer,  which  is  an  important  heat-transfer  mechanism  within 
furnaces  and  boilers.  Also,  the  large  agglomerate  size  and  non¬ 
volatile  nature  of  the  carbon  may  lead  to  incomplete  combustion  [142]. 

On  the  other  hand,  it  has  also  been  suggested  by  Law  et  al. 
[144]  that  agglomeration  may  be  either  prevented  or  minimized  through 
the  bubbling  of  Some  volatile  component  trapped  within  the  droplet 
interior.  The  resultant  internal  pressure  buildup  will  shatter  the 
agglomerate  into  much  smaller  pieces.  This  so-called  "micro- 
explosion"  event  has  been  theoretically  predicted  by  Law  [134,135]  and 
experimentally  verified  by  Lasheras  et  al.  [145-147]  for  miscible  and 
immiscible  liquid  mixtures  consisting  of  oil,  water  and  alcohol. 

A  preliminary  experiment  was  conducted  by  Law  et  al.  [142]  at 
room  temperature  air  on  the  burning  of  a  single  coal-oil  mixture 
droplet  suspended  at  the  end  of  a  quartz  fiber.  A  mixed  behavior  was 
observed  showing  the  occasional  ejection  of  small  masses  from  the 
droplet  and  the  eventual  formation  of  an  irregularly  shaped  coal- 
particle  agglomerate.  The  agglomerate  was  formed  when  burning  was 
terminated  upon  complete  depletion  of  the  oil  component.  Addition  of 
small  quantities  of  water  intensified  the  fragmentation. 
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Law  [L48]  also  studied  free-falling  droplets  in  a  hot  environment 
simulating  the  furnace  interior.  The  use  of  free-falling  droplets 
allowed  the  agglomerates  to  be  formed  without  uncertainties  due  to  the 
presence  of  the  suspension  fiber.  It  was  determined  [148]  that  the 
coal  burns  in  the  form  of  agglomerates  instead  of  individual 
particles.  Therefore,  from  combustion  considerations  it  is 
unnecessary  to  fine-crush  the  coal  unless  the  particles  can  be  reduced 
to  micron-size.  Law  [148]  also  determined  that  attempts  to  make  the 
coal-oil  mixture  more  volatile  in  the  hope  of  facilitating  combustion 
nay  in  fact  prolong  the  total  burning  time,  since  a  more  volatile  oil 
produces  delays  in  agglomeration,  heat-up  and  ignition  of  the 
agglomerate. 

1.4  Specific  Objectives  of  the  Study 

As  discussed  in  the  last  section,  many  investigators  have  studied 
the  combustion  process  of  both  single-suspended  and  free-flowing 
carbon  particles.  Combustion  of  coal-slurry  droplets  has  also  been 
observed.  However,  no  method  has  been  demonstrated  which  can  provide 
predictions  of  carbon-slurry  combustion  properties  under  typical 
combustion-chamber  conditions.  The  main  objective  of  this 
investigation  was  to  develop  a  comprehensive  theory  for  the  combustion 
of  carbon-slurry  drops  and  to  evaluate  the  theory  by  comparison  with 
measurements. 

The  experimental  portion  of  this  investigation  involved  observing 
the  combustion  of  slurry  drops  supported  in  a  turbulent-diffusion 
flame  and  free  drops  in  the  post-flame  region  of  a  flat-flame  burner. 
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These  test  conditions  were  examined  in  order  to  simulate  the  local 
environment  of  drops  within  a  combustion  chamber.  In  order  to  study 
the  nature  of  agglomeration  of  the  carbon  particles,  slurry  droplets 
were  extinguished  at  various  stages  in  its  combustion  process.  The 
resulting  agglomerate  structure  was  studied  using  a  scanning  electron 
microscope. 

The  specific  objectives  of  each  experimental  test  apparatus  is 
discussed  in  the  following  sections. 

1.4.1  Combustion  of  Large,  Suspended  Slurry  Droplets  in  ja  Turbulent- 

Diffusion  Flame 

The  specific  objectives  of  the  combustion  tests  of  large, 
suspended  slurry  droplets  in  a  turbulent-diffusion  flame  can  be 
summarized  as  follows: 

1.  Identify  the  conditions  under  which  a  slurry  droplet  will 
ignite  and  completely  burn. 

2.  Determine  if  carbon  particles  are  ejected  from  the  carbon- 
slurry  droplet  during  the  liquid-evaporation  stage  as  has 
been  observed  for  coal-slurry-droplet  combustion. 

3.  Determine  the  effect  of  a  catalyst  on  the  carbon  slurry 
droplet  combustion  process. 

4.  Examine  the  structure  of  the  carbon  residue  at  various  times 
during  the  combustion  process. 

5.  Determine  the  different  stages  of  the  combustion  process  of 
a  slurry  droplet  and  their  relative  duration. 
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6.  Develop  and  evaluate  a  liquid-evaporation  model  to  predict 
the  heat-up  and  evaporation  of  the  liquid  phase  of  a  slurry 
droplet.  The  model  should  be  able  to  accurately  predict  the 
droplet  diameter  and  temperature  during  the  liquid 
evaporation  process. 

7.  Develop  and  evaluate  a  carbon  agglomerate  combustion  model 
to  predict  the  heat-up  and  combustion  of  the  particle.  The 
model  should  be  able  to  accurately  predict  the  particle 
diameter  and  temperature  during  the  combustion  process. 

8.  Determine  which  gas  species  most  readily  reacts  with  carbon. 

1.4.2  Combustion  of  Small ,  Free-Falling  Carbon  Particles  in  a  Laminar— 

Premixed  Flame 

The  specific  objectives  of  the  combustion  tests  of  small,  free- 
falling  carbon  particles  in  a  laminar- premixed  flat-flame  burner  can 
be  summarized  as  follows: 

1.  Study  the  combustion  process  of  carbon  particles  typical  of 
the  size  that  would  be  encountered  in  a  combustion  chamber. 

2.  Develop  an  experimental  apparatus  which  eliminates  the 
effect  of  the  support  probe  from  the  carbon-particle 
combustion  process. 

3.  Study  the  effect  of  different  equivalence  ratios  of  the 
flame  gas  at  constant  temperature,  and  different  temperature 
levels  for  a  given  equivalence  ratio. 

4.  Identify  the  conditions  which  lead  to  the  combustion  of  the 
carbon  particle  by  either  diffusion-controlled  combustion  or 
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kinetically  controlled  combustion.  Most  of  the  regions  of 
the  turbulenc-dif fusion  flame  test  approached  the 
diff us tonally— controlled  limit  or  the  quenching  limit,  which 
reduced  the  sensitivity  of  the  reaction  rate  evaluation. 
Therefore,  more  kinetically  controlled  test  conditions  will 
be  studied. 

5.  Examine  the  variation  of  particle  density  during  the 
agglomerate  combustion  process,  by  measuring  particle  mass 
as  well  as  particle  diameter. 

6.  Evaluate  the  theoretical  model  for  these  test  conditions. 


CHAPTER  II 


EXPERIMENTAL  APPARATUS  AND  PROCEDURES 

2. 1  Introduction 

The  first  experimental  arrangement  involved  a  single  slurry 
droplet  suspended  from  a  thermocouple  probe  in  an  atmospheric 
pressure,  turbulent-diffusion  flame.  This  arrangement  was  considered 
since  slurry  drops  placed  at  various  locations  in  the  flame  would 
experience  conditions  similar  to  those  in  a  combustion  chamber.  By 
mounting  the  droplet  on  a  thermocouple,  the  apparatus  also  afforded  an 
opportunity  to  measure  the  temperature  of  the  carbon-slurry  droplet 
during  combustion. 

The  second  experimental  apparatus  involved  injecting  slurry 
droplets  through  a  laminar-premixed  flame.  This  arrangement  provided 
a  Deans  of  observing  the  combustion  of  very  small  slurry  droplets, 
while  maintaining  control  of  the  surroundings  of  the  drops,  and 
eliminating  any  effects  due  to  the  presence  of  a  support  probe. 

Both  apparatus  are  discussed  in  detail  in  the  remainder  of  this 
chapter. 

2.2  Turbulent  Diffusion  Flame  Apparatus 
2.2.1  Test  Apparatus 

The  test  apparatus  involved  a  modification  of  an  earlier 
arrangement  employed  for  a  study  of  a  gaseous  turbulent-diffusion 
flame  [149,150,151].  A  sketch  of  the  apparatus  appears  in  Figure  6. 
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Figure  6.  Schematic  diagram  of  the  turbulent  diffusion  flame  apparatus. 
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The  arrangement  consists  of  an  injector  flowing  propane  gas,  which  is 
burned  as  a  diffusion  flame  in  stagnant  room  air.  The  carbon 
agglomerate  was  supported  from  a  fine-wire  thermocouple  probe  at 
various  locations  in  the  flame,  photographed  with  a  motion-picture 
camera,  and  the  particle  temperature  was  recorded  with  an  analog-to- 
digital  signal  converter. 

The  gas  injector  was  positioned  near  the  bottom  of  the  test  stand 
and  oriented  vertically  upward.  The  test  stand  was  an  area  1.2  m 
square  and  3  m  high  enclosed  with  a  single  thickness  of  16  mesh 
screen.  The  test  stand  was  located  within  a  room  having  dimensions 
4mx7mx4m  high.  Combustion  products  vrere  removed  through  an 
exhaust  fan  located  near  the  ceiling  of  the  test  cell.  The  injector 
was  mounted  on  a  support  which  allowed  a  vertical  movement  of  1.2  m. 
Two  traversing  mechanisms  were  used  to  position  the  nozzle  in  the 
horizontal  plane. 

The  nozzle  was  a  full-cone  air-atomizing  injector,  with  no  swirl, 
manufactured  by  the  Spraying  Systems  Company  (model  1/4  J  2050  fluid 
nozzle  and  model  67147  air  nozzle).  The  outlet  diameter  of  the  nozzle 
was  1.194  mm. 

The  gas  flame  was  stabilized  near  the  exit  of  the  injactor  by 
employing  an  array  of  four  hydrogen  capillary  flames.  The  flame  tubes 
were  mounted  in  a  symmetrically  opposed  pattern  injecting  toward  the 
centerline  of  the  injector  passage.  The  flow  passages  for  the 
hydrogen  were  slots  0.4  mm  wide  and  2  mm  long,  with  the  long  axis 
parallel  to  the  injector  centerline.  The  lower  edges  of  the  slots 
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were  1.8  nn  above  the  face  of  the  injector.  The  exit  planes  of  the 
slots  were  6.4  nan  apart. 

The  propane  gas  flow  rate  was  metered  with  a  Matheson  model  604 
rotameter,  and  the  pressure  was  controlled  by  a  Harris  model  2515 
pressure  regulator  with  an  output  capacity  of  0.069  MPa.  The  pressure 
upstream  of  the  rotameter  was  monitored  with  a  Heise  absolute-pressure 
gauge  with  a  range  of  0.0-2.07  MPa,  and  the  downstream  pressure  was 
monitored  with  a  Heise  absolute-pressure  gauge  with  a  range  of 
0.0-0.41  MPa.  The  flow  race  of  the  hydrogen  gas  was  metered  with  a 
Matheson  model  601  rotameter,  and  the  pressure  was  controlled  with  a 
Matheson  model  1H  pressure  regulator  having  a  0.138  MPa  output 
capacity.  The  rotameters  for  both  the  propane  and  hydrogen  flows  were 
calibrated  with  a  Precision  Scientific  Company  wet-test  meter 
(233  ml/rev). 

Slurry  droplets  were  mounted  on  fine-wire  platinum/platinura-lOX 
rhodium  thermocouple  with  a  50  ym  wire  diameter.  The  thermocouple 
junction  was  encased  in  a  400  ym  bead  of  Sauereisen  cement  to  help 
support  the  drop.  The  thermocouple  was  mounted,  in  turn,  on  a  metal 
bracket  which  extended  into  the  flame. 

A  hypodermic  syringe  was  employed  to  mount  the  slurry  droplets  on 
the  thermocouple  probe.  The  syringe  was  capable  of  producing  drops  in 
the  500-5000  ym  diameter  range.  The  flame  was  deflected  from  the 
mount  region  until  the  droplet  was  in  place.  The  deflector  was 
constructed  from  a  200  x  305  mm  sheet  of  stainless  steel  3  mm  thick. 
The  carbon-particle  combustion  process  was  initiated  by  removing  the 
deflector  with  a  pneumat ic-cyLinder  arrangement. 
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2.2.2  Instrumentation 

2.2.2. 1  Drop  Environment  Measurements 

The  drop  environment  at  each  test  location  within  the  flame 
was  characterized  in  previous  studies  [149,150,151].  These  studies 
involved  measurements  of  mean  velocity,  temperature  and  concentration. 

2. 2. 2. 1.1  Velocity  Measurements 

Velocity  measurements  along  the  centerline  of  the  propane 
flame  were  conducted  with  a  single-channel  laser-Doppler  anemometer 
(LDA).  An  equipment  list  for  the  LDA  is  supplied  in  Table  1.  A 
schematic  of  the  LDA  arrangement  is  shown  in  Figure  7.  This  unit 
employed  a  helium-neon  laser  operating  in  the  dual-beam,  forward- 
scatter  mode.  ■  The  laser  power  was  50  raW.  The  sending  and  receiving 
optics  had  a  focal  length  of  241  mm  with  an  11.6°  angle  between  the 
beams.  The  aperture  diameter  of  the  photodetector  was  0.256  mm.  The 
signal  was  focused  on  the  photodetector  with  a  200  mm  lens.  This 
produced  an  ellipsoidal  measuring  volume  3.0  mm  in  length  with  a 
diameter  of  0.31  mm.  A  frequency  shifter  was  used  so  that  flow 
reversals  could  be  detected.  Interference  of  the  light  from  the  flame 
was  eliminated  with  a  band-pass  filter. 

The  LDA  actually  measures  the  velocity  of  small  seeding  particles 
in  the  flow.  The  soot  that  formed  in  the  flames  provided  a  portion  of 
the  seeding  particles  in  the  region  of  combustion.  In  order  to  avoid 
biasing  the  signal,  seeding  was  also  added  to  the  surrounding  air  in 
the  region  where  measurements  were  being  made.  These  seeding 


Table  1. 

Component 

Helium-Neon  Laser 

Integrated  Optics 

Frequency  Shifter 

Photodetector 

Tracker 

RMS  Voltmeter 

Dual  Beam  Oscilloscope 

Integrating  Voltmeter 


Bandpass  Filter 


Equipment  List 


Manufacturer 

Model 

Spectra  Physics 

125A 

Thermo-Systems 

900 

Thermo-Systems 

980 

Thermo-Systems 

960 

Thermo-Systems 

1090 

Thenno-Sys  terns 

1076 

Tektronix 

912 

Hewle  1 1 -Packard 

240IC 

Thermo-Sy stems 
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particles  were  prod1'''*'*  as  oil  vapor  condensed  in  an  air  stream 
passing  through  an  oil  heater.  The  seeding  rate  was  controlled  by 
varying  the  heat  flux  and  air  flow  rate.  The  air-vaporized  oil 
mixture  was  passed  through  a  condenser  tube  upon  leaving  the  oil 
heater.  The  length  of  the  condenser  tube  could  be  varied  so  that  the 
air-oil  droplet  mixture  entered  the  test  air  at  approximately  the 
ambient  temperature.  The  oil  heater-condenser  tube  apparatus  was 
positioned  so  that  the  particles  approached  the  combusting  jet  at  the 
approximate  entrainment  velocity. 

As  the  seeding  particles  pass  through  the  measuring  volume  they 
scatter  light.  The  photomultiplier  picks  up  signals  and  transfers 
them  to  the  frequency  shifter  (which  also  supplies  power  for  the  Bragg 
cell)  and  then  to  the  tracker.  The  tracker  process  and  analyzes  the 
signals  and  then  sends  the  verified  signals  to  the  integrating  digital 
voltmeter,  the  true  RMS  meter,  and  the  oscilloscope  (which  monitors 
the  quality  of  signals). 

The  frequency  tracker  provides  a  signal  proportional  to  the 
instantaneous  flow  velocity  in  the  measuring  volume.  This  signal  can 
be  processed  similar  to  a  hot-wire  detector  signal.  Mean  velocities 
were  found  by  integrating  the  signal,  fluctuating  velocities  could  be 
found  with  a  true  RMS  meter.  Various  velocity  components  and  the 
Reynolds  stress  could  be  found  by  rotating  the  plane  of  the  laser 
beams.  The  data  reported  at  each  location  was  the  average  of  at  least 
three  measurements  each  integrated  over  a  period  of  one  minute.  At 
axial  locations  close  to  the  injector,  the  tracking  signals  were  too 
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infrequent  to  maintain  a  continuous  signal.  The  occurrence  of  this 
discontinuous  signal  ruled  out  the  use  of  the  integrating  voltmeter. 
Therefore,  150  individual  signals  were  recorded  from  the  tracker  and 
their  average  was  used  to  determine  the  velocity  quantities. 

2.2.2. 1.2  Temperature  Measurements 

The  mean  gas  temperatures  along  the  centerline  of  the 
turbulent-diffusion  flame  were  measured  using  the  thermocouple  probe 
illustrated  in  Figure  8.  The  temperature  probe  was  constructed  from 
25  Vim  diameter  platinum/platinum-10%  rhodium  wires  which  were  spot- 
welded  onto  750  ym  diameter  lead  wires  of  the  same  material.  The 
fine-wire  thermocouples  were  manufactured  by  Omega  Engineering, 
Incorporated.  The  25  ym  thermocouple  wire  was  used  when  possible; 
however,  at  positions  close  to  the  injector,  the  temperature  was  high 
enough  to  cause  the  thermocouple  probe  to  fail.  When  this  occurred, 
thermocouple  wire  with  a  50  ym  diameter  was  used. 

The  maximum  radiation  correction  for  this  temperature  probe  was 
less  than  35°C  over  the  test  range.  The  reference  junction  was  at  the 
ambient  temperature,  shielded  from  the  flame.  An  integrating  digital 
voltmeter  was  used  to  average  the  signal  over  a  two-minute  time 
period. 


2.2.2. 1.3  Concentration  Measurements 

Mean  concentration  measurements  were  conducted  along  the 
centerline  of  the  turbulent-diffusion  flame.  The  gas  samples  were 
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extracted  from  the  flow  at  nearly  Isokinetic  conditions  using  the 
water-cooled  probe  illustrated  in  Figure  9.  Both  the  inside  and  the 
outside  tubes  were  stainless  steel,  and  the  inside  tube  had  a  2  mm 
internal  bore.  The  tip  of  the  probe  tapered  back  from  3  mm  to  6  mm  in 
external  diameter  over  a  distance  of  about  10  mm.  The  internal  tube 
was  cooled  by  water  which  gave  substantial  cooling  within. 2  mm  of  the 
tip.  The  temperature  of  the  cooling  water  was  maintained  constant  at 
about  339  K  (+  3  K)  in  order  to  avoid  condensation  of  fuel  and  water 
vapor  in  the  probe.  The  flows  were  isokinetically  sampled  by  applying 
a  suction  to  the  sampling  probe  and  measuring  the  sampling  flow  rate. 
The  flow  rate  through  the  sampling  probe  was  adjusted  to  match  the 
local  flow  rate  determined  from  the  velocity  measurements.  The 
sampling  flow  rate  was  measured  with  a  Precision  Scientific  wet-test 
meter  (28.3  ml/rev.). 

The  samples  were  analyzed  with  a  Perkin  Elmer,  Model  880,  gas 
chromatograph  using  a  hot-wire  detector.  A  1.0  ml  gas-tight  syringe 
was  used  to  transfer  the  sample  from  the  sampling  port  to  the  gas 
chromatograph.  The  response  of  the  hot-wire  thermal-conductivity 
detector  was  recorded  with  a  Leeds  and  Northrup  Speedomax  Strip  Chart 
Recorder. 

A  Carbosieve,  type  B,  60-80  mesh  column  was  used  to  measure  the 
propane  concentration  in  the  combust ing-gas  jet.  The  column  was 
.92  m  x  3.2  mm  O.D.  stainless  steel  and  was  manufactured  by  Supelco, 
Incorporated.  The  flow  rate  of  the  carrier  gas,  helium,  was 
maintained  at  50  ml/min.  The  column  temperature  was  varied  from 
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Figure  9.  Sketch  of  the  probe  used  for  gas -concent re c J cm 

measurements  In  the  turbulent  diffusion-flame  test. 
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175°C  co  225°C  with  temperature  programing  at  30oc/min.  The  column 
was  calibrated  using  calibration  gases  manufactured  by  the  Scott 
Company.  The  sample  size  used  for  the  propane  concentration 
measurements  was  0.5  ml. 

The  permanent  gases  (oxygen,  nitrogen,  carbon  monoxide,  carbon 
dioxide,  hydrogen  and  water  vapor)  were  analyzed  using  a' Carbosieve, 
type  S,  100-120  mesh,  2.1  m  x  3.2  ran  O.D.  stainless-steel  column. 
The  gas  chromatograph  was  programmed  for  a  four  minute  hold  at  30°C, 
followed  by  temperature  programming  from  35°c  to  175°C  at  24°C/min. 
Helium  was  used  as  the  carrier  gas  at  a  flow  rate  of  25  ml/min.  The 
column  was  calibrated  for  oxygen,  nitrogen,  carbon  monoxide  and  carbon 
dioxide  using  calibration  gases  manufactured  by  the  Scott  Company. 
The  column  was  calibrated  for  water  vapor  by  extracting  samples  over  a 
liquid  water  surface  maintained  at  a  constant  temperature  at 
atmospheric  pressure.  The  saturation  pressure  corresponding  to  the 
liquid  temperature  provided  the  mole  fraction  of  water  vapor  in  the 
air-vapor  mixture.  The  column  was  calibrated  for  hydrogen  by 
preparing  samples  of  hydrogen  and  air  whose  composition  was  known  on  a 
volume  basis.  In  each  case,  the  sample  size  used  was  0.5  ml. 

2. 2. 2. 2  Drop  Measurements 

Two  slurry  fuels,  provided  by  Suntech  Group,  were  examined 
during  the  tescs.  The  fuels  consisted  of  a  medium-thermal  carbon 
black,  having  an  ultimate  particle  size  of  0.3  ym,  with  JP-1C  as  the 
iiquid  fuel,  c. f . ,  Reference  2  for  more  complete  specifications.  The 


bulk  of  the  testing  was  conducted  with  a  noncatalyzed  slurry 


with 
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50.41  dispersed  carbon  by  weight.  A  catalyzed  s  Lurry,  which  had  49 .  IZ 
dispersed  carbon  by  weight  which  contained  a  proprietary  lead  compound 
as  a  cacalyst,  was  also  considered.  The  properties  of  the  carbon 
slurry  fuels  are  summarized  in  Table  2.  While  these  fuels  are 
representative  of  carbon  slurries,  fuel  development  efforts  are 
continuing  and  they  do  not  necessarily  represent  the  most  attractive 
formulations . 

Samples  of  the  carbon  agglomerates  were  obtained  by  quenching  the 
process  by  redeflecting  the  flame.  Carbon  residue  left  on  the  probe 
after  combustion  was  examined  with  a  scanning  electron  microscope. 
International  Scientific  Instruments,  model  H-7 .  Permanent  records  of 
appearance  of  the  residue  were  obtained  with  Polaroid,  type  52,  fine- 
grain  film. 

2.2. 2. 2.1  Diameter  Measurements 

The  drop  combustion  process  was  observed  with  a  16  mm, 
Photosonics,  model  16-B,  motion-picture  camera.  The  camera  was 
powered  by  a  Kepco,  SM  36-5  AM  D.C.  power  supply.  The  film  speed  was 
indicated  with  a  timing  light  on  the  camera,  activated  with  an  Adtrol 
Electronic  pulse  generator,  model  501.  Kodak  Plus— X  reversal  film  was 
used  foe  all  tests. 

Backlighting  for  shadowgraph  measurements  was  provided  by  a  Pek, 
model  401A,  75W  mercury  arc  lamp.  The  light  from  the  arc  was 
collimated  and  directed  toward  the  drop  location.  A  diffuser  screen 
was  employed  behind  the  drop  to  equalize  the  light  intensity  of  the 
background.  The  background  intensity  was  adjusted  so  that  envelope 
flames  around  the  drop  could  be  observed. 
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Table  2.  Properties  of 

the  Carbon  Slurry 

Fuels* 

Noncatalyzed 

Catalyzed*3 

Designation 

790-928 

790-942 

Liquid 

JP-10 

JP-10 

c 

Dispersed  Carbon  (wt  %) 

50.4 

49.2 

aFuels  supplied  by  R.  S.  Stearns,  Suatech,  lac., 

P.  0.  Box  1135,  Marcus  Hook,  PA  15061. 

^The  catalyst  is  a  proprietary  lead  compound  [2]. 

£ 

The  dispersed  carbon  is  a  medium  thermal  carbon  (carbon 
black)  having  an  ultimate  particle  size  of  0.3  pm. 
Further  properties  of  this  material  are  provided  in 
Reference  2. 
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Vanguard  notion  picture  analyzer-  Photographs  or  objects  of  known 
size  at  the  drop  location  provided  a  calibration  of  distances  on  the 
f  i  in. 


2. 2. 1.1. 2  Temperature  Measurements 

The  carbon-agglomerate  temperature  was  measured  during  the 
combustion  process  with  the  thermocouple  temperature  probe  discussed 
in  the  previous  section  and  a  Nicolet,  model  206,  analog-to-digitai 
recording  oscilloscope.  The  Nicolet  has  a  memory  capacity  of  4090 
data  points.  The  combustion  process  of  a  large  carbon  agglomerate  was 
completed  within  10  seconds.  This  gave  a  time  between  samples  less 
than  2.5  milliseconds  which  proved  to  be  more  chan  adequate  to  record 
any  temperature  variation  of  the  carbon  particle  during  the  combustion 
process. 

2.2.3  Experimental  Conditions 

The  bulk  of  the  testing  was  undertaken  at  locations  along  the 
centerline  of  the  gaseous  propane  diffusion  flame,  between  75  and  600 
injector  diameters  from  the  injector.  No  external  ignition  source  was 
employed  for  these  tests.  A  single  operating  condition  was  employed 
for  all  tests  and  is  summarized  in  Table  3.  A  photograph  of  Che  flame 
appears  in  Figure  10. 

A  summary  of  mean  velocities,  temperatures  and  compositions  in 
the  flow  is  presented  '  Figure  il .  A  complete  analysis  of  this  flame 
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Table  3.  Turbulent  Diffusion-Flame 

Characteristics3 

Fuel : 

propane 

Fuel  f.Tou  rate: 

176  mg/s 

Initial  fuel  jet  velocity 

88.7  m/s 

Initial  fuel  jet  diameter: 

1.194  mm 

Jet  Reynolds  number: 

23600 

Orientation: 

vertical  (upward) 

Flame  Height : 

460  mm  (visual) 

Injector  Thrust: 

15.6  mN 

Hydrogen  flow  rate: 

0.126  mg/s 

Ambient  and  injector  inlet  temperature: 

296  K 

Ambient  pressure: 

97  kPa 

aThe  nozzle  is  Spraying  Systems  Company  model  1/4J  2050  fluid 
nozzle  and  a  model  67147  air  nozzle,  air  atomizing  injector. 
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Figure  10.  Photograph  of  th<i 


turbulent  sas-d if fusion  flaac. 
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X  /  d 


Figure  11.  Flow  properties  of  the  turbulent  propane  diffusion  flame. 

Note:  x/d  indicates  injector  diameters  downstream  from 
injector,  and  uc/u  indicates  centerline  velocity 
normalized  by  injector  velocity. 


TEMPERATURE  (K) 


66 


appears  elsewhere  [149,150,131].  The  chosen  operating  condition  gave 
a  variation  of  equivalence  ratio  of  approximately  0.2  to  1.4,  and  a 
temperature  range  of  approximately  600  to  1900  K.  along  the  flame 
centerline.  Initial  slurry-drop  sizes  were  in  the  range  of 
500-1500  ym. 

2. 3  Laminar  Premixed  Flame  Apparatus 

2.3.1  Test  Apparatus 

A  sketch  of  the  laminar-premixed  flame  apparatus  appears  in 
Figure  12.  The  arrangement  basically  consists  of  a  slurry-droplet 
generator,  a  convergent-flow  section,  a  liquid  evaporation  section, 
and  a  flat-flame  burner.  The  monodispersed  drops  produced  in  the  drop 
generator  are  accelerated  in  the  convergent-flow  section,  and  the 
liquid  fuel  evaporated  from  the  drop  in  the  heated  section  leaving 
behind  an  agglomerate  of  solid-carbon  particles.  The  carbon 
agglomerates  then  passed  through  a  flat  flame  where  they  reacted.  A 
sampling  probe  placed  above  the  flat-flame  burner  collected  the  carbon 
agglomerates  for  further  analysis. 

The  laminar-premixed  flat  flame  was  fueled  with  a  mixture  of 
gaseous  nitrogen,  oxygen,  carbon  dioxide,  carbon  monoxide,  methane  and 
hydrogen  using  rotometers  to  blend  in  the  desired  amount  of  a  specific 
gas.  The  mixing  pressure  was  monitored  with  a  Heise  absolute-pressure 
gauge  having  a  pressure  range  of  0.0-0.41  MPa. 

A  Matheson  rotameter,  model  605,  was  used  to  meter  the  nitrogen 
flow  rate.  The  oxygen,  carbon  dioxide  and  methane  gas  flow  rates  were 
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Figure  12.  Schematic  diagram  of  the  particle  generator 
and  laminar  premixed  flame  apparatus. 
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all  meter^c  with  Matheson  model  604  rotameters,  while  carbon  monoxide 
and  hydrogen  were  metered  with  a  Matheson  model  603  rotameter.  All  of 
the  rotameters  were  calibrated  with  a  283  ml/rev  Precision  Scientific 
Company  wet-test  meter. 

The  pressures  of  the  gases  entering  the  rotameters  were 
controlled  using  model  8H  Matheson  two-stage  pressure  regulators  for 
nitrogen,  carbon  monoxide,  methane  and  hydrogen.  The  oxygen  pressure 
was  controlled  using  a  model  92-100  Harris  two-stage  regulator,  and 
the  carbon  dioxide  pressure  was  controlled  using  a  model  1H  Matheson 
two-stage  regulator.  Each  pressure  regulator  had  an  output  pressure 
capacity  of  0-1400  kPa. 

Gas  mixing  was  aided  by  flowing  the  gases  through  a  length  of 
10  mm  I.D.  tubing  equal  to  approximately  150  tube  diameters.  The  gas 
mixture  then  entered  the  mixing  chamber  through  a  6.4  n an  diameter 
inlet  on  the  side  of  the  mixing  chamber  25  mm  from  the  bottom.  The 
gaseous  mixing  chamber  was  constructed  from  a  51  ram  nominal  diameter 
schedule  40  stainless-steel  welded  pipe  and  had  a  length  of  300  mm. 
An  endcap  over  the  lower  opening  and  a  152  mm  diameter  flange  around 
the  upper  opening  were  constructed  from  13  mm  thick  stainless-steel 
plates. 

To  provide  a  thoroughly  mixed  fuel  gas  with  a  uniform-velocity 
profile,  several  layers  of  small  beads  separated  by  a  void  space  were 
provided  in  the  gaseous  mixing  chamber.  Upon  entering  the  nixing 
chamber,  the  gas  encountered  a  50  mm  high  void  space  followed  by  a 


50  ran  layer  of  3  mm  diameter  solid  spherical  glass  beads. 


The  gas 
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mixture  then  encountered  another  50  am  high  void  space,  another  50  no 
layer  of  the  saae  cype  of  glass  bead,  and  another  50  am  high  void 
space.  Finally  the  gas  encountered  a  50  nun  layer  of  2  am  diaaeter 
stainless-steel  balls.  The  two  layers  of  glass  beads  and  the  layer  of 
stainless-steel  balls  were  supported  by  a  20  mesh  size  stainless-steel 
screen  which  was  spot-welded  to  the  sides  of  the  mixing  chamber.  A 
much  finer  stainless  steel  screen  with  a  mesh  size  of  100  was  placed 
immediately  on  top  of  the  20  mesh  size  screen.  The  glass  beads  or 
stainless  steel  balls  were  then  in  turn  placed  upon  the  100  mesh  size 
screen. 

After  passing  through  the  stainless-steel  balls,  the  fuel  gas  was 
ignited  and  a  flat  flume  was  stabilized  near  the  top  surface  of  the 
balls.  The  flat  flame  was  measured  to  be  approximately  1-2  mm  thick. 
The  flame  appeared  completely  uniform,  and  was  stable  at  all  operating 
conditions. 

The  ability  of  mixing  the  fuel  gases  together  in  ary  combination 
enables  the  system  to  operate  under  a  wide  range  of  conditions. 
Equivalence  ratios  between  0.2  and  1.4  were  obtained  during  the  tests. 
The  temperatures  of  these  flames  could  be  varied  between  1600  K.  and 
the  softening  temperature  of  stainless  steel.  By  adding  either  carbon 
monoxide  and  carbon  dioxide,  or  hydrogen  along  with  the  methane,  the 
C/H  ratio  can  be  varied  over  a  very  wide  range.  Finally,  the 
stainless-steel  balls  inhibit  flash-tack  and  aid  in  flame  attachment 
such  that  the  gas  velocity  above  the  flat-flame  burner  can  be  varied 
between  0.5  and  3.0  m/s. 
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An  ignitor  was  constructed  over  the  flat-flame  burner  from 
0.25  vsr  diameter  nicrome  heating  wire.  The  heating  wire  was  powered 
from  an  Ohmite  120  volt,  A. 75  amp.  Variac;  and  a  12.6  volt,  3  amp. 
transformer.  Once  the  flame  was  ignited,  a  53  mm  I.D.  quartz  tube  was 
placed  on  the  burner  surface  using  asbestos  coated  forceps.  The 
quartz  tube  surrounded  the  flame  to  prohibit  entrained  air.  from  mixing 
with  the  prooucts  of  combustion  and  lowering  the  flame  temperature. 
The  quartz  tube  also  helped  maintain  a  uniform  high  temperature  of  the 
products  of  combustion  far  above  the  burner.  The  quartz  tube  was  cut 
into  sizes  ranging  from  50  mm  to  200  tma  in  50  mm  increments  to 
facilitate  the  access  of  measuring  probes  to  various  locations  above 
the  flat  flame. 

The  droplet  generator  consisted  of  a  TSI  model  3050  Berglund-Liu 
vibrating-orif ice  monodisperse-aerosol  generator,  a  Sage  model  355 
continuously-variable  syringe  pump,  and  a  Hewlett-Packard  model  331QB 
function  generator.  Three  different  orifices  were  used  in  the 
Berglund-Liu  drop  generator:  a  TSI  model  10356  100  pm  orifice,  a  TSI 
model  10355  50  ym  orifice  and  a  TSI  model  10353  20  ym  orifice.  Both 
a  10  and  50  ml  reusable  mult  if it  syringe  with  metal  luer  tips  were 
used  with  the  syringe  pump. 

A  mixture  of  methane,  oxygen,  carbon  monoxide  and  nitrogen  was 
used  in  the  Berglund-Liu  droplet  generator  as  the  dilution  and 
dispersion  flows.  The  fuels  were  mixed  in  proportions  to  yield  the 
same  equivalence  ratio  end  C/H  ratio  as  the  fuel  mixture  in  the  flit- 
flame  burner  assembLy.  The  dispersion  flow  separates  the  droplets  as 
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they  are  famed  and  transports  them  away  from  the  orifice  to  prevent 
clogging.  The  dilution  flow  is  used  to  vary  the  droplet 
concentration.  The  dispersion  flow  makes  up  a  very  small  fraction  of 
the  total  flow  and  the  dilution  flow  the  majority  of  the  total  flow. 
The  flow  rates  of  the  methane,  oxygen,  carbon  monoxide  and  nitrogen 
were  all  metered  with  Matheson,  model  602,  rotameters,-  which  were 
calibrated  using  a  Precision  Scientific  Company  wet-test  meter 
(28.3  ml/rav).  The  pressures  of  the  methane,  carbon  monoxide  and 
nitrogen  flows  were  controlled  with  Matheson,  model  8H,  two-stage 
pressure  regulators.  The  oxygen  flow  pressure  was  controlled  with  a 
Harris,  model  92-100,  two-stage  pressure  regulator. 

Only  the  noncatalyzed  slurry  with  50. 42  dispersed  carbon  by 
weight  in  JP-10  liquid  fuel  was  used  for  laminar-premixed  flame  tests. 
In  order  to  generate  slurry  droplets  with  the  Berglund-Liu  droplet 
generator  without  clogging  the  orifice,  the  slurry  had  to  be  thinned. 
N-pentane  was  mixed  with  the  slurry  fuel,  and  it  was  determined  that 
the  carbon  particles  remained  suspended  in  the  mixture  for  several 
hours.  Thus,  mixtures  of  carbon  slurry  and  n-pentane  ranging  from  one 
to  five  parts  n-pencane  per  part  of  volume  of  slurry  were  used  in  the 
droplet  generating  system. 

The  syringe  pump  was  adjusted  to  the  minimum  velocity  necessary 
to  establish  a  liquid  jet  emerging  from  the  orifice  in  the  Berglund- 
Liu  droplet  generator.  The  frequency  of  the  signal  generator  was  set 
at  a  value  between  1000  and  4000  Hz.  Thus,  between  1000  and  4000 
drops  per  second  were  formed  ranging  in  initial  diameter  between  24 
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and  270  ym.  Using  the  50  ml  syringe  allowed  continuous  operation  for 
up  to  3-5  hours. 

The  Berglund-Liu  droplet  generator  exhausted  the  slurry  droplets 
into  a  50  am  length  of  a  115  mm  diameter  plexiglass  tube.  The 
plexiglass  tube  was  connected  to  a  stainless-steel  45°  convergent 
cone,  which  was  attached  to  the  plexiglass  with  silicon  sealant.  The 
cone  was  40  ssn  high,  had  a  115  an  diameter  opening  at  the  base  and  a 
3.4  nm  diameter  opening  at  its  tip.  The  outlet  of  the  cone  was  welded 
to  a  tube  having  an  inner  diameter  of  3.4  mm  and  a  length  of  610  mm. 
The  top  end  of  this  tube  was  welded  to  a  tube  with  a  1.8  mm  inner 
diameter  and  a  length  of  460  mm.  This  tube  passed  through  a  Conax 
Corporation  thermocouple  sliding  seal  at  the  bottom  of  the  mixing 
chamber.  The  top  end  of  the  1.8  am  diameter  tube  was  flush  with  the 
surface  of  the  flat-flame  burner.  Great  care  was  taken  in  the 
construction  of  the  droplet  guide  to  assure  a  smooth  flow  transition 
between  the  individual  components. 

Nichrome  heater  wire  and  fiberfrax  ceramic  insulation  were 
wrapped  around  the  outside  of  the  3.4  mm  diameter  section.  Thus,  when 
a  particle  left  the  Berglund-Liu  droplet  generator  it  was  accelerated 
through  the  convergent  core  section  of  the  droplet  guide.  Next,  the 
droplet  entered  the  3.4  mm  diameter  section  which  was  heated  to  600  K 
by  the  nichrome  heating  wire  and  a  120  volt,  4.75  amp.  Ohraite  variac. 
In  this  region  the  liquid  fuel  evaporated  leaving  behind  an 
agglomerate  of  solid-carbon  particles.  The  evaporating  liquid  aiso 
served  to  further  accelerate  the  agglomerate.  When  the  gas  flow 
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emerged  from  Cop  of  the  1.8  as  diameter  tube  it  had  the  same 
equivalence  ratio,  CH  ratio  and  velocity  as  the  surrounding  gas  of  the 
flat-flame  burner.  The  particle  velocity  was  somewhat  lower  than  the 
flow  velocity  due  to  the  gravity  effects  on  the  particle.  The  initial 
size  of  the  carbon  agglomerates  leaving  the  1.8  n*a  diameter  tube 
ranged  between  10  and  140  m.  The  particles  continued'  to  flow  up 
through  the  quartz  tube  as  they  reacted. 

Figure  13  is  a  time  exposure  photograph  of  the  carbon 
agglomerates  passing  through  the  flat— flame  burner  and  reacting.  The 
quartz  tube  described  as  shown  in  Figure  12  was  not  in  place  for  this 
photograph,  and  as  a  result  the  particles  began  to  cool  upon  exiting 
the  hot  potential  core  of  the  flame.  Although  the  particles  do 
disperse  somewhat,  the  photograph  show  that  the  particles  travel  in  a 
straight  line  producing  a  streak  approximately  3  mm  wide  at  the  base 
for  an  initial  particle  diameter  of  73.4  a. 

2.3.2  Instrumentation 

2.3.2. 1  Particle  Environment  Measurements 

The  particle  environment  was  characterized  by  measuring  the 
velocity,  temperature  and  concentration  of  the  gas  along  the 
centerline  of  the  laminar— premixed  flame. 

2.3.2. 1.1  Velocity  Measurements 

Velocity  measurements  were  made  along  the  centerline  of  the 
lamina r-preaixed  flat  flame  with  basically  the  same  laser-Doppler 
anemometer  as  for  the  turbulent-diffusion  flame  velocity  measurements. 
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However,  for  the  velocity  measurements  of  the  laainar-premixed  fla-e  a 
smaller  laser  was  employed  (Spectra-Physics  Helium-die  or.  laser,  model 
120,  with  a  power  output  of  5.0  mW),  the  laser  was  operated  in  the 
back-scatter  mode,  and  the  tracker  was  replaced  by  a  burst-counter 
(TSI  model  1934  input  conditioner,  TSI  model  1985  timer,  TSI  model 
1942  readout  and  TSI  model  1988  analog  output). 

The  laainar-prenixed  flame  did  not  provide  a  sufficient  seeding 
particle  density  for  the  LDA  to  give  gas  velocity  measurements.  In 
order  to  increase  the  seeding  concentration  of  the  flow,  a  particle 
generator  based  on  a  reversed-cyclone  separator  described  by  Glass  and 
Kennedy  [152]  was  constructed.  Micron-sized  aluminum-oxide  particle 
(TSI  model  10083)  were  generated  since  they  have  acceptable  refractory 
properties  to  permit  their  use  in  combusting  flows. 

The  reverse-cyclone  separator  particle  generator  employed  in  this 
investigation  is  illustrated  in  Figure  14.  The  overall  height  of  the 
seeder  was  380  sga.  During  the  particle  generator  operation  a  30—60  am 
layer  of  aluminum-oxide  seeding  particles  were  stored  at  the  bottom  of 
the  particle  generator.  The  dry  filtered  gas  mixture  entered  the 
seeder  horizontally  and  tangentially  through  a  5  mra  diameter 
stainless— steel  tube  280  ram  from  the  bottom  of  the  particle  generator. 
The  gas  traveled  tangentially  around  the  particle  generator  and  mixed 
with  the  seeding  particles.  The  air  stream  and  entrained  particles 
then  left  the  particle  generator  through  a  20  mm  inner  diameter 
plexiglass  tube  located  at  the  top  of  the  particle  generator.  The 
entrance  to  the  seeding  particle  guide  was  at  the  sane  height  as  the 
tangential  gas  inlet. 


Figure  14. 


Schematic  diagraa  or  reverse-cvc lone-separator 
coEoustion-t laae  particle  seeder. 
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The  reverse-flow  cyclone  particle  seeder  was  operated  with  an  air 

o  q 

flow  rate  of  540  car/s,  which  gave  a  seeding  rate  of  1*1  x  10 

'y 

particles/m  .  Between  3.5  and  5.0  cmJ/s  of  the  seeded  flow  was  bled 
from  the  main  flow  and  injected  through  a  1.8  mm  inner  diameter  tube 
to  the  centerline  of  the  flat-flame  burner.  The  exit  velocity,  C/H 
ratio  and  the  equivalence  ratio  of  the  flow  was  matched  to  the  fuel- 
gas  flow  entering  into  the  flat-flame  burner.  The  remaining  portion 
of  the  seeded  flow  was  routed  to  a  collection  chamber  with  a  filter  to 
salvage  the  unused  aluminum-oxide  particles. 

2.3.2. 1.2  Temperature  Measurements 

The  gas  temperature  along  the  centerline  of  the  laminar- 
premixed  flat  flame  were  measured  using  the  thermocouple  probe 
illustrated  in  Figure  15.  The  temperature  probe  was  constructed  from 
25  pm  diameter  platinum/ platinum- 10%  rhodium  wires  which  were  spot- 
welded  onto  50  Um  diameter  lead  wires  of  the  same  material.  The 
thermocouple  wires  were  manufactured  by  Omega  Engineering, 
Incorporated.  The  temperature  probe  was  mounted  on  a  320  mm  Unislide 
traversing  mechanism  which  allowed  the  probe  to  be  positioned  at 
various  locations  along  the  centerline  of  the  laminar-premixed  flat 
flame. 

The  eraissivity  of  the  thermocouple  was  maintained  at  a  low  level 
of  approximately  0.22  +  0.02  by  coating  the  thermocouple  with  a  layer 
of  silica  [153].  This  was  accomplished  by  flowing  the  silicon  oil 
ladened  exhaust  of  a  model  1402  Sargent-Welch  Scientific  Company'  duo- 
seal  vacuum  pump  over  the  thermocouple  while  the  thermocouple  at 
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15.  Sketch  of  the  temperature  probe  used  in 
the  laminar  premixed  flame  apparatus. 


1800  K  was  protruding  from  a  1.8  mm  diameter  tube  into  the  flat-flame 
burner.  A  silica  layer  with  a  thickness  of  approximately  3  ym  was 
coated  onto  the  thermocouple  junction.  The  silica  layer  also 
eliminates  the  catalytic  decomposition  heating  of  the  platinum 
thermocouple  [153]. 

The  platinum/ plat inure- 10%  rhodium  thermocouple  -probe  was 
calibrated  with  a  50  ym  chromel/alumel  thermocouple  in  a  furnace 
environment  where  the  radiation  error  would  be  minimized.  The  furnace 
was  a  Lindberg  "Hevi-Duty",  .ype  54032-A,  furnace  with  a  power 
capacity  of  960  watts.  The  furnace  was  powered  by  a  Lindberg,  type 
59445  power  supply.  Both  thermocouples  agreed  to  within  6°C  of  each 
other  up  to  a  temperature  of  1500  K,  which  was  the  maximum  temperature 
of  the  furnace. 

The  maximum  radiation  correction  for  this  temperature  probe  was 
39°C  over  the  test  range.  The  reference  junction  was  at  the  ambient 
temperature,  shielded  from  the  flame.  A  Hewlett-Packard,  model  240 IC, 
integrating  digital  voltmeter  was  used  to  average  the  signal  over  a 
two-minute  period. 

2. 3.2. 1.3  Concentration  Measurements 

The  concentration  measurements  were  made  using  a  gas 
chromatograph  (Varian  model  3720,  incorporating  a  thermal-conductivity 
detector).  The  gas  samples  were  extracted  from  the  flow  at  nearly 
isokinetic  conditions  using  a  small  diameter  probe  shown  in  Figure  16. 
Both  the  inside  and  the  outside  tubes  were  stainless  steel,  and  the 
inside  tube  had  a  l  mm  internal  bore.  The  tip  of  the  probe  tapered 
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back  from  2  mm  to  6  mm  in  external  diameter  over  about  5  mm.  The 
internal  tube  was  cooled  by  water  which  gave  substantial  cooling 
within  2  a  of  the  tip.  The  temperature  of  the  cooling  water  was 
maintained  constant  about  338  K  (+  2  K)  to  avoid  the  condensation  of 
fuel  and  water  vapor  in  the  probe  prior  to  the  sampling  point. 

The  samples  were  analyzed  using  a  carbosieve  column  (Supelco, 
Carbosieve,  type  S,  100/120  mesh,  2.1  m  x  3.2  mm  stainless  steel; 
4  min.  hold  at  35°C;  programming,  25°C  to  175°C  at  25°C/min.)  for  the 
analysis  of  hydrogen,  oxygen,  nitrogen,  carbon  monoxide,  carbon 
dioxide  and  methane.  The  concentration  of  water  vapor  was  obtained  by 
a  conservation  of  element  analysis.  Helium  at  a  flow  rate  of 
30  nl/tnin.  was  used  as  carrier  gas.  The  output  of  the  hot-wire 
thermal-conductivity  detector  was  recorded  with  a  Varian  strip  chart 
recorder. 

The  Carbosieve  column  was  calibrated  for  hydrogen,  oxygen, 
nitrogen,  carbon  monoxide,  carbon  dioxide  and  methane  using  calibrated 
gases  manufactured  by  the  Scott  Company.  The  sample  size  used  for  the 
laminar-premixed  flat  flame  composition  measurements  was  0.5  ml. 

Since  the  products  of  combustion  from  a  premixed  methane  flame 
contain  a  large  mole  fraction  of  water,  the  water  must  be  condensed 
prior  to  gas  analysis  to  avoid  damaging  the  thermal-conduct ivi  ty 
detectors  in  the  gas  chromatograph.  The  sampling  system  that 
accomplished  this  water  condensation  process  is  illustrated  in  Figure 
17.  The  water-cooled  gas-sampling  probe  is  shown  in  position  above 
the  flat-flame  burner.  The  gas-sampling  probe  was  mounted  on  a  320  mm 


Unislide  traversing  mechanism  which  allowed  the  probe  to  be  positioned 
at  various  locations  along  the  centerline  of  the  laminar-premixed  flat 
flame.  The  sample  entered  the  probe  where  the  temperature  was 
approximately  338  K  and  the  pressure  60.8  kPa,  both  of  which  help  to 
quench  any  further  gas  reactions.  The  exhaust  from  the  gas  sampling- 
probe  was  directed  to  a  glass  cold  trap  (Supelco,  model'  2-2422)  to 
condense  a  portion  of  the  water  vapor  from  the  gas  flow.  The  flow  was 
then  directed  to  a  125  mi  glass  sampling  bulb  (Supelco,  model  2-2161). 
The  sampling  bulb  was  then  connected  to  a  Sargent-Welch  Scientific 
Company  duo-seal  vacuum  pump,  model  1402.  The  flow  rate  was 
controlled  with  a  Whitey  forged-body  regulating  valve,  model  SS-1RS4, 
and  metered  with  a  Varian  aerograph  bubble  meter. 

The  sampling  procedure  was  to  match  tha  gas  velocity  at  the 
sampling- probe  inlet  with  the  velocity  of  the  the  surrounding  flow. 
After  allowing  the  system  to  operate  for  approximately  5  minutes  to 
purge  the  system,  the  stopcock  on  the  downstream  end  of  the  glass 
sampling  bulb  was  closed.  After  a  few  seconds  the  system  reached 
pressure  equilibrium,  and  the  stopcock  on  the  upstream  end  of  the 
sampling  bulb  was  closed.  A  5  ml  sample  was  then  drawn  from  the 
septum  on  the  glass  sampling  bulb  with  a  0-1.0  ml  pressure-lok  series 
"c"  gas  syringe,  model  030033.  The  gas  sample  was  then  immediately 
taken  to  the  gas  chromatograph  for  analysis. 
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2.3.2. 2  Particle  Measurements 

2. 3. 2. 2.1  Temperature  and  Velocity  Measurements 

The  temperature  of  the  carbon  particles  as  they  reacted 
above  the  flat-flame  burner  assembly  was  measured  using  an  optical 
pyrometer.  The  pyrometer  used  was  a  Leeds  and  Northrup,  model  8622, 
with  a  temperature  range  of  1000-3150  K.  The  optical  pyrometer  was 
calibrated  using  the  platinum/platinua-IOZ  rhodium  thermocouple 
temperature  probe  illustrated  in  Figure  15. 

Still  photographs  of  the  reacting  carbon  particles  were  taken 
with  a  Graf lex  Optar  f/7.5,  203  mm  camera,  model  N-G55586.  Polaroid 
high-speed  instant  film,  type  57,  was  used.  The  camera  lens  was 
placed  1.6  m  from  the  film  face  which  gave  a  magnification  of 
7.5  to  1. 

In  order  to  obtain  the  still  photographs,  a  monodispersed  aerosol 
was  established  at  the  flat-flame  burner  face  in  a  completely  dark 
environment.  The  aperture  on  the  camera  optics  was  held  open  while  a 
light  pulse  illuminated  the  particles.  The  light  pulse  was  generated 
by  a  General  Radio  strobotac,  model  1538-A.  The  photograph  was  then 
analyzed  on  a  Spencer  microscope,  manufactured  by  the  American  Optical 
Companv,  to  observe  particle  spacing  and  carbon-reaction 
characteristics. 

Tlie  particle  velocity  measurements  were  made  with  the  LDA  system 
used  to  measure  the  velocity  of  the  gas  above  the  flat-flame  burner. 
The  particle  spacing  was  also  measured  with  the  LDA  system.  A 
Nicoiet,  model  206,  analog-to-digi tai  recording  oscilloscope  was 
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employed  to  measure  the  time  between  the  pedestal  signals  generated  by 
individual  particles.  Since  the  particle  velocity  was  known,  the 
distance  between  particles  could  be  determined. 

2. 3. 2. 2. 2  Diameter  Measurements 

The  probe  used  to  collect  the  carbon-agglomerated  particles 
for  analysis  is  shown  in  Figure  13.  The  probe  was  positioned  along 
the  flame  centerline  at  various  locations  above  the  flat-flame  burner 
with  a  320  mm  Unislide  traversing  mechanism.  The  collection  filter 
was  constructed  by  modifying  a  Gelman,  model  2220,  stainless-steel  in¬ 
line  filter  holder.  The  upstream  end  of  the  filter  holder  was 
modified  to  include  two  concentric  stainless-steel  cones  which  were 
Silver  soldered  to  the  upstream  face  of  the  filter  holder.  The  inner 
cone  had  a  6.4  mm  inlet,  a  15°  angle,  and  a  length  of  23.7  mm.  The 
outer  cone  had  a  12.7  mm  inlet,  a  30°  angle,  and  a  length  of  30.3  mm. 
Both  cones  had  an  approximate  thickness  of  1.6  am.  A  stainless-steel 
tip  was  silver  soldered  to  the  outer  core.  The  tip  was  12.7  mm  in 
diameter,  6.4  mm  in  length  and  tapered  to  a  4.7  mm  inlet  with  an  angle 
of  45°. 

A  nitrogen-quenching  flow  passed  through  the  outer  core  of  the 
filter  holder  from  a  6.4  mm  diameter  tube.  The  nitrogen  flowed 
between  the  two  concentric  cones  to  the  inlet  of  the  collection  probe 
to  quench  the  reacting  particles  as  they  entered  the  probe.  The 
particles  then  were  collected  on  47  mm  diameter  Alpha  metricei,  grade 
Alpha-iOU,  Gelman  membrane  filters,  model  00333.  These  filters  have  a 
pore  size  of  U.2  pm,  a  thickness  of  65  pm,  a  very  high  chemical 
resistance,  ana  a  maximum  operating  temperature  of  450  K. 


connected  to  a  Sargenc- 


The  ^article  collection  probe  was  than 
ViVLc  ,  Scientific  Company  duo-seal  vacuum  pump,  model  1402.  The 
velocity  at  the  inlet  to  the  collection  probe  was  matched  with  the  gas 
Velocity  at  the  sampling  location.  The  flow  rate  was  contrullec  with 
a  Whitcy  forged  body  regulating  valve,  aodel  SS-lRS-t,  and  teterta  witn 
a  Varian  aerograph  bubble  aetec.  A  sufficient  amount  of  nitrogen  gas 
was  raised  in  the  collection  filter  to  lower  the  gas  temperature  to 
4oU  n  and  quench  the  reacting  carbon  particles.  The  nitrogen  was 
metered  with  a  Macheson  o03  rotoneter,  and  the  pressure  was  controlled 
with  a  Matheson,  model  8H,  two-stage  pressure  regulator  with  an  output 
capacity  of  G-14J0  kPa.  The  operating  pressure  in  the  filter  was 
subataospheric  with  a  value  of  approximately  71  kPa  which  also  tended 
to  help  quench  the  reacting  carbon  particles. 

A  llupro  BUT  series  quick  response  toggle  operated  bellows  valve 
was  placed  in  the  flow  line  between  the  collection  probe  and  the 
vacuum  pump  in  order  to  be  able  to  operate  the  collection  filter  for  a 
short  time  period  of  a  few  seconds.  This  small  sampling  time  resulted 
in  widely  spaced  single  particles  on  the  collection  filter.  After  the 
sample  was  collected  the  filter  medium  was  removed  and  the  sices  of 
the  particles  collected  were  analyzed  on  a  scanning  electron 
microscope,  International  Scientific  Instruments,  model  Super  III— A. 
Permanent  records  of  the  particle  sizes  were  recorded  with  fine-grain 
Polaroid  film,  type  52. 
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The  diacetsr  quoted  is  the  arithaatic  average  of  six  diameter 
measurements  per  particle  each  spaced  30°  apart  and  passing  through 
the  center  of  the  particle.  At  least  eight  particles  were  examined  in 
order  to  obtain  an  average  particle  size. 

2. 3. 2.7.3  Mass  Measurements 

The  filter  collection  probe  illustrated  in  Figure  18  was 
also  employed  to  collect  particles  for  particle  mass  measurese-ts. 
However,  the  sampling  period  far  the  mass  neasurasents  was  increased 
to  between  10  seconds  and  2  hours,  and  a  stainless-steel  filter  medium 
was  substituted  for  the  Alpha  snetricel  filter  describee  earlier.  The 
stainiess-steel  filter  was  a  Dynalloy  filter,  model  X4,  manufactured 
by  the  Fluid  Dynamics  Branch  of  Brunswick  Technetics.  The  filters  had 
a  diameter  of  47  nan,  a  thickness  of  300  pa,  and  a  pore  size  of  3  -a. 
The  filters  were  preconditioned  by  baking  them  for  three  hours  at  a 
temperature  in  excess  of  1000  K  to  burn  off  any  contaminants  which  nay 
have  adhered  to  the  filter. 

After  the  carbon  particles  were  collected,  the  filter  was  removed 
and  placed  in  an  evacuated  reaction  chamber.  The  temperature  was 
elevated  to  500  K,  and  the  pressure  lowered  to  less  Chan  4000  Pa. 
These  conditions  were  maintained  for  several  hours  tc  remove  any 
condensibles  entrapped  within  the  particles.  The  dried  carbon 
particles  were  then  exposed  to  an  oxygen  and  hydrogen  environment. 
The  temperature  was  elevated  above  1000  K  and  the  pressure  increased 
to  35  Kjpa.  These  conditions  were  maintained  for  several  hours  until 
all  of  the  carbon  was  reacted.  By  measuring  the  carbon  dioxide  and 
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carbon  monoxide  concentrations  in  the  reaction  chamber,  the  mass  of 
carbon  initially  present  on  the  filter  was  determined.  Knowing  the 
total  number  of  particles  collected  during  the  collection  period  {from 
the  droplet  production  frequency)  gave  the  mass  of  carbon  per 
par'-icle.  The  pressure  of  the  reaction  chamber  was  kept 
subatmospheric  to  monitor  whether  the  system  developed  leaks,  since 
nitrogen  would  show  up  in  the  gas  analysis  of  the  reacted  gas  if  a 
le.ik  was  present. 

The  carbon  slurry  mixture  incorporated  several  additives  to 
assure  stability.  In  order  to  determine  whe . her  these  additives 
altered  the  carbon  mass  measurement,  a  slurry  mixture  consisting  of 
the  same  ultimate  sized  carbon  particles  (0.35  m)  and  nearly  pure 
JP-IO  (greater  than  99.95%)  was  prepared  in  this  laboratory.  The 
resulting  slurry  was  then  injected  through  the  flat-flame  burner  and 
combusted.  Particles  were  collected  at  various  axial  locations  for 
several  operating  conditions,  and  their  mass  analyzed.  The  mass 
measurements  of  the  carbon  agglomerate  with  and  without  the  additives 
present  were  round  to  agree  well  within  the  experimental  accuracy.  It 
is  therefore  assumed  in  the  analysis  that  the  additives  did  not 
significantly  alter  the  mass  measurements.  ' 

A  diagram  of  the  particle  dryer  and  reaction  chamber  system  is 
shown  in  Figure  19.  The  reaction  chamber  was  constructed  from  a  95  cc 
clear  fused  quartz  test  tube  with  an  I.D.  of  25  ram.  The  lower  120  mm 
of  the  test  tube  was  wrapped  with  nicroae  heating  wire  and  high 
temperature  fiberfrax  ceramic  fiber  insulation.  The  upper  80  mm  of 


Figure  19.  Schematic  diagram  of  the  particle-dryer  reaction-chamber  system 
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the  test  tube  was  not  insulated  which  kept  it  at  a  lower  temperature. 
The  tube  was  capped  with  a  two— hole  size  No.  5  rubber  stopper.  A 
septum  holder  and  septum  were  connected  to  one  of  the  holes  in  the 
rubber  stopper.  The  other  hole  was  connected  to  a  three-way  stopcock 
which  allowed  the  test  tube  to  be  connected  to  either  a  vacuum  pump  or 
an  oxygen  and  hydrogen  supply.  The  vacuum  pump  was  a  Sargent-Welch 
Scientific  Company  duo-seal  vacuum  pump,  model  1402.  The  oxygen 
pressure  was  controlled  with  a  Harris  two-stage  regulator,  model 
92-100;  and  the  hydrogen  pressure  was  controlled  with  a  Matheson  two- 
stage  pressure  regulator,  model  8H. 

A  pressure  gauge  was  also  incorporated  into  the  particle  dryer- 
reaction  chamber  system.  The  pressure  gauge  was  a  Wallace  and  Tiernan 
absolute-pressure  gauge,  model  FA129,  with  a  range  of  0  to  339  kPa. 
After  the  particles  were  dried,  as  described  earlier,  the  system  was 
evacuated.  Then  oxygen  was  added  until  the  pressure  reached  83  kPa, 
hydrogen  was  then  added  until  the  pressure  reached  85  kPa.  The 
reaction  chamber  temperature  was  elevated  above  1000  K,  and  after  all 
of  the  carbon  reacted  (approximately  two  hours),  the  temperature  was 
reduced  to  300  K.  A  5  ml  sample  was  extracted  from  the  reaction 
chamber  through  the  septum  with  a  0-1.0  ml  Pressure-Lok,  series  "c", 
gas  syringe,  model  030033.  The  gas  sample  was  then  analyzed  using  a 
Varian,  model  3700,  gas  chromatograph. 

2.3.3  Experimental  Conditions 

The  experimental  conditions  for  the  laminar  flat-flame  burner 


tests  are  summarized  in  Table  4, 


Four  equivalence  ratio  conditions 


Table  4.  Summary  of  Test  Conditions  for  the  Laminar  Flat-Flame  Burner  Tests 
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were  examined:  <£  *  0.2,  0.6,  1.0  and  1.4.  Three  temperature  ranges 
were  also  studied:  a  low  temperature  range,  »  1660  K;  an 
intermediate  temperature  range,  =  1815  K;  and  a  high  temperature 
range,  T-  =  1950  K.  The  temperature  at  a  given  equivalence  ratio  was 
increased  by  replacing  a  portion  of  the  combustible  mixture  with 
additional  nitrogen.  This  substitution  provided  an  increase  in  flame 
temperature  without  additional  thermal  decomposition  of  nitrogen 
[154]. 

The  73.4  ym  particles  burning  in  the  -  0.2  and  0.6  conditions 
were  completely  burned  within  the  200  mm  axial  length  limit  of  the 
quartz  tube.  However,  burning  was  incomplete  for  the  <p  -  1.0  and 
1.4  conditions.  Thus,  particles  with  smaller  initial  diameters  were 
tested  at  the  <p  =  1.0  test  condition  in  order  to  examine  variation 
of  the  combustion  process  with  particle  diameter. 


CHAPTER  III 


THEORETICAL  CONSIDERATIONS 

The  life  history  of  a  single  slurry  droplet  reacting  in  a 
combusting  environment  was  modeled  during  the  present  investigation. 
As  pointed  out  in  the  first  chapter,  the  liquid  fuel  in  a  slurry 
droplet  is  much  more  volatile  than  the  suspended  carbon  particles. 
Therefore,  it  is  reasonable  to  expect  that  the  liquid  fuel  will  be 
vaporized  first  and  that  a  single  carbon-agglomerated  particle  will  be 
left  behind.  Thus,  the  slurry-combustion  process  will  be  modeled  as 
two  separate  stages:  (1)  the  heat— up  and  gasification  of  the  i iquia 
fuel,  and  (2)  the  heat-up  and  reaction  of  the  carbon  agglomerate.  The 
discussion  will  begin  with  the  liquid-gasification  theory  and  its 
evaluation  since  this  analysis  is  relatively  simple  and  the  results 
serve  to  establish  methods  for  specifying  flow  properties  for  the  more 
complex  agglomerate-reaction  model. 

3.1  Liquid  Drop-Life-History  Model 

The  objective  of  the  liquid-gasification  analysis  is  to  predict 
the  variation  of  the  drop  temperature,  the  mass  of  the  liquid  and  the 
drop  diameter  as  a  function  of  time  for  the  period  where  liquid  is 
present  in  the  slurry.  The  analysis  was  relatively  conventional  since 
the  carbon  particles  remained  within  the  liquid  droplet-  The  approach 
generally  follows  methods  developed  during  earlier  work,  on  modeling 
drop  gasification  in  this  laboratory  [11,150,135,156]. 
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3.1.1  Assumptions 

Present  test  conditions  were  considered  involving  spherical 
particles  with  no  particle  interaction.  An  illustration  of  the  liquid 
gasification  model  is  presented  in  Figure  20. 

The  following  assumptions  were  employed  For  the  analysis  of  gas- 
phase  transport: 

1.  The  flow  field  around  the  drop  was  assumed  to  be  quasisteady, 
i.e.,  at  each  instant  of  time  the  flow  is  equivalent  to  a 
Steady  flow  for  the  same  boundary  conditions,  and  the 
velocity  of  the  surface  was  neglected.  Ambient  properties 
were  taken  to  be  local-mean  properties  ignoring  turbulent 
fluctuations. 

2.  Predicting  flow  around  spheres  exactly  is  impractical  for  use 

in  modeling  sprays;  therefore,  the  conventional  film  theory 
approximation  was  employed  to  treat  convection  [150].  This 
implies  that  gas-phase  transport  can  he  represented  as  a 
stagnant,  spherically-syrometric  layer  whose  outer  radius  is 
determined  from  empirical-convection  correlations  for 

spheres . 

3.  Only  diffusion  of  mass  by  concentration  gradients  was 

considered,  employing  an  effective-binary  diffusivity.  All 

species  were  assumed  to  have  equal  molecular  weights  and 
specific  heats,  and  constant  average  gas-phase  properties 
were  computed  at  a  mean  state  defined  as  follows: 

avg  p  1  ' 


(3.1) 
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Figure  20.  Illustration  of  transport  terms  used  in  the  theoretical  model. 


where  9  represents  both  temperature  and  species  mass 
fractions.  Predictions  vary  significantly  as  a  is  changed; 
a  single  value  was  chosen  to  provide  best  agreement  between 
predictions  and  measurements  over  the  test  range.  Methods 
developed  during  past  work.  were  employed  to  compute 
properties  [11,150,155,156].  A  summary  of  the  property 
evaluation  method  is  given  in  Appendix  A. 

4.  The  liquid  was  assumed  to  evaporate  without  the  presence  of 
an  envelope  flame  since  this  condition  is  most  representative 
of  drops  in  combusting  sprays  [11,150,155,156].  Thus, 
reaction  effects,  such  as  fuel  decomposition  or  oxidation, 
were  neglected. 

5.  The  gas  phase  was  assumed  to  be  transparent  to  radiation,  and 
radiation  from  the  drop  to  the  surroundings  was  neglected  due 
to  the  relatively  low  temperature  of  the  drop  liquid. 

6.  The  liquid  surface  was  assumed  to  be  in  thermodynamic 
equilibrium,  with  the  Clausius-Clapeyron  equation  providing 
the  relationship  between  liquid-surface  temperature  and  the 
concentration  of  vapor  at  the  surface. 

7.  The  solubility  of  gases  in  the  liquid  was  neglected  since 

this  effect  is  small  at  atmospheric  pressure  [156],  JP-10  is 

a  pure  hydrocarbon;  therefore,  the  analysis  only  considers 
the  evaporation  of  a  single  component.  Only  the  fuel  was 
considered  to  diffuse  in  the  gas  phase. 


8.  Radial  motion  of  the  liquid  surface  was  neglected. 
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8.  Xii ii  prsssurt  of  the  flow  field  surrounding  che  drop  was 
constant. 

10,  Dufour  and  Soret  effects  were  neglected. 

These  assumptions  are  generally  similar  to  those  employed  for 
drop-transport  analysis.  Their  justifications  are  discussed  more 
completely  in  Chapter  1  and  elsewhere  [11,1571. 

3.1.2  Vapor-Phase  Conservation  Equations 

Employing  the  assumptions  listed  in  the  previous  section,  the 
gas-phase  conservation  equations  are  as  follows  [11,150]: 

Conservation  of  Mass 


_d 

dr 


(Pr2v)  =  0 


(3.2) 


Conservation  of  Species 


h  •)  dYi 

—  (pr2(vY  -  D — - 
dr  v  F  siF 


})  =  0 


(3.3) 


Conservation  of  Energy 

fjnVcpd-y-jgn-o  0.4) 

The  boundary  conditions  for  these  equations  are  as  follows: 

r  =  dp/2,  Yf  =  YFgp,  T  =  Tp,  Pv  =  i*!  (3.3) 

dY 

r»dp/2,  Pv(l  -  YFgp)  +PD(-^)p=  0 


(3.*>) 
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rg* 


■«T> 


C3.7) 


where 


r. 


Sh 

Sh-2 


(3.8) 


r  . 

C  5 


r 


p 


Nu 

Nu-2 


(3.9) 


Nu  or  Sh  =  2  + 


[0.552  Re1/2  (Pr  or  Sc)1/3] 

[1  +  1.232/(Re  (Pr  or  Sc)4/3)]1/2 


(3.10) 


and 


Re  -  dpl««,  “  “p'p/W*  (3.11) 

The  terms  given  in  Equations  (3.8)  to  (3.11)  are  incorporated  to 
account  for  convection  effects.  The^e  expressions  were  obtained  by 
Faeth  and  Lazar  [41]  and  are  cited  in  References  11  and  157. 

The  boundary  condition  of  Equation  (3.6)  follows  from  the  fact 
that  the  nonfuel  gases  are  assumed  to  be  insoluble  in  the  liquid 
phase.  A  conventional  heat-transfer  coefficient  is  also  defined,  as 
follows: 

a(S>p  ■ h  it«  -  v-  <3a2> 

Integrating  Equation  (3.2)  and  applying  Equation  (3.5)  yields 
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pr-v  *  constant  =  d^  ft"/4- 
P  1 


(3.13) 


Integrating  Equation  (3.3)  and  applying  the  boundary  conditions  yields 


m'f  dp/p->  =  Sh  In  (1  +  S) 


(3.14) 


where 


B  » 


«Fg,  -  W 


F«P 

(*  -  YF8p> 


(3- 15) 


Integrating  Equation 


(3.4)  in  a  similar  manner  results 


Nu  In  1(1  +  B)1^  1 1 

1(1  +  B)Le_i  -  11  , 


in 


(3.16) 


where  the  Lewis  number,  Le,  is  based  only  on  the  specific  heat  of  the 
gaseous  fuel,  since  only  the  gaseous  fuel  diffuses  away  from  the 
particle. 

Equations  (3.14)  and  (3.16)  provide  the  basic  relationships  for 
mass  and  heat  transfer  in  the  presence  of  mass  diffusion  from  the 
drop.  These  equations  are  limited  to  sphericaliy-symmetric  flow  and 
the  effect  of  convection  is  handled  by  the  terms  Nu  and  Sh. 


3.1.3  Droplet  Evaporation  Model 

The  conservation  equations  for  the  drop  are  obtained  under  the 
assumptions  that  all  drop  properties  are  uniform  at  each  instant  of 
time.  Dufour  and  Sorer  effects  are  also  neglected.  With  these 
assumptions,  the  conservation  equations  are  as  follows  [ 1 1 J : 
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The  carbon  particles  are  assumed  to  remain  suspended  in  the 
liquid  fuel,  becoming  agglomerated  as  the  liquid  evaporates  and  the 
droplet  diameter  becomes  smaller.  When  the  solid-carbon  particle 
fores  art  agglomerate  which  is  as  tightly  packed  as  possible,  the 
lieu  id  continues  to  evaporate  through  the  agglomerate  pore  structure 
until  all  of  the  liquid  has  evaporated. 

A  uniform  temperature  at  each  instant  of  time  was  assumed  for  the 
liquid  droplet.  The  mass  fraction  of  the  liquid  fuel  was  also  assumed 
to  be  uniform  throughout  the  droplet  interior.  The  temperature  and 
concentration  profiles  assumed  by  the  model  are  sketched  in  Figure  21. 

3.1.4  Calculation  Procedure 

The  drop-life-history  calculations  were  begun  by  initializing 
the  particle  displacement,  velocity,  temperature,  diameter  and  carbon 
mass  fractions.  The  ambient-gas  mass  fractions,  temperature  and 
velocity  were  also  initialized  at  this  point.  Time  was  set  equal  to 
zero,  and  the  specific  heats  for  the  carbon  and  liquid  fuel  and  the 
liquid-fuel  density  were  calculated  for  the  particle  temperature. 
Employing  both  the  particle  and  ambient  temperature,  the  gas-phase 
thermophysical  and  transport  properties  were  calculated  as  discussed 
in  Appendix  A.  A  group  of  non-dimensional  parameters  were  then 
calculated  from  the  gas-phase  properties:  they  included:  the  Prandtl 
number,  the  Lewis  number  (based  on  the  vaporized  liquid-fuel  specific 
heat),  the  Schmidt  number,  the  Reynolds  number  from  Equation  (3.11), 
the  Nusselt  and  Sherwood  numbers  from  Equation  (3.10),  and  the  drag 
coefficient  from  Equation  (3.21). 
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EVAPORATION 


COMBUSTION 


Figure  21 


Li<  aid  gasification  and  combustion  models 
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The  droplet  mass  fraction  for  the  liquid  fuel,  solid  carbon,  and 
support  probe  were  then  calculated.  Next,  the  average-droplet  density 
and  specific  heat  and  the  droplet  diameter  were  calculated  from 
Equations  (3.22),  (3.23)  and  (3.24),  respectively.  The  mass  transfer 
driving  potential  was  obtained  from  Equation  (3.15),  and  the  mass  flux 
of  liquid  from  the  drop  surface  was  obtained  from  Equation  (3.14). 
The  convective  heat-transf er  coefficient  was  then  calculated  from 
Equation  (3. 1C,.  Finally,  the  time  derivatives  of  the  particle  mass, 
temperature,  velocity  and  displacement  were  obtained  from  Equations 
(3.17)  through  (3.20). 

The  minimum-droplet  diameter  (i.e.,  the  droplet  diameter  when  the 
carbon  particles  are  packed  around  the  support  probe  as  tightly 
together  as  possible)  was  computed  with  the  carbon  agglomerate  packed 
density  quoted  by  Bruce  et  al.  [2].  If  the  current  droplet  diameter 
was  smaller  than  this  computed  minimum  diameter,  the  droplet  diameter 
was  set  equal  to  this  minimum  diameter.  It  must  be  pointed  out  that 
even  though  the  minimum  diameter  has  been  reached,  there  still  may  be 
liquid  fuel  present  in  the  space  between  the  individual  carbon 
particles. 

If  the  mass  of  liquid  remaining  in  the  droplet  is  greater  than 
zero,  time  is  given  an  incremental  step  and  the  entire  calculation 
procedure  was  repeated.  Gear's  method  of  integration  was  employed  to 
estimate  the  new  values  for  particle  mass,  temperature,  velocity  and 
location.  In  the  case  of  a  suspended  drop  the  particle  location  would 
be  fixed  and  the  particle  absolute  velocity  would  be  zero  which  means 
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that  Equations  (3.19)  and  (3.20)  may  be  eliminated  from  the  analysis. 
For  a  free-flowing  drop,  the  support-probe  mass  becomes  zero  and  the 

terms  Y  and  in  would  be  excluded  from  the  analvsis. 

S  s 

3.1.5  Droplet  Combustion  Model 

Drop  transport  in  the  presence  of  envelope  flames  has  been  the 
subject  of  numerous  investigations.  In  the  following,  a  brief  account 
of  envelope-flame  theory  will  be  presented  and  the  results  will  be 
employed  to  indicate  the  effect  of  envelope  flames  on  drop-transport 
rates  within  a  turbulent-diffusion  flame  environment.  Further 
discussion  of  the  properties  of  envelope  flames  can  be  found  in  recent 
articles  [11,48]. 

The  analysis  of  envelope  flames  will  employ  the  same  assumptions 
as  the  drop-evaporation  analysis  of  Section  3.2.3,  except  as  noted  in 
the  following.  It  was  assumed  that  fuel  and  oxidizer  react  in 
stoichiometric  proportions  within  an  infinitely  thin  flame  sheet  which 
completely  surrounds  the  drop.  Reaction  rates  within  the  flame  sheet 
are  assumed  to  be  infinitely  large,  so  that  the  concentrations  of  fuel 
and  oxidant  are  zero  at  the  flame.  In  order  to  reduce  the  complexity 
of  the  transport  equations,  the  same  average  properties  are  employed 
on  both  sides  of  the  flame  sheet  and  the  Lewis  number  is  taken  to  be 
unity.  These  approximations  are  typical  of  most  models  of  envelope 
flames  [11,48].  The  property  assumptions  are  easily  removed  if 
required . 

The  structure  of  the  gas  phase  near  evaporating  and  combusting 
drops,  under  present  assumptions,  is  sketched  in  Figure  21.  The  sane 
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surface  and  ambient  conditions  are  shown  for  the  two  cases.  It  is 
evident  that  the  presence  of  an  envelope  flame  increases  temperature 
and  concentration  gradients  at  the  drop  surface,  enhancing  drop- 
transport  rates.  The  objective  of  the  analysis  is  to  provide  drop 
heat-  and  mass-transfer  rates  given  conditions  at  the  drop  surface  and 
in  the  surrounding  gas.  These  results  can  then  be  incorporated  in 
drop-life-history  computations. 

Following  Law  et  al.  [48],  the  overall  stoichiometry  for  a 
single-component  liquid  fuel  was  taken  to  be 

£flFf]  +  (£f/vf°)[0]  ^  Products  (3.25) 

where  is  the  mass-flux  fraction  the  liquid  phase,  defined  as: 

£  m3L  m  Yfs(|.  +  ^Yc}  (3.26) 


In  these  equations,  f  refers  to  the  liquid-phase  fuel  evaporating  from 
the  drop,  which  also  contains  carbon  particles  and  the  support  probe. 

The  quasisteady-transport  equations  can  be  solved  using  the  film 
theory  approximation,  to  yield: 


=  \Cjl  =  sh  ln  1(1  +  V°  +  V1 


(3.27) 


wht  r  i* 


(1 


Y  )/Y 
as  as 


(3.28) 


(3.29) 


B0  =  Y0  /<Vvf°> 


£f  =  YfS  o  +  M/By  . 


(3.30) 


The  quantity  Y  is  the  mass  fraction  of  nonfuel  species  at  the  drop 
surface.  This  quantity  is  obtained  from  the  equation  of  state  for 
liquid-vapor  equilibrium;  knowing  the  liquid  temperature  and 
composition,  and  the  ambient  pressure.  Basing  the  heat-transfer 
coefficient  on  the  flame  temperature. 


h  =  A (41)  ! fT  -  t  > 

ar  s  f  s  ’ 


(3.31) 


solution  of  the  transport  equations  yields 


i„  Id  +  by)(,  +  bd)! 


(3.32) 


Neglecting  dissociation,  the  flame  temperature  is 


Tf  I  Ts 
T.  •  -  Ts 


(q  0r  /  f 


(3.33) 


where  is  the  heat  released  per  unit  mass  of  oxidizer  consumed  at 
the  flame  front 


Qr  -  ef  Qf/(ef/vf°) 


(3.34) 


If  an  envelope  flame  is  present.  Equations  (3.27)  to  (3.34) 
replace  Equations  (3.14)  and  (3.16)  when  determining  drop  heat-  and 
mass-transfer  rates.  Otherwise,  the  computation  of  drop-life 
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histories  is  the  same  as  for  evaporating  drops.  In  cases  where  the 
surroundings  of  a  drop  contain  oxygen,  but  a  stable  envelope  flame  is 
not  present,  drop  transport  rates  can  be  approximated  by  setting 
Y  =  0  in  Equations  (3.27)  to  (3.34)  which  retrieves  the  formulation 

Uoa 

for  an  evaporating  drop. 


At 

moderate  pressures, 

single-component 

drops  in  a 

fixed 

surroundings  reach  a 

steady 

wet-bulb  temperature 

(aside  from 

minor 

effects 

of  non-unity 

Lewis 

number  which  have  been  ignored  in 

this 

analysis).  In  this  circumstance,  the  rate  of  evaporation  is  given  by 
m"  C 

=  Sh  In  [1  +  (C  (T  -  Ts)  +  Qr  Y0oo)/ifg]  (3.35) 

where  Tg  can  be  approximated  by  the  boiling  temperature  of  the  fuel 
with  little  error  [29] . 

Numerous  comparisons  between  predictions  of  envelope-flame  models 
and  measurements  have  been  reported  [11].  In  general,  models  similar 
to  Equations  (3.27)  to  (3.35)  provide  a  satisfactory  method  for 
correlating  data  once  the  results  are  calibrated  by  selecting  a 
suitable  reference  state  for  computing  average  properties.  The 
present  formulation  has  been  tested  for  pure  fuels  in  combusting  gas 
environments  having  particular  relevance  to  analysis  of  drop— life 
histories  in  spray  flames.  This  includes  measurements  for  supported 
drops  in  laminar  and  turbulent  flames  for  pressures  in  the  range 
0. 1-6.8  MPa  [11]. 
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3.2  Carbon-Particle  Life-History  Model 

Once  all  the  liquid  has  evaporated,  the  carbon  agglomerate  begins 
to  heat-up  from  the  liquid  wet-bulb  temperature.  The  objective  of 
this  portion  of  the  analysis  is  to  predict  the  agglomerate-life 
history,  i.e. ,  the  variation  of  agglomerate  temperature,  mass, 
diameter,  velocity  and  displacement  as  a  function  of  time.  The 
agglomerate  reaction  analysis  follows  the  general  approach  developed 
during  recent  studies  of  carbon-particle  combustion  [85,38-90] ,  but 
differs  in  some  details.  An  illustration  of  the  particle-combustion 
model  is  presented  in  Figure  20. 

3.2.1  Assumptions 

Assumptions  1-3  were  adopted  from  the  liquid-gasification 
analysis.  Additional  assumptions  for  this  portion  of  the  analysis  are 
as  follows : 

4.  The  particle  was  pure  carbon  having  a  uniform  temperature  at 
each  instant  of  time.  The  gas  mixture  in  the  flow  was 
approximated  by  C>2 ,  N2,  C02,  CO,  H20,  H2  and  OH. 

5.  Carbon  reaction  was  limited  to  the  apparent  surface  of  the 

particle.  The  effects  of  pores  (causing  reaction  within  the 
particle).  intrinsic  variations  of  carbon-black  reactivity 
and  catalysts  were  treated  by  introducing  empirical 

area/react ivity  multiplication  factors  of  the  basic  carbon 
reaction  rate. 

6.  Agglomerates  were  observed  to  react  for  equivalence  ratios 
both  greater  and  less  than  unity,  with  the  former  conditions 


involving  low  oxygen  concentrations.  Two  surface  reaction 
aodels  were  examined,  both  of  which  attempt  to  predict 
reaction  rates  over  the  full  range  of  equivalence  ratios 
where  agglomerate  reaction  rates  are  significant.  The  first 
method  employed  an  approach  suggested  by  Neoh  et  al.  [121] 
where  carbon  reaction  with  OH  is  assumed  to  dominate  the 
high-equivalence  ratio  region  while  reaction  with  0^  becomes 
more  significant  at  low-equivalence  ratios.  The  second 
procedure  adopted  the  approach  used  by  Libby  and  Blake 
[88,89]  for  reaction  of  carbon  with  oxygen  and  carbon 
dioxide,  and  water  vapor. 

7.  The  gas  phase  was  assumed  to  be  in  thermodynamic  equilibrium 
at  the  particle  surface. 

8.  As  before,  the  gas  phase  was  assumed  to  be  transparent  to 

radiation;  however,  since  agglomerates  reach  high 

temperatures,  their  radiation  to  solid  surfaces  surrounding 
the  flow  was  considered.  Particle  to  particle  interaction, 
however  was  ignored.  A  unity  value  was  used  for  the 

emissivity  since  the  emissivity  of  carbon  is  fairly  large. 
In  addition,  the  porous  structure  of  the  particle  would 
enable  the  emissivity  to  approach  the  black-body  emissivity 
value  of  unity. 

9.  The  total  gas  pressure  was  constant  throughout  the  region 
surrounding  Che  carbon  particle.  Faeth  [82]  has  shown  that 
this  assumption  is  valid  except  for  particles  much  smaller 
than  the  ones  studied  during  the  present  investigation. 


Ill 


These  assumptions  are  generally  similar  to  those  employed  for 
transport  to  particles  during  combustion  processes  [11,85-39].  The 
main  distinctions  between  the  present  analysis  and  that  of  Libby  and 
Blake  [88,89]  for  carbon-particle  reaction  involves  consideration  of 
different  reaction  mechanisms,  allowance  for  convection  effects,  and 
the  use  of  empirical  area/reactivity  multiplication  factors. 


3.2.2  Conservation  Equations 

The  governing  equations  of  conservation  of  mass,  conservation 
of  energy  and  conservation  of  chemical  elements  are  presented  in  this 
section.  These  equations  were  obtained  from  the  general  equations 
presented  by  Williams  [25]  and  are  applicable  in  all  regions 
surrounding  the  particle.  Employing  spherical  symmetry,  these 
equations  are: 


Conservation  of  mass,  as  given  by  Equation  (3.2)  or  as: 


r2v  =  (Pv)  r  2  =  &"r  2  =  constant 
p  P  P  P 


(3.36) 


which  indicates  that  the  radial  mass  flow  rate  ner  unit  solid  a  is  a 

P 

constant. 


Conservation  of  energy: 


r  di 
p  p  dr 


—  <L_(  r2  di) 

C  dr  dr 
P 


(3.37) 


where 


1X2 


Co<T  ’  Tref)  +  £ 


i=l 


ui 


ref  i 


(3.38) 


Conservation  of  elements: 


2ft”  4Y.- 
P  P 


r  fi-  41L-  PD  4-(r2 


dr 


dr 


1=1,2. 


(3.39) 


The  conservation  of  radial  momentum  for  the  gas  phase  yields  the 
pressure  distribution  external  to  the  particle  and  is,  therefore,  not 
essential  for  the  present  study. 

Assuming  no  gases  are  absorbed  by  the  carbon  particle  there  is  no 
net  mass  flux  of  the  elements  oxygen,  nitrogen  and  hydrogen  at  the 

particle  surface.  Furthermore,  the  net  mass  flux  of  carbon  at  the 

surface  is  equal  to  m^.  Therefore,  the  surface  boundary  conditions 

for  Equations  (3.37)  and  (3.39)  are: 


r  =  r  ;  i  -  i 
P  P 


PD(£Ii)  ,  _  (1  .  y  ) 

dr  p  P  lp' 


(3.40) 


dYi 

PD( — -)  =  y.  i=2  3  4 

dr  P  P  J-P  ’  ’  ’ 


where  carbon  is  denoted  as  species  1.  The  ambient  boundary  conditions 
are : 


r  =  .  i  ~  igp  ,  T  =  Tq 


jj  *  *  3  ^  > 


r  =  r 


1-1.2, 3. 4  . 


(3.41) 
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Equations  (3.41)  treat  the  general  case  where  the  film  radii  for  mass 
and  heat  transfer  are  not  the  same.  These  radii  are  found  from  the 
usual  correlations  of  film  theory,  and  are  given  in  Equations  (3.8) 
and  (3.9). 

The  solution  of  Equation  (3.39),  subject  to  the  boundary 
conditions  of  Equations  (3.40)  and  (3.41)  yields: 

Yt  =  1  -  (1  -  Y1<B)  exp(-2K/Sh)  (3.42) 

P 

Yj  =  YicD  exp(-2K/Sh)  ,  i=2.3,4  (3.43) 


K  =  mj;  fp/pD 


(3.44) 


Employing  the  sign  convention  th~.t  convection  beat  transfer  from 
the  surface  of  the  particle  is  positive,  we  have: 


-X(-H~> 

dr  p 


^  <di\ 

Vdrp 


(3.45) 


where  the  second  part  of  Equation  (3.45)  follows  by  neglecting  the 
slight  departure  of  the  Lewis  number  from  unity.  The  solution  of 
Equation  (3.37)  yields  the  convective  heat  flux,  as  follows: 


9c  Cp 


AK/rr 


Ti  -  ioo)  1  ~  exp(TK7LeNuT 
P 


(3.46) 


3.2.3  Gas-Phase  transport 

This  portion  of  tne  analysis  generally  follows  Libby  and  Blake 
188.893.  and  Will  only  be  briefly  described.  Since  binary 
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dirfusivities  are  equal,  mass  transport  is  conveniently  expressed  in 
terms  of  element  mass  fractions.  Four  elements,  C,  0,  S  and  H,  and 
seven  species,  02,  S2.  C02,  CO.  H2O.  H2  and  OH,  aDpear  in  the 

analysis-  Mote  that  this  analysis  follows  Ubhayakar  and  Williams 
[85],  Libby  and  Blake  [88,89]  and  others  in  neglecting  atomic  oxygen 
on  the  grounds  that  below  4000  K  its  concentration  will  be  small, 
insufficient  to  contribute  significantly  to  the  carbon  oxidation. 

The  element  and  species  mass  fraction  are  related  as  follows; 

7 

Yi  "  Pij  *2,3,4  (3.47) 

where  U._  is  the  mass  fraction  of  the  i'th  element  in  the  j'th 
species.  The  utilization  of  the  above  element  mass  fractions  greatly 
simplified  the  theoretical  analysis  since  elements  are  conserved 
through  the  particle  boundary  layer. 

Given  the  at  a  particular  location.  Equation  (3.47)  provides 
four  equations  to  determine  the  seven  unknown  Yj.  Three  additional 
equsticns  were  obtained  from  the  assumption  of  thermodynamic 
equilibrium.  The  equilibrium  relations  employed  in  the  analysis  were: 

V  5-1  h2  +  1/2  o2 

Kp2 

H20  -  1/2  Hi  +  OH 

K  .> 

CO 2  tJ  CO  +  1/2  0, 


(3.48) 
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where  the  equilibrium  constants  for  the  above  reactions  are  functions 
of  temperature  and  given  as: 


YH.  <Y0  > 


1/2 


,  .  ,  -Jl.  -  2 _  pl/2 

pl  yh,o 


Kn2  =  (^2)I/2  YoH 
v  - y — - 

H-,0 


pl/2 


p,/2 


'p3 


C02 


(2.49) 


Where  the  equilibrium  constants  K  .  are  taken  from  the  JANAF 

pi 

Thermochemical  Tables  [159]. 


3.2.4  Surface  Reactions 

There  are  numerous  uncertainties  in  modeling  the  reaction  of  a 
porous  carbon  agglomerate  in  a  flame-  The  flame  environment  contains 
several  species  which  are  potential  oxidants  of  carbon,  e.g. ,  Ho0. 

C. 

C0„,  0,  C>2  and  OH,  and  the  reaction  mechanism  is  not  well  established 

L 12 1 1  -  Since  the  material  is  porous,  reaction  is  not  limited  to  the 
apparent  outer  surface  of  the  agglomerate  [160].  Carbon  substances 

exhibit  intrinsic  variations  in  reactivity  due  to  their  surface 

structure  [160].  Finally,  the  presence  of  a  catalyst  was  found  to 

influence  reaction  rates. 

The  effect  of  various  gaseous  reactants  was  treated  by 
considering  a  mechanism  recently  proposed  by  Seoh  et  ai.  [121]  as 
well  as  an  extension  of  an  earlier  approach  employed  by  Lihbv  and 
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Blake  [38,89].  The  effects  of  pores,  surface  reactivity  and  catalysts 
were  treated  by  identifying  empirical  area-reactivity  multiplication 
factors,  selected  to  best  match  the  present  measurements.  (It  was 
found  that  a  fixed  value,  for  a  given  carbon  black  and  reaction,  could 
correlate  the  measurements  over  the  range  of  the  data.) 

Neoh  et  al.  [121]  concluded  that  OH  is  the  dominant  carbon 
oxidant  under  fuel-rich  conditions;  that  0  is  of  secondary  importance 
at  temperatures  below  20C0  K;  and  that  02  becomes  a  significant  factor 
for  fuel-lean  conditions,  particularly  at  lower  temperatures.  Present 
test  conditions  correspond  to  these  circumstances ,  therefore,  carbon 
reaction  with  both  OH  and  02  was  considered,  while  ignoring  the 
remaining  species.  Neoh  et  al.  [121]  found  that  reaction  of  carbon 
with  OH  could  be  represented  by  assuming  a  constant  collision 
efficiency,  yielding  the  following  reaction  rate  expression 

Rl  “  KrlpOH  .  (3.50) 

The  results  of  Nagle  and  Strickland-Cons table  [109]  were  employed  to 
determine  the  rate  of  reaction  of  carbon  with  These  rate 

expressions  are  also  in  reasonable  agreement  with  later  measurements 
by  Park  and  Appleton  [120].  The  carbon  reaction  rate  is  given  by 

R2  =  Kr2p02X/O  +  Kr3P0?)  +  Kr4P0?(l  -  X)  (3.51) 

where 

X  =  f  1  +  *r5^^ripO?^ 


(3.52) 
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The  measurements  of  Neoh  et  al .  (121)  suggest  that  the  larger  of  or 
should  represent  the  reaction  rate  at  any  condition.  Therefore, 
the  dimensionless  mass  burning  rate  was  determined  from 


where  the  a^  appearing  in  Equation 
reactivity  multiplication  factors. 

The  carbon-reaction  mechanism 
considers  reaction  with  0^  and 
expressions : 


(3.53) 
empirical  area- 

used  by  Libby  and  Blake  [88] 
CC>2  employing  the  following 


K.  =  (rp/pD)  max  [ajRj,  a2R2) 

(3.53)  are  the 


R3  = 


K  ✓ 
r6  02 


(3.54) 


R 


4 


K 


r7fco2  _ 


(3.55) 


This  approach  was  extended  during  the  present  study  to  include  the 
reaction  of  carbon  with  H^O,  by  employing  the  results  of  Johnstone  et 
al.  [103].  The  reaction  rate  expression  for  this  case  is 


R5  “  Kr8PH20/(1  +  Kr9PH2  +  *rl()l>H  0>  •  0-56) 

For  this  second  approach,  the  total  reaction  rate  of  carbon  was 
obtained  as  the  sum  of  the  rates  for  0^,  CO2  and  R2O,  ignoring 
potential  interactions  between  reactants.  This  yields: 

5 

K  =  ( r p/PD)  l  a.R. 


(3.57) 
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The  specific  reaction  rate  parameters  used  in  the  computations 
are  summarized  in  Table  5,  where  the  are  assumed  to  have  the 

following  general  form 

ni 

Kri  =  AiT  exp(-ei/RT)  .  (3-58) 

In  all  cases,  P^  was  found  from  the  mass  fractions  at  the  particle 

surface  as  follows: 

Pf  -  M  YiP/Mi  .  (3.59) 

The  area-reactivity  multiplication  factors  which  provided  the 
best  match  of  the  present  data  are  summarized  in  Table  6  for  the 
turbulent-diffusion  flame  tests.  Neoh  et  al.  [121]  estimate  collision 
efficiencies  for  OH  in  the  range  0,13-0.28,  the  latter  value  was 
employed  for  the  present  estimate  of  a^.  The  area-reactivity 
multiplication  factors  for  R^  and  could  not  be  separated  for 
turbulent-diffusion  flame  test  conditions,  since  the  relative 
proportions  of  CO^  and  H^O  were  roughly  the  same  throughout  the  flame: 
therefore,  the  values  were  taken  to  be  identical  for  both  reactions. 
As  might  be  expected,  the  presence  of  pores  in  the  agglomerate  result 

in  values  greater  than  unity  since  the  actual  reaction  area  is 

greater  than  the  apparent  surface  area.  The  a^  for  a  catalyzed  slurry 
were  only  evaluated  for  the  second  reaction  approximation.  The  use  of 
a  catalyst  increases  the  intrinsic  surface  reactivity  of  the  carbon, 
yielding  higher  a^  than  for  the  noncatalyzed  agglomerate. 
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Table  5-.  Summary  of  Reaction  Rate  Parameters 


i 

A 

n 

E  (kcal/gmol) 

ia 

1/2  2 

361  kg  K  /ms  atm 

-1/2 

0 

2 

2 

2400  kg/m  s  atm 

0 

30.0 

3 

21.3  atm-1 

0 

-  4.1 

4 

2 

0.535  kg/m  s  atm 

0 

15.2 

5 

18.1  x  10^  kg/m2s 

0 

97.0 

6 

2 

87100  kg/m  s  atm 

0 

35.8 

7 

2 

2470  kg/m  s  atm 

0 

41.9 

eQ 

00 

15.15  x  10  3  kg/m2s  atm 

0 

32.7 

9 

9.42  x  lO'11  atm-1 

0 

-60.8 

10 

7.07  x  10~16  atm"1 

0 

-79.3 

Assuming  a  collision  efficiency  of  0.28  [121]. 

jj  2 

Assuming  a  surface  area  of  reaction  of  1.15  m  /g  of  carbon, 
which  is  the  average  surface  area  over  the  period  of  reaction 
[ 103] . 


Table  6.  Summary  of  Area-Reactivity  Factors 
for  Turbulent  Diffusion  Flame  Tests 


Reaction  Noncatalyzed  Catalyzed 


Neoh,  et  al. ,  £121]  reaction  approximation: 

1  54.3 

2  758. 

Extended  Libby  and  Blake  [88]  reaction  approximation: 


58.2 

70.6 

112.9 

140.0 

112.9 

140.0 

5 
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3.2.5  Particle-Life  History 

The  results  of  the  previous  sections  will  now  be  applied  to 
determine  the  particle-life  history  of  a  cold  particle  suddenly 
immersed  in  a  hot-oxidizing  ambient.  To  simplify  the  analysis,  the 
particle  will  be  assumed  to  have  a  uniform  temperature. 

Assuming  that  the  density  of  the  particle  is  a  constant, 
conservation  of  particle  mass  yields: 


Jl£  -  -PDK/Pcrp 


(3.60) 


Conservation  of  energy  at  the  particle  surface  yields: 


<p^-rMt'p(1p  -  V  +  <  + 


'P  C  . 
c  pc  p 


(3.61) 


where  qr  is  the  particle  surface  heat  flux  for  radiation  to  the 
enclosure  surrounding  the  flow: 

q"  *  O  e  (t£  -  tJ)  .  (3.62) 


Tne  particle  velocity  and  displacement  were  obtained  by  integrating 
Equations  (3.19)  and  (3.20)  after  replacing  the  average  droplet 
density  with  the  carbon  density. 

The  initial  conditions  for  Equations  (3.60)  and  (3.61)  were  based 
on  conditions  at  the  end  of  the  liquid-vaporization  period.  At  this 
point,  all  the  carbon  initially  in  the  slurry  was  present  in  the 
agglomerate,  the  agglomerate  density  is  known  [2]  and  the  particle  was 
at  the  wet-bulb  temperature  for  liquid  evaporation. 
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As  mentioned  earlier,  the  model  assumes  a  uniform  temperature  for 
the  carbon  particle  as  illustrated  in  Figure  22.  The  temperature 
profile  of  the  gas  phase  decreases  from  the  particle  temperature  at 
the  particle  surface;  to  the  ambient  temperature,  at  a  radius 
predicted  by  the  film  theory.  Figure  22  also  illustrates  that  the 
oxidant  mass  fraction  decreases  and  .the  product  species  mass  fraction 
exponentially  increases  from  their  ambient  values  to  their  values  at 
the  particle  surface. 

3.2.6  Calculation  Procedure 

Initially  the  time,  particle  radius,  temperature,  velocity  and 
displacement  were  set  equal  to  their  values  at  the  end  of  the  liquid- 
gasification  period.  The  values  for  the  mass  fraction  at  the  particle 
surface  are  initially  set  equal  to  the  ambient  values.  The  transport 
and  thermophysical  properties  of  the  gas  phase  surrounding  the 
particle  were  evaluated,  as  discussed  in  Appendix  A,  considering  the 
mass  fraction  and  temperature  of  the  particle  surface  and  of  the 
surroundings. 

The  equilibrium  constants:  Kpl ,  Kp2  and  Kp3,  from  Equations 
(3.48)  and  (3.49),  were  evaluated  using  the  particle  temperatures, 
ambient  pressures  and  regression  equations  curve  fitting  the  chemical 
equilibrium  data  of  the  JANAF  tables  [159J.  Values  for  the  kinetic 
rates  given  in  Equation  (3.58)  were  then  determined  at  the  particle 
temperature.  The  film  parameters  were  determined  next  from  Equations 
(3.8)-(3. 10).  The  non-dimensional  carbon  mass-loss  parameter  was  then 
calculated  by  employing  Equations  (3. 50)-( 3. 53)  or  Equations 
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Figure  22.  Carbon-agglomerate  reaction  model. 
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( 3. 54)-(3 . 57 ) ,  depending  upon  whether  the  OH-O2  mechanism  or  the 
O2-CO2-H2O  mechanism  was  chosen.  Next,  the  gas-phase  mass  fractions 
(O2,  CO2,  CO,  H2,  OH,  H2O  and  N2)  at  the  particle  surface  were 
calculated  by  simultaneously  solving  Equations  (3.39)  and  (3.47)  by 
either  the  Newton-Raphson  method,  or  the  interval  bisection  method. 
The  Newton-Raphson  method  was  normally  used,  while  the  interval 
bisection  method  was  employed  for  the  more  unstable  cases. 

The  properties  of  the  gas  phase  were  then  reevaluated  with  the 
updated  mass  fractions.  The  procedure  was  continued  until  convergence 
of  the  mass  fractions  was  obtained.  At  this  point.  Equations  (3.20), 
(3.60)  and  (3.61)  were  solved  for  the  values  of  particle  displacement, 
radius  and  temperature  at  the  next  time  increment,  using  Gear's  method 
of  integration.  Particle  velocity  was  obtained  from  experimental 
measurements.  The  entire  procedure  was  continued  until  the  particle 
was  consumed  or  extinguished.  If  the  particle  was  suspended,  the 
particle  location  is  fixed  and  Equations  (3.19)  and  (3.20)  may  be 
dropped  from  the  analysis. 


CHAPTER  IV 


RESULTS  AND  DISCUSSION 

The  main  objective  of  che  present  study  was  to  compare 
measurements  and  predictions  fc r  the  heat-up  and  combustion  processes 
of  a  slurry  droplet  und-nr  conditions  which  simulate  the  environment  in 
a  combustion  chamber.  The  theoretical  model  discussed  in  Chapter  III 
was  used  to  predict  the  diameter,  temperature,  mass,  velocity  and 
position  of  a  slurry  droplet  burning  in  both  a  turbulent-diffusion 
flame  and  a  lamina r-premixed  flat  flame. 

4 . 1  Slurry  Droplet  Combustion  in  the  Turbulent-Diffusion  Flame 

For  the.  turbulent  diffusion  flame  tests,  a  slurry  droplet  was 
supported  on  a  probe  in  the  flame.  The  droplets  were  observed  to  have 
various  degrees  of  reaction  in  the  different  flame  regions.  The  drop- 
life-hi story  and  the  particle-life-history  models  were  employed  tc 
predict  the  heat-up  and  combustion  process  of  the  slurry  droplet. 
Each  of  these  results  are  discussed  in  the  following  sections. 

4.1.1  Flame  Regions 

Figure  23  is  an  illustration  of  the  variation  in  mean 
properties  along  the  centerline  of  the  turbulent  flame  fueled  with 
gaseous  propane.  In  the  region  near  the  nozzle  exit,  gas  velocities 
and  the  concentrations  of  propane  are  relatively  high,  while  the 
temperature  is  low  and  oxygen  is  absent.  Moving  downstream,  the  gas 
velocity  and  concentration  of  propane  decrease  monotonically.  The 


MASS  FRACTION 


10.0 


NONCATALYZED  A' a  NC 

SLURRY  1 _ E _ 1  i  FC  1  E  i 


CATALYZED 

SLURRY 


£ _ ff  *  ft  E  , 


Figure  23.  Droplet  environment  diagram  modified  to 

show  the  various  droplet  burning  regions. 
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temperature  increases  at  first,  reaching  a  maximum  at  the  flane  tip, 
where  the  fuel  has  disappeared  and  combustion  product  concentrations 
are  highest.  Beyond  the  tip  of  the  flame  the  flow  decays  with 
increasing  concentrations  of  oxygen  and  nitrogen  and  decreasing 
temperatures  and  combustion  product  concentrations. 

The  structure  of  the  burner  flame  used  in  this  study  is  typical 
of  other  gas  and  spray  fueled  diffusion  flames  [11,149-151,157].  The 
region  near  the  injector  is  highly  fuel  rich  and  has  a  relatively  low 
temperature.  Significant  concentrations  of  oxygen  and  high 
temperature  levels  only  appear  as  the  turbulent  reaction  zone  is 
approached. 

Drop  behavior  varied  substantially  with  position  in  the  flame. 
Table  7  is  a  description  of  the  various  types  of  behavior  that  were 
observed:  fragmentation  (F),  noncombusting  (NC),  full  combustion  (FC), 
glowing  (G),  and  evaporatioi  (E>.  The  regimes  are  listed  in  the  order 
in  which  they  appear  with  increasing  distance  from  the  injector. 
Table  3  is  a  summary  of  the  portion  of  the  flame  where  the  various 
regions  were  observed,  for  both  the  catalyzed  and  noncatalyzed 
slurries. 

The  various  drop  burning  regimes  are  also  indicated  on  the  lower 
part  of  Figure  23.  The  observations  in  each  of  these  regions  will  be 


considered  in  turn. 


In  order  of  increasing  distance  from  the  Injector. 
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i.  1 . 1 . 1  Fragmentation  Revlon 

Th^  fragmentation  region  is  closest  to  the  injector, 
involving  relatively  high  gas  velocities  and  propane  concentrations, 
noderate  combustion  product  concentrations,  virtually  no  oxygen  and 
temperatures  from  300  to  roughly  1600  K.  In  this  region,  the  liquid 
fuel  evaporated  with  no  envelope  flame  present.  A  luminous  wake  was 
visible  behind  the  drop  as  the  liquid  evaporated.  This  behavior  is 
typical  of  a  liquid  fuel  evaporating  in  a  high-temperature  gas  in  the 
absence  of  oxygen  [11,41].  The  liquid  drop  primarily  evaporates 
without  combustion.  As  the  fuel  g>ases  become  heated  they  decomposed 
forming  soot.  When  the  temperature  is  high  enough,  the  soot  glows 
yielding  a  diffuse-luairswus  wake. 

When  the  bulk  of  the  liquid  had  evaporated,  most  of  the  carbon 
agglomerate  was  blown  off  the  probe  in  three  to  five  large  flakes. 
The  mechanical  breakdown  of  the  agglomerate  results  from  the 
relatively  high  gas  velocities  in  this  region.  Figure  24  is  a  SE>! 
photograph  of  the  residue  remaining  of  the  probe  after  exposure  of  a 
catalyzed-slurry  drop  to  the  fragmentation  region  of  the  flame.  The 
appearance  of  the  residue  from  the  noncatalyzed  slurry  was  similar. 
The  surface  is  very  irregular  with  waves  frozen  into  the  structure  by 
the  high  gas  velocity.  The  surface  structure  is  similar  to  that 
observed  when  these  drops  are  burned  in  air,  c.f..  Figure  3.  At 
higher  magnifications,  sub-elements  having  a  size  comparable  to  the 
ultimate  carbon  particle  size,  0.3  m,  were  observed. 
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4. 1.1.2  Noncombusting  Region 

The  noncombusting  region  was  relatively  narrow  for  both 
fuels.  It  was  located  on  the  fuel-rich  side  of  the  reaction  zone; 
where  velocities  and  oxygen  concentrations  were  relatively  low,  while 
temperatures  and  product  concentrations  approached  their  maximum 
values.  The  liquid  fuel  evaporated  without  an  envelope  flame,  similar 
to  behavior  in  the  fragmentation  region.  However,  there  was  little 
flaking  of  the  residue  and  greater  amounts  of  carbon  remained  on  the 
probe.  After  the  liquid  had  evaporated,  a  luminous  wake  of  reduced 
intensity  was  still  observed.  This  suggests  continued  shedding  of 
small  carbon  particles  which  eventually  glowed  as  they  heated  (and 
began  to  react)  in  the  flow. 


Figure  25  is  a  SEM  photograph  of  the  surface  of  the  residue  from 


the  noncombusting  region.  The  residue  from  both  nontcatalyzed  and 
catalyzed  slurries  was  similar.  The  structure  is  similar  to  the 


structure  seen  in  the  fragmentation  region.  The  major  difference 
between  the  two  regimes  is  that  the  lower  velocities  do  not  cause  the 


residue  to  flake.  Rates  of  carbon  reaction  are  also  probably  somewhat 
higher,  although  the  rates  are  still  relatively  low.  The  support  of 
particles  is  weakened  due  to  reaction,  which  caused  small  portions  to 
flake  off. 


4. 1.1. 3  Full-Combustion  Region 

The  full-combustion  region  has  low  gas  velocities,  and 
relatively  high  temperatures  and  high  concentrations  of  oxygen  and 
combustion  products.  An  envelope  flame  was  observed  around  the  drop 
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as  the  liquid  evaporated.  Similar  to  combustion  in  air,  a  dark  period 
was  observed  between  the  time  when  the  envelope  flame  was  extinguished 
and  the  residue  began  to  glow.  The  dark  period  involves  transient 
heating  of  the  surface  from  the  low  temperature  levels  characteristic 
of  liquid  evaporation  to  the  high  temperatures  characteristic  of 
carbon  reaction.  In  contrast  to  combustion  in  air,  both  the 
noncatalyzed-  and  catalyzed-slurry  residues  exhibited  glowing  in  the 
burner  flame.  The  glowing  persisted  until  very  little  carbon  residue 
remained. 

Figure  26  shows  two  SEM  photographs  of  the  residue  from  a 
catalyzed-slurry  drop  in  the  full-combustion  region.  The  residue  from 
the  noncatalyzed  slurry  was  similar.  The  structure  is  thin  and  very 
porous.  In  this  region,  the  residue  was  a  cenosphere;  which  is 
probably  the  result  of  convective  and  radiative  heat  losses,  which 
tend  to  quench  the  reaction  at  the  outer  surface.  The  shielding 
effect  of  the  outer  structure  allows  higher-temperature  levels  in  the 
interior  of  the  residue  and  the  reaction  completes  the  gasification  of 
the  solid.  The  greater  extent  of  reaction  on  the  outer  surface 
results  in  a  structure  that  is  less  related  to  the  original  slurry 
particles  than  was  the  case  for  the  region  nearer  to  the  injector. 

4. 1.1.4  Glowing  Region 

The  glowing  region  was  only  observed  for  the  catalyzed 
slurry,  and  is  relatively  narrow.  Oxygen  concentrations  are  high  in 
this  region  (mass  fractions  in  the  range  0.1j-0.20)  and  temperatures 
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are  relatively  low  (800-1000  K).  The  liquid  fuel  evaporated  with  no 
envelope  flame  present  and  no  luminous  wake.  Some  time  after  the 
liquid  had  evaporated,  2  to  15  s  for  present  test  conditions,  the 
residue  began  to  glow.  Glowing  persisted  until  little  carbon  was  left 
on  the  probe. 

The  surface  structure  of  the  residue  from  this  region  was  similar 
to  the  full-combustion  region.  SEM  photographs  of  the  surface  appear 
in  Figure  27.  The  structure  is  very  porous. 

4. 1.1.5  Evaporation  Region 

The  evaporation  region  was  farthest  from  the  injector,  where 
the  flow  decays  to  approach  the  properties  of  room  air.  In  this 
region,  the  liquid  evaporated  leaving  a  carbon  agglomerate-  There  was 
little  evidence  of  subsequent  reaction  of  the  agglomerate. 

SEM  photographs  of  the  residue  In  this  region  appear  in  Figure 
28.  The  surface  appears  smooth  in  some  regions  and  coarse  in  others 
(due  to  fracture  of  the  sample).  The  appearance  is  generally  similar 
to  evaporation  in  air  with  little  reaction  of  the  agglomerate,  c.  f. , 
Figures  4  and  5. 

4. 1.1.6  Effect  of  Catalyst 

The  appearance  of  the  residue  from  both  slurry  fuels  was 
similar  in  each  region.  The  main  effect  of  the  catalyst  was  to  extend 
the  region  where  reaction  of  the  residue  could  be  sustained.  The 
noncatalyzed  residue  reacted  in  the  range  x/d  =  159-383  while  the 
catalyzed  residue  reacted  in  the  range  x/d  =  159-500.  This 


SE>1  photograph  of  the  carbon  residue  from  a  catalyzed 
slurry  droplet  rro.-n  a  glowing  region,  x/d  =  439. 
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substantially  increases  the  residence  time  where  reaction  of  residue 
occurs  in  Che  flow.  For  example,  if  the  particles  are  small  enough  to 
move  with  the  gas  velocity 


fxe 

tr  *  j  dx/u  (4.1) 

where  x.  and  xg  are  the  beginning  and  end  of  the  reaction  region,  tr 
is  the  residence  time  in  the  reaction  zone,  x  is  distance  from  the 
injector  and  u  is  the  gas  velocity.  Completing  the  integration  of 
Equation  (4.1)  for  present  burner  conditions  yields  tr  =  41  ms  for  the 
noncatalyzed  residue  and  tr  =»  76  as  for  the  catalyzed  slurry — almost 
twice  as  long.  Naturally,  within  an  actual  combustion  chamber, 
conditions  will  be  different.  It  seems  likely,  however,  that  the 
catalyzed-slurry  residue  will  react  over  a  greater  portion  of  the 
flow. 

4.1.2  Liquid  Evaporation 

Results  during  the  steady-evaporation  period  of  the  drop  were 
employed  for  initial  evaluation  of  the  liquid-gasification  theory. 
During  steady  evaporation,  with  fixed-ambient  conditions,  the  drop  was 
stabilized  at  its  wet-bulb  temperature  where  virtually  all  the  heat 
transferred  to  the  drop  is  utilized  for  the  heat  of  vaporization  of 
the  evaporating  liquid.  Predicted  and  measured  evaporation  rate 
constants,  defined  as 


K 


e 


(4.2) 
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were  compared.  Since  Kg  varies  with  drop  diameter  due  to  convection, 
values  were  compared  at  various  drop  sizes. 

Predicted  and  measured  evaporation-rate  constants  for  the  steady 
evaporation  period,  are  illustrated  in  Figure  29.  The  evaporation- 
rate  constant  is  plotted  as  a  function  of  particle  diameter  with  the 
equivalence  ratio  of  the  flame  taken  as  a  parameter.  The  predictions 
employed  a  =  0.9,  which  is  the  value  that  has  generally  yielded  the 
best  correlation  between  predictions  and  measurements  during  past  work 
in  this  laboratory  [11.150,151,155,156]. 

Predictions  and  measurements  in  Figure  29  are  in  good  agreement, 
except  for  results  at  equivalence  ratios  of  0.761  and  0.870  where  the 
predictions  underestimate  the  measurements.  Envelope  flames  were 
observed  at  these  conditions  which  accounts  for  the  enhancement  of 
evaporation  rites  over  present  preaictions  which  ignore  envelope 
flames.  Farther  downstream,  where  the  equivalence  ratio  was  lower,  an 
envelope  flame  would  cause  an  even  grater  enhancement  of  evaporation 
rates  over  present  predictions;  however,  the  drops  could  not  sustain 
stable  envelope  flames  at  these  positions.  At  higher  equivalence 
ratios,  =  0.939  and  1.350,  mean  oxygen  concentrations  are  finite 
around  the  drop  (c.f..  Figure  12)  and  envelope  flames  may  also  have 
been  present  although  they  could  not  be  observed  due  to  flame 
luminosity.  This  would  have  little  consequence,  however,  since  the 
presence  of  envelope  flames  near  equivalence  ratios  of  unity  does  not 
have  a  large  influence  on  drop  evaporation  rates  [111.  The  model 
could  easily  be  extended  to  accommodate  envelope  flames.  However, 
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this  was  not  done  here  since  liquid  evaporation  in  combusting  sprays 
is  generally  confined  to  regions  where  envelope  flames  are  unlikely, 
even  though  they  were  observed  when  drops  were  supported  in  a 
diffusion  flame  [11,150,151,155,156]. 

To  illustrate  the  effect  of  droplet  combustion  a  simple 
calculation  was  performed  for  a  435  m  initial  diameter  combusting 
liquid  JP-10  droplet  over  the  range  of  fuel  equivalence  ratios  of  the 
bulk  gas.  These  results  are  compared  to  the  theoretical  prediction 
for  an  evaporating  droplet,  and  the  experimental  measurements  in 
Figure  30  for  the  same  droplet  size  and  conditions.  For  equivalence 
ratios  greater  than  unity,  the  ambient  oxygen  concentration  is  small 
and  both  theories  yield  the  same  burning-rate  constant,  which  is  in 
good  agreement  with  the  data.  As  the  equivalence  ratio  decreases 
below  unity,  the  theories  begin  to  diverge.  In  the  region  where 
envelope  flames  were  observed  the  combusting  drop  model  yields  bast 

v. 

agreement  with  the  data.  As  the  equivalence  ratio  continues  to 
decrease,  however,  envelope  flames  are  no  longer  stable  and  the 
measurements  shift  to  the  predictions  for  evaporating  drops.  The 
results  indicate  that  the  models  give  reasonable  predictions  for  both 
evaporating  and  combusting  drops  in  sprays,  as  long  as  the  appropriate 
branch  of  the  solution  can  be  selected  for  equivalence  ratios  less 
than  unity. 

4.1.3  Carbon  Combustion 

Equation  (3.10)  was  developed  from  theoretical  and  experimental 
results  for  smooth  spheres  and  drops  and  provided  a  good  correlation 


14 


of  convective  effects  for  drops  during  this  and  earlier  investigations 
{ 7—  ill-  The  use  of  this  equation,  however,  underestimated  transport 
races  of  the  present  agglomerates  for  all  reasonable  selections  of 
physical  properties.  This  behavior  was  attributed  to  the  fact  that 
the  agglomerates  have  an  open  porous  structure  near  the  surface 
providing  roughness  effects  and  allowing  the  flow  to  penetrate  the 
apparent  surface  to  some  extent.  It  was  found  that  this  enhancement 
of  convective  transport  rates  could  be  accommodated  in  the  model  by 
multiplying  the  Xusselt  and  Sherwood  numbers  given  by  Equation  (3.10) 
by  a  constant  convective  enhancement  factor  of  6,7  throughout  tr.e 
carbon-combustion  analysis  for  the  turbulent-diffusion  flame. 

Predicted  and  measured  burning  rates  For  noneataiyzed 
agglomerates  are  illustrated  in  Figures  31  and  32.  Burning  rate  is 
plotted  as  a  function  of  agglomerate  diameter  with  flame  position 
(equivalence  ratio)  as  a  parameter.  The  data  for  the  turbulent- 
diffusion  flame  tests  are  tabulated  in  Appendix  B.  Equivalence  ratios 
of  0.870-1.350  are  illustrated  in  Figure  31,  while  lean  conditions, 
corresponding  to  equivalence  ratios  of  0.272-0.761,  are  illustrated  in 
Figure  32.  The  predictions  employ  the  reaction  parameters  summarized 
in  Table  3.  For  this  selection  of  parameters,  the  two  reaction 
approximations  gave  nearly  the  same  results  and  both  methods  can  be 
represented  by  a  single  line.  Clearly,  the  use  of  a  single  set  of 
area-reactivity  multiplication  factors,  a^,  provides  a  good 
representation  of  the  data  over  the  entire  range  of  the  tests.  The 
low  reaction  rates  at  low  equivalence  ratios  were  also  predicted 
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PARTICLE  DIAMETER  (Mm) 


Figure  31.  Quasisteady  burning  rates  for  noncataiyzed 

carbon  agglomerates  at  'nigh  equivalence  ratios. 


reasonably  well,  suggesting  that  the  model  has  some  potential  for 
predicting  quenching  of  the  agglomerate  reaction. 

There  are  two  main  combustion  limits  for  carbon  agglomerates,  the 
reaction-controlled  limit  and  the  diffusion-controlled  limit  [159]. 
For  reaction  control,  reaction  rates  are  low  in  comparison  to 
diffusion  rates,  and  there  is  no  mass  transfer  limitation  on  the 
supply  of  gaseous  reactants  to  the  particle  surface.  This  limit  is 
characterized  by  the  burning  rate  being  relatively  independent  of 
particle  diameter,  e.g.,  the  results  for  x/d  382.5 
in  Figure  32  approach  the  reaction  controlled  limit. 

The  diffusion-controlled  limit  is  characterized  by  high  reaction 
rates  in  comparison  to  diffusion  rates,  leading  to  small 
concentrations  of  reactant  gases  at  the  particle  surface.  In  this 
case,  the  burning  rate  is  largely  determined  by  the  mass-transfer 
properties  of  the  flow  around  the  particle.  The  magnitude  of  the 
burning  rate  is  then  influenced  by  particle  diameter,  as  follows: 


ddP  -  A  -(0.5  to  1.0) 

<rr~  p 


(4.3) 


The  power  of  d^  in  Equation  (4.3)  tends  toward  -1  as  the  Reynolds 
number  of  the  particle  decreases.  From  this  argument,  it  is  apparent 
that  results  in  the  range  255.0  <_  x/d  340,  in  Figures  31  and  32  are 
approaching  the  diffusion-controlled  limit.  The  remainder  of  the 
measurements  involve  both  reaction  and  diffusion  effects  to  some 
extent.  Diffusion  rates  become  large  for  small  particles;  therefore, 
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all  regions  of  the  flow  would  eventually  exhibit  reaction  control  for 
sufficiently  small  particles. 

Figure  33  is  an  illustration  of  predicted  and  measured  burning 
rates  for  the  catalyzed  agglomerate.  These  test  results  were  limited 
to  x/d  =  297.5,  340.0  and  382.5  and  were  analyzed  only  using  the 
extended  Libby  and  Blake  [88,89]  reaction  mechanism.  Reaction  of  the 
agglomerate  is  nearly  diffusion  controlled  at  the  first  two  positions, 
resulting  in  little  change  in  the  burning  rate  when  the  catalyst  was 
used.  At  x/d  *  382.5,  however,  adding  catalyst  increased  the  burning 
rate  20-25%,  since  the  rate  of  gasification  of  carbon  is  nearly 
reaction  controlled  at  this  position.  Naturally,  the  effect  of 
catalyst  would  be  even  greater  at  positions  farther  downstream,  e.g. , 
the  presence  of  catalyst  reduced  the  lean  limit  for  quenching  from  an 
equivalence  ratio  of  0.37  to  0.21. 

Another  property  of  catalyzed  slurries  was  their  tendency  to 
spontaneously  shatter  from  time-to-time.  When  observed,  shattering 
occurred  after  all  the  liquid  had  apparently  evaporated.  This 
behavior  could  be  influenced  by  the  presence  of  the  drop  support,  and 
was  not  consistently  observed.  If  shattering  does  persist  for  small 
unsupported  particles,  however,  this  would  provide  another  (physical) 
mechanism  for  the  improved  combustion  properties  observed  for 
catalyzed  slurries  during  combustion-chamber  tests  [2,3].  Tests  with 
unsupported  particles,  having  sizes  more  typical  of  sprays,  will  be 
required  in  order  to  establish  the  importance  of  this  effect. 
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PARTICLE  DIAMETER  (^m) 


Figure  33.  Quasisteady  burning  rates  for  catalyzed  agglomerates . 
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The  final  step  in  evaluating  the  model  for  the  turbulent- 
diffusion  flame  tests  involved  comparison  of  predicted  and  measured 
slurry  drop-life  histories.  Life  histories  were  computed  for 
noncatalyzed  slurries  over  the  test  range.  The  results  are 
illustrated  in  Figures  34-41  which  covers  equivalence  ratios  of 
0.272-1-350.  Predicted  and  measured  particle  diameters  and 
temperatures  are  plotted  as  a  function  of  time,  with  results 
illustrated  for  both  reaction  approximations.  The  local  temperature 
of  the  flow  is  also  shown  on  each  plot,  for  reference  purposes.  The 
time  scale  of  Figures  34-37  has  been  expanded  in  the  period  where 
liquid  is  present,  since  this  period  is  short  in  comparison  to  the 
total  lifetime  of  the  slurry  drop.  The  liquid  evaporation  period  is 
not  shown  in  Figures  38-41. 

The  two-stage  combustion  process  of  slurry  drops  is  quite  evident 
from  the  results  pictured  in  Figures  34-37.  The  first  stage  involves 
heat-ur  of  the  particle  to  the  wet-bulb  temperature  for  liquid 
gasification.  Some  evaporation  occurs  throughout  this  period; 
however,  the  rate  is  highest  when  the  particle  is  near  the  wet-bulb 
temperature.  Since  JP-10  is  a  pure  hydrocarbon,  its  wet-bulb 

temperature  is  nearly  constant  at  a  fixed  location  in  the  flame 
(having  a  value  somewhat  below  the  boiling  temperature  of  the  liquid) 
(llj.  The  second  stage  involves  heat-up  and  reaction  of  the 
agglomerate,  once  all  the  liquid  has  evaporated.  Initially,  the 
variation  of  agglomerate  diameter  is  relatively  slow  in  this  period, 
until  the  particle  temperature  approaches  the  gas  temperature. 
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Figure  34.  Noncatalyzed  slurry  drop-life  history  at  x/d  =  170.0. 
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35.  Noncatalyzed  slurry  drop-life  history  at  x/d  =  212.5. 
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Figure  36-  Noncatalyzed  slurry  drop-life  history  ac  x/d  -  340.0- 
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Figure  37.  Noncatalyzed  slurry  drop-life  history  at  xl d  =  467.5. 
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Figure  38.  Noncataly^ed  slurry  drop-life  history  at  x/d 
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Figure  39-  Noncatalyzed  slurry  drop-life  hiscory  at  x/d  ~  29,-5 
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Figure  AO.  Noncatalyzed  siurry  drop-life  history  at  x/d 
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Figure  41.  Noncatalyzed  slurry  drop-life  history  at  x/d  =  510.0. 
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Maxinun  parci.de  temperatures  can  be  above  or  below  the  gas 
temperature,  depending  on  the  energy  release  by  reaction,  and  heat 
transport  by  convection  and  radiation.  Near  the  lean  and  rich  limits, 
particle  temperatures  tend  to  be  below  the  flame  temperature  due  to 
radiative  heat  losses.  Intermediate  equivalence  ratios,  e.g.,  0.557 
for  Figure  36,  yield  relatively  high  reaction  rates  and  particle 
temDeratares  exceed  the  flame  tempe rat jre.  The  relative  rates  of 
reaction,  convection,  and  radiation,  however,  vary  with  particle  size. 
Therefore.  there  is  no  fixed  particle  temperature  during  the  latter 
part  of  agglomerate  reaction,  even  at  a  fi„ed  location  in  the  flan°, 
unlike  the  liquid  evaporation  period.  As  the  particle  becomes  small, 
however,  convection  heat  transfer  rates  biccme  large  in  comparison  to 
reaction  and  radiation  effects,  and  the  particle  eventually  approaches 
the  local  gas  temperature.  This  is  not  observed  in  Figures  34-41  due 
to  the  finite-sized  bead  used  to  mount  the  slurry  drop. 

The  comparison  between  the  predictions  and  measurements 
illustrated  in  Figures  34-41  is  quite  good.  This  is  partly  expected, 
since  the  model  was  capable  of  correlating  the  burning-rate 
measurements  for  this  range  of  flame  conditions.  However,  it  is 
encouraging  that  the  model  is  also  capable  of  predicting  heat-up, 
where  reaction  effects  are  small,  as  well  as  final  temperature  levels, 
where  radiation  becomes  important.  There  is  little  difference  in  the 
predictions  using  the  two  reaction  approximations,  similar  to  the 
agglomerate  burning-rate  results. 
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4 . 2  Carbon-Agglomerat e  Combustion  in  the  Laminar  Premixed  Flame 

4.2.1  Empirical  Parameters 

Figure  42  is  a  scanning  electron  microscope  photograph  of  four 
carbon  agglomerated  particles  collected  near  the  injector  exit.  Since 
the  particles  were  irregular,  the  average  of  at  least  eight  particle 
diameters  is  employed  for  the  diameters  quoted  in  the  following.  The 
agglomerate  diameters  of  the  particles  in  Figure  42  vary  within  15%  of 
each  other,  which  was  typical  for  all  of  the  test  conditions. 

Figure  43  shows  two  SEM  photographs  of  particles  collected  at  the 
same  test  condition,  one  just  after  injection  and  the  second  after  it 
had  partially  reacted  by  passing  8. 1  mm  through  the  burner  gases.  The 
nonreacted  particle  has  some  irregularities,  but  is  basically 
spherical  in  shape  and  has  a  relatively  smooth  surface.  After  some 
reaction,  however,  the  particle  becomes  very  porous.  There  appears  to 
be  a  less  dense  pore  layer  surrounding  a  more  dense  central  core. 
This  yielded  an  apparent  average  density  reduction  from  1820  kg/m^  to 
1640  kg/m^.  The  porous  structure  allows  gas  to  percolate  through  the 
particle  and  provides  substantially  more  surface  area  for  reaction 
than  the  surface  area  of  a  smo  th  sphere  having  the  same  outer 
diameter  —  supporting  the  need  for  separate  convection  enhancement 
factors  and  area/reactivity  factors  employed  in  the  present  analysis. 

The  upper  photograph  of  Figure  43  had  a  reacted  mass  fraction  of 
0.0,  and  the  lower  photograph  had  an  0.i22  reacted  mass  raction. 
°articles  burning  in  other  flame  conditions,  which  had  similar  reacted 
mass  fraction  values,  appear  analogous  to  these  photouraohs.  The 


SEM  photograph  or  four  particles  collected  from 
the  5  =  1.0,  =  1683  K,  d^  =  52.6  dm  test 

condition  at  a  distance  of  1.8  rsa  from  the  injector. 
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upper  photograph  of  Figure  44  had  a  reacted  mass  fraction  of  0.902  and 
the  lower  photograph  0.999.  Again  particles  burned  in  other  flame 
conditions  which  had  equivalent  reacted  mass  fraction  values  appeared 
very  much  the  same  as  the  photographs  in  Figure  44.  From  these  four 
photographs,  it  appears  that  as  the  particle  burns,  a  surface  pore 
structure  evolves.  Initially  the  central  core  is  unaffected,  however, 
as  the  pore  structure  develops  even  the  center  of  the  particle 
develops  void  spaces.  Near  the  end.  of  the  combustion  process,  the 
pore  structure  has  grown  to  the  point  where  the  particle  has  an  almost 
lacey  structure  with  a  near  zero  apparent  density,  while  still  having 
a  substaintial  diameter. 

Figure  45  is  an  illustration  of  the  variation  of  particle  density 
with  the  extent  of  reaction,  measured  for  the  various  flame 
conditions.  The  density  variation  is  relatively  independent  of 
initial  particle  diameter  and  flame  conditions  yielding  the 
correlation: 


Qp  «  1820.(1  -  £)0'6  (4.4) 

where  £  is  the  reacted  mass  fraction.  The  greatest  discrepancies 

between  Equation  (4.4)  and  the  measurements  occurred  for  either  the 

very  small  diameter  cases  or  the  first  test  case  investigated.  The 

small  particle  test  conditions  (d  *  10  um)  deviated  from  the  bulk  of 

po 

Che  data  since  experimental  accuracy  has  deteriorated.  The  deviation 

of  the  first  test  condition  examined  (  $  *  0.6,  =  1663  K, 

d  =  73.4  Um)  from  the  bulk  of  the  data  mav  possibly  be  attributed 
po  ' 
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to  experimental  inaccuracies  since  the  data  collection  procedure  had 
not  yet  been  fully  developed,  and  perhaps  the  required  prudence  was 
not  exercised. 

In  addition  to  the  particle  density  variation,  the  pore  structure 
also  gave  rise  to  transport-enhancement  factors  since  the  ambient  gas 
could  percolate  through  the  porous  particle.  The  experimental  values 
of  the  transport-enhancement  factor  for  each  experimental  data  point 
examined  is  shown  in  Figure  46. 

Also  shown  in  the  figure  are  the  transport-enhancement  factors 
for  the  large  particle  test  which  was  taken  to  be  a  constant  as  a 
result  of  insufficient  data.  In  all  the  test  cases,  the  particle 
diameters  ranged  betveea  0  and  75  urn.  As  for  the  large  particle 
tests,  both  the  energy  and  mass  transport-enhancement  factors  were 
found  to  be  equal.  The  empirical  curve-fit  of  the  transport- 
enhancement  factor  which  best  correlated  the  data  is  given  as  a 
function  of  reacted  mass  fraction  and  particle  diameter: 

Y  =  -0.00858dp  +  18. 3e  -13. 2e2  +  2.02  (4.5) 

where  the  particle  diameter,  dp,  is  expressed  in  microns.  Also 
plotted  in  Figure  46  is  the  transport-enhancement  factor  empirical 
curve-fit  for  a  constant  particle  diameter  of  0  and  75  pm.  The 
empirical  curve-fit  for  the  two-constant  particle  diameter  cases  bound 
the  experimental  data. 

Note,  at  the  particle  injector  exit,  where  the  particle  is  nearly 
smooth  with  a  uniform  density;  a  near  unity  enhancement  factor  is 
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indicated.  In  addition,  during  the  first  2051  of  the  combustion 
process  -here  the  particle  is  less  porous,  a  smaller  transport- 
enhancement  factor  is  predicted.  Finally,  the  transport-enhancement 
correlation  based  on  the  small  particle  test  gives  a  reasonable 
prediction  for  the  value  obtained  from  the  large  particle  test  over 
the  last  802  of  the  combustion  process. 

The  increase  in  particle  reactivity  due  to  the  pore  structure  is 
indicated  in  Figure  47  for  the  0^— H2O-H2O  reaction  mechanism.  The 
area-reactivity  multiplication  factors  on  a  zero  suppressed  scale  are 
plotted  in  this  figure  vs.  reacted  mass  fraction  for  the  oxygen, 
carbon  dioxide  and  water  vapor  reaction  mechanisms.  Again  for  each 
experimental  test  condition  the  particle  diameter  ranged  between  0  and 
75  Urn. 

The  empirical  factors  were  established  in  the  following  manner: 
the  values  of  the  transport-enhancement  factor  were  determined  first, 
under  the  constraint  that  the  mass  and  energy  transport-enhancement 
factors  must  be  equal.  The  variation  of  the  transport-enhancement 
factor  with  reacted  mass  fraction  was  determined  by  matching  the 
initial  temperature  increase  period  as  well  as  conditions  where  the 
reaction  rate  was  dif fusionally  controlled.  The  results  were 
correlated  to  give  an  expression  for  the  transport-enhancement  factor 
based  on  particle  diameter  and  reacted  mass  fraction.  Once  the 
transport-enhancement  factor  relationship  was  determined,  an 
optimization  routine  was  employed  uhich  examined  each  individual  data 
point  and  determined  the  appropriate  area— react ivity  multiplication 
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factors  values  to  match  the  race  of  change  of  particle  properties. 
The  results  were  correlated  to  yield  the  relationship  between  the 
area-reactivity  multiplication  factor  values  with  particle  diameter 
and  reacted  mass  fraction. 

The  empirical  curve-fits  of  the  area-reactivity  multiplication 
factors  for  the  0-— GC^-l^O  reaction  mechanismare  given  as: 

a-  =  -0.206d  +  23. 7e  +  68.8 

Cb  p 

aC02  =  -0.290dp  +  43.2£  +  128.1  (4.6) 

a„  .  =  -0.513d  +  50. 1£  +■  178.6 

H20  p 

where  the  particle  diameter,  dp,  is  expressed  in  microns.  The 
empirical  curve-fits  for  the  multiplication  factors  are  plotted  for 
each  reaction  mechanism  for  a  particle  diameter  of  0  and  75  pm.  The 
empirical  curve-fits  bound  the  experimental  data  for  all  three 
reaction  mechanisms. 

Also  shown  in  Figure  47  are  the  large  particle  values  used  for 
the  area-reactivity  multiplication  factors.  The  large  particle  test 

did  not  afford  sufficient  data  to  examine  the  multiplication  factors' 

dependence  on  particle  diameter  or  fraction  of  original  mass  consumed, 
therefore,  constant  valued  multiplication  factors  were  used.  For  each 
reaction  mechanism,  the  area-reactivity  multiplication  factors  for  the 
smaller  particle  test  show  an  overall  increase  of  at  least  20%  above 
the  large  particLe  values. 


AfttA-RtACnvm  MUllIPl IC/dlON  FaCIORS 


igure  AS.  Area-reactivity  multiplication  factors  for  the  OH-O^ 

reaction  ^echaniae.,  fur  all  the  cest  conditions  examined 
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Figure  43  plots  the  area-reactivity  multiplication  factors  vs. 
reacted  nass  fraction  on  a  zero  suppressed  scale  for  the  OH-CH 
reaction  mechanism.  As  before,  the  particle  diameter  in  each 
experimental  test  condition  ranged  between  0  and  75  pin. 

fhe  empirical  curve-fits  of  the  area-rectivity  multiplication 
factor  for  the  hvdroxyl-oxygen  reaction  mechanists  have  a  form  si  7Z1.  X5-T 
to  the  multiplication  factors  for  the  oxygen-carbon  dioxide-water 
vapor  reaction  mechanists: 


a0H  =  -0.234dp  +23. 2t  +  65.3 

an  »  -0.307d_  +  79. 3£  +  902.3 
°2  P 


(4.7) 


where  the  particle  diameter  is  again  expressed  in  microns.  As  before, 
the  multiplication  factor  empirical  curve-fit  for  a  constant  diameter 
of  0  and  73  '^a  bound  the  experimental  data.  Upon  examining  the 
multiplication  factor  values  for  both  sets  of  reaction  mechanisms  for 
each  test  condition  studied,  there  does  not  appear  to  be  a 
relationship  between  the  multiplication  factor  and  either  the  flame 
equivalence  ratio  or  flame  temperture. 

As  earlier.  Figure  48  also  plots  the  large  particle  area- 
reactivity  multiplication  factor  for  both  the  OH  and  0^  reaction 
mechanisms.  Again  the  large  particle  values  were  chosen  to  be  a 
constant  due  to  lack  of  experimental  data.  As  before,  the  small 
particle  multiplication  values  exhibit  an  overall  increase  of  at  least 
20.”  above  the  constant  large  particle  values. 
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4.2.2  Carbon  Agglomerate  Combustion 

Predicted  and  measured  burning  rates  for  the  noncatalyzed 
agglomerate  are  illustrated  in  Figures  49-53.  In  these  figures 
particle  burning  rate  (based  on  rate  of  change  of  the  particle 
diameter  with  time)  is  plotted  as  a  function  of  particle  diameter  with 
either  flame  temperature,  equivalence  ratio  or  initial  particle 
diameter  taken  as  a  parameter.  The  data  for  these  tests  are  tabulated 
in  Appendix  B.  For  these  tests, the  burning  rate  based  on  particle 
diameter  had  an  estimated  error  of  less  than  15%  for  a  75  ym  particle 
and  less  than  30%  for  a  2  in  particle. 

Although  the  maximum  experimental  error  for  the  laminar-premixed 
flame  test  was  less  than  5%,  the  possible  experimental  errors  for  the 
laminar-premixed  flame  test  are  substantially  larger.  Conventional 
error  analysis  gives  the  following  possible  average  experimental  error 
range:  diameter  +10%,  velocity  +  5%,  axial  location  +  20%,  mass  +  28%, 
time  +  20%  and  temperature  +  10%. 

The  theoretical  and  experimental  variation  of  burning  rate  as  a 
function  of  particle  diameter,  with  flame  temperature  as  a  parameter 
is  illustrated  in  Figures  49  and  50.  The  predictions  employ  the 
empirical  correlations  for  the  reaction  parameters  discussed  in  the 
previous  section.  The  agreement  between  theory  and  experiment  is 
excellent,  except  in  the  small  diameter  region  where  experimental 
error  is  large.  Note  that  the  theoretical  and  experimental  results 
for  the  various  diameter  cases  at  a  constant  equivalence  ratio,  as 
shown  in  Figure  50,  line  along  a  smooth  continuous  curve. 
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Figure  50*  Variation  of  burning  rate  vs,  particle  diameter  with 
temperature  at  a  constant  equivalence  ratio  of  1.0. 
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Figure  53-  Variation  of  burning  rate  bs.  particle  diameter  with 

equivalence  ratio  at  a  constant  temperature  of  1950  K. 
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Recalling  that  Che  diffusion-controlled  limit  is  characterized  by 
an  increasing  burning  rate  with  decreasing  particle  diameter,  and  that 
the  kinetically  controlled  limit  is  characterized  by  a  burning  rate 
which  is  relatively  indpendent  of  particle  diameter,  it  is  evident 
from  the  results  in  Figures  49  and  50  that  increasing  the  temperature 
at  a  constant  equivalence  ratio  increases  the  burning  rate,  and  tends 
to  shift  the  reaction  into  the  diffusiunal-controlled  regime. 

The  predictions  employing  either  reaction  mechanism  are  in  exact 
agreement  in  all  diffusion-controlled  regimes,  since,  the  burning  rate 
is  independent  of  reaction  mechanism  in  the  diffusional  controlled 
regime.  The  predictions  of  the  two  reaction  mechanisms  for  the 
kinetically  controlled  cases  agreed  within  10Z,  which  is  comparable  to 
experimental  error.  Therefore,  either  reaction  mechanism  seems  quite 
capable  of  correlating  the  agglomerate  reaction  properties. 

In  order  to  evaluate  the  usefulness  of  the  model  which  employed  a 
more  simple  version  of  the  empirical  parameter  presented  in  this 
investigation,  a  set  of  constant  empirical  parameters  were  chosen. 
The  constant  values  chosen  for  the  02-C07-H20  empirical  parameters  are 
tabulated  in  Table  9.  The  parameters  used  were  obtained  from  Equation 
(4.4)-(4.6)  using  a  particle  diameter  of  35  pm  and  a  reacted  mass 
fraction  of  0.90.  It  was  felt  that  these  sets  of  parameters  were 
representative  and  would  give  a  reasonable  prediction  for  most  test 
conditions.  In  an  attempt  to  further  simplify  the  model,  the  drag 
coefficient  for  a  smooth  spherical  particle  was  used  to  estimate  the 
particle  velocity.  This  simplified  model  was  then  used  to  calculate 
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the  burning  rates,  using  the  02-C02-H2°  reaction  mechanism,  and  the 
results  are  presented  in  Figures  49-53.  For  all  the  conditions 
examined,  the  maximum  error  between  predictions  and  measurements  was 
20".  These  results  indicate  that  the  model  has  some  ability  to  yield 
meaningful  results  even  in  this  simplified  form. 

The  variation  of  agglomerate  burning  rate  with  particle  diameter 
at  constant  temperature,  with  equivalence  ratio  taken  as  a  parameter, 
is  shown  in  Figures  51,  52  and  53.  The  results  indicate  that 

increasing  oxygen  concentrations  (or  lower  equivalence  ratio) 
increases  the  burning  rate,  tending  toward  dif fusionally-controlled 
conditions,  for  these  temperature  levels. 

Predicted  and  measured  reacted  mass  fraction  as  a  function  of 
time  are  illustrated  in  Figures  54  and  55,  with  temperature  as  a 
parameter.  The  agreement  between  predictions  and  measurements  is 
good,  suggesting  adequate  correlation  of  the  empirical  parameters.  In 
general,  increasing  temperatures  cause  increasing  mass  burning  rates, 
with  the  greatest  increase  occurring  for  equivalence  ratios  near 
unity. 

The  variation  of  reacted  mass  fraction  as  a  function  of  time, 
with  particle  diameter  as  a  parameter,  is  illustrated  in  Figure  56  for 
a  fixed  temperature  and  equivalence  ratio.  As  expected,  the  results 
indicate  that  smeller  pcrticles  burn  more  quickly,  with  the  larger 
decrease  in  combustion  time  with  reduced  size  occurring  in 
di f f usionai ly-controlled  regime  —  at  high  temperatures. 
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The  reacted  mass  fraction  as  a  function  of  time,  at  a  constant 
temperature  with  equivalence  ratio  as  a  parameter,  is  plotted  in 
Figures  57-59.  Again  it  is  observed  that  the  agreement  between  theory 
and  experiment  is  good.  For  these  conditions,  decreasing  the 
equivalence  ratio  increases  the  burning  rate;  however,  low  equivalence 
ratios  at  these  high  temperatures  are  difficult  to  achieve  for 
airbreathing  combustion  systems. 

The  final  step  in  evaluating  the  model  for  the  laminar-premixed 
flame  tests  involves  comparison  between  the  predicted  and  measured 
particle-lif e  histories.  These  results  are  presented  in  Figures 
60-73.  Predicted  and  measured  particle  diameter  {normalized  by 
initial  particle  diameter),  particle  mass  (normalized  by  initiel 
particle  mass)  and  particle  temperature  (normalized  by  the  flame 
temperature)  a/e  plotted  as  a  function  of  time  for  both  reaction 
mechanisms.  The  data  used  in  preparing  the  plots  are  tabulated  in 
Appendix  8.  The  agreement  between  theory  and  experiment,  in  general 
is  excellent  for  both  reaction  mechanisms.  An  interesting  feature  of 
these  results  is  that  the  particle  mass  decreases  much  more  quickly 
than  the  particle  diameter.  The  remaining  trends  are  similar  to  the 
behavior  discussed  earlier  for  supported  particles. 

Experimental  uncertainties  were  too  large  to  yield  an  accurate 
determination  of  particle  drag  properties.  Drag  coefficients  were  of 
the  same  order  of  magnitude  as  the  values  for  smooth  spheres  -  but 
were  widely  scattered. 


reacted  mass  fraction 


187 


TIME  (S) 


Figure  58.  Variation  of  reacted  mass  fraction  vs.  time  witii 
equivalence  ratio  at  a  constant  temperature  of 
1815  K  and  a  constant  particle  diameter  cf  73.4  pm. 
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Figure  59.  Variation  of  reacted  mass  fraction  vs.  time  with 
equivalence  ratio  at  a  constant  temperature  of 
1950  K  and  a  constant  particle  diameter  of  73.4  am- 
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Figure  60.  Parcicle- Li f e  hiscory  for  t  =  0.2, 

T,  =  1667  K  and  d  *  73.4  ym. 
f  po 
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Figure  61.  Particle-life  histor; 
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Figure  62.  Particle -l i fe  history  for  $  =  0.6, 
Tp  =  181*  K  and  d  -  73.4  um. 
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CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 


5 . 1  Summary 

The  overall  objective  of  the  present  study  was  to  experimentally 
observe  and  theoretically  predict  the  heat-up  and  combustion  processes 
of  a  slurry  droplet.  The  carbon  slurry  fuels  were  formulated  by 
mixing  noncatalyzed  and/or  catalyzed  medium- thermal  carbon  black 
(roughly  502  by  weight)  with  JP-10  liquid  fuel. 

The  combustion  of  individual-supported  slurry  drops,  400-1000  ym 
initial  diameter,  was  observed  in  an  open  turbulent-diffusion  flame  at 
atmospheric  pressure.  The  environment  of  the  drops  within  the  flame 
was  known  from  measurements  of  mean  velocity,  temperature  and 
composition.  The  drops  were  suspended  from  the  junction  of  a 
thermocouple.  Shadowgraph  motion  pictures  were  also  obtained  for  the 
process.  The  data  provided  the  variation  of  particle  temperature  and 
diameter  as  a  function  of  time. 

A  two-stage  combustion  process  was  observed.  The  first  stage 
involved  heat-up  and  evaporation  of  the  liquid,  leaving  a  solid 
agglomerate  of  the  carbon  particles  on  the  probe-  The  second  stage 
involved  heat-up  and  reaction  or  quenching  of  the  agglomerate. 

A  model  was  constructed  for  both  the  liquid-evaporation  process 
and  the  agglomerate-reaction  process.  This  involved  analysis  of  heat- 
and  mass-transfer  processes  in  the  convective  environment  surrounding 
the  droplet  or  agglomerate,  and  the  chemical  reaction  of  carbon  with 
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oxygen,  carbon  dioxide,  and  water  vapor;  or  oxygen  and  hydroxyl  at  the 
surface  of  the  particle.  Radiative  heat  transfer  between  the 
agglomerate  and  its  surroundings  was  also  considered.  The  predictions 
of  both  of  the  models  were  compared  with  measurements  in  the  flame 
considering  equivalence  ratios  in  the  range  of  0.272-1.350.  The 
liquid-evaporation  model  predictions  for  droplet  si2e  end.  temperature 
agreed  very  well  with  the  experimental  observations,  except  in  regions 
where  an  envelope  flame  was  observed  (equivalence  ratios  of  0.761  and 
0.870),  where  an  extended  model  which  considered  the  presence  of  an 
envelope  flame  was  required  for  the  actual  predictions.  Good 
agreement  between  predicted  and  measured  agglomerate  reaction 
properties  was  also  obtained  throughout  the  test  range  for  both  the 
O2-GO2-H2O  and  the  O2-OH  agglomerate  reaction  model.  Both  models 
accurately  predicted  both  particle  temperature  and  sice  variations. 
Use  of  a  catalyzed  slurry  was  found  to  increase  burning  rates  and 
extended  the  lean  limit  for  agglomerate  reaction. 

A  laminar- premixed  flame  apparatus  was  also  constucted,  in  order 
to  examine  the  variation  of  the  empirical  factors  over  a  broader  range 
of  test  conditions.  These  tests  emphasized  smaller  freely  moving 
particles,  more  representative  of  practical  slurry  sprays,  and  test 
conditions  in  the  kinitically-controlled  regime.  The  apparatus 
consisted  of  a  particle  genorator  which  produced  a  stream  of  carbon 
agglomerate  particles,  with  initial  initial  diameters  in  the  10-100  pm 
range,  which  passed  through  the  post  flame  region  of  a  flat-flame 
burner.  Measurements  were  made  of  particle  size,  mass,  temperature 
and  velocity  as  a  function  of  position  above  the  burner. 
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A  relationship  between  apparent  particle  density  and  reacted-mass 
fraction  was  obtained  which  was  relatively  independent  of  particle 
size  and  local  particle  environment.  The  variation  of  both  the 
transport-enhancement  factors  and  the  area-reactivity  multiplication 
factors,  for  both  reaction  mechanisms,  were  found  to  depend  on 
particle  diameter  and  reacted-mass  fraction  to  some  degree.  However, 
constant  values  similar  to  those  found  for  the  large  supported 
particles,  also  yielded  a  reasonable  correlation  of  the  measurements. 
The  burning  rate  increased  with  increasing  flame  temperature  and  with 
decreasing  equivalence  ratio. 

5.2  Conclusions 

The  major  findings  of  the  study  can  be  summrized  as  follows: 

1.  The  life  history  of  a  carbon  slurry  drop  can  be  broken  down 
into  four  major  stages: 

i.  Heat  up  of  the  particle  to  temperature  levels 
appropriate  for  evaporation  of  the  liquid.  Possible  ignition 
of  liquid  fuel  vapors. 

ii.  Evaporation,  and  possibly  combustion,  of  the  liquid 
fuel  in  a  relatively  conventional  manner  accompanied  by 
agglomeration  of  the  carbon  particles  origianlly  present 
within  the  droplet  into  a  solid  residue. 

iii.  Further  heating  of  the  carbon  residue  once  all  the 
liquid  has  evaporated. 

iv.  Ignition,  or  a  gradually  increasing  reaction  rate  of 
the  carbon,  followed  by  combustion  of  the  residue. 
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2.  The  behavior  of  the  particles  varied  significantly  throughout 
the  flame.  The  various  combustion  regimes  are  identified  in 
Table  6.  They  are  characterized  by  Che  present  or  absence  of 
envelope  flames,  residue  reaction  and  residue  breakup. 

3.  In  general,  the  solid  carbon  in  the  slurry  remained  in  the 

agglomerate,  and  the  heat-up  and  reaction  of  the  agglomerate 
required  the  bulk  of  the  life-time  of  the  particle  (requiring 
90-95 %  of  the  life-time  even  in  locations  where  carbon 
reaction  rates  were  fast).  Agglomerate  reaction  was  only 
observed  for  flame  equivalence  ratios  of  0.37-1.4  for 

uoncatalyzed  slurries  and  0.21-1.4  for  catalyzed  slurries. 

4.  The  carbon  residue  broke  up  into  smaller  fragments  in  the 
high  relative  velocity  region  near  the  injector.  Otherwise, 
most  of  the  carbon  formed  a  single  residue  particle  whose 
size  was  controlled  by  the  initial  drop  size.  In  this 
circumstance,  the  desirability  of  fine  atomization  for 
practical  combustors  is  obvious. 

5.  Addition  of  the  catalyst  significantly  increased  the  size  of 
the  region  of  the  flame  where  reaction  of  the  residue 
occjrred.  Catalyzed-slurry  residue  also  reacted  more 
completely. 

6.  Following  active  carbon  combustion,  the  residue  had  an  open 
lacey  structure  indicating  that  reaction  was  occurring  to 
some  depth  in  the  particle.  The  sub-elements  of  this 
structure  were  approximately  the  same  size  as  the  ultimate 
particle  size  in  the  carbon  black  (0.3  pm). 
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7.  Cenospheres  were  observed  as  carbon  residue  in  some  cases. 
It  is  postulated  that  they  result  from  quenching  of  the 
reaction  at  the  outer  surface  of  the  particle,  due  to 
excessive  radiative  and  convective  heat  losses  to  the 
environment . 

8.  Similar  to  earlier  work  [11,150,151,155,156],  good 
predictions  of  liquid-gasification  rates  and  life  histories 
were  obtained  by  selecting  a  fixed  reference  condition  for 
evaluating  average  properties.  Drop  gasification  could  be 
treated  as  a  simple  evaporation  process,  except  for  a  narrow 
region  where  envelope  flames  were  observed.  The  presence  of 
an  envelope  flame  resulted  in  an  underestimation  of 
gasification  rates,  when  the  drop  was  assumed  to  evaporate; 
however,  a  conventional  extention  of  the  model  to  consider 
envelope  flames  rectified  this  difficultly. 

9.  Good  predictions  of  agglomerate  reaction  rates  and  life 

histories  were  obtained  using  either  the  Neoh  et  al.  [121], 
or  extended  Libby  and  Blake  [88,89]  reaction  approximations 
over  a  wide  range  of  test  conditions,  including  particle 
diameters  between  0  and  1000  jjra,  flame  equivalence  ratios 
between  0.2  and  1.4,  and  flame  temperatures  between  500  and 
2000  K.  This  agreement  was  obtained  by  selecting 

empirical  area-reactivity  factors  which  allow  for  effects  of 
pores,  intrinsic  surface  reactivity  and  catalyst,  for  each 


reaction;  and  transport-enhancement  factors  which  accounted 


for  the  ambient  gas  porcolating  throught  the  porous  particle. 
These  factors  were  significantly  greater  than  unity, 
suggesting  that  pore  diffusion  and  reaction  within  pores  was 
important  for  present  test  conditions.  The  area-reactivity 
multiplication  and  the  transport  enhancement  factors  were 
found  to  depend  on  both  particle  diameter  and  reacted-mass 
fraction  to  some  degree.  The  values  for  the  transport- 
enhancement  factors  for  the  small  particle  test  had  values 
reasonably  close  to  the  large  particle  values.  The  vallues 
for  the  area-reactivity  multiplication  factors  for  the  small 
particle  tests  exhibit  approximately  a  20Z  overall  increase 
above  the  values  used  for  the  large  particle  tests. 

10.  Adding  catalyst  increased  agglomerate  burning  rates,  except 
at  the  diffusion-controlled  limit,  and  also  extended  the  lean 
limit  for  quenching.  Agglomerate  shattering  in  the  lean 
region  of  the  flame  was  also  occasionally  observed  for 
catalyzed  slurries. 

5-3  Recommendations  for  Further  Study 

The  present  model  employs  an  oversimplified  treatment  of 
agglomerate  structure;  involving  the  assumption  of  a  smooth  surface 
with  empirical  convection  enhancement  and  area-reactivity  factors  to 
aliow  for  effects  of  flow  penetrating  the  porous  surface,  pore 
diffusion,  intrinsic  surface  reactivity  and  catalyst.  These  empirical 
factors  were  found  to  be  nearly  constant  over  the  present  test  range. 
Since  the  effects  of  surface  structure  probably  vary  with  particle 
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size,  percent  reactivity,  properties  of  the  carbon  black,  and  flow 
conditions,  it  is  unlikely  that  the  empirical  factors  found  in  the 
present  study  are  generally  appropriate. 

In  order  to  further  expand  the  level  of  understanding  of  the 
carbon  agglomerate  combustion  process  the  following  items  are 
recommended  as  possible  areas  of  further  study: 

1.  Evaluate  the  effect  of  ultimate  carbon  black  particle  size  on 
agglomerate  reaction.  In  particular,  a  bi-modal  blend  of 
ultimate  particle  particle  size  in  the  slurry  should  be 
considered  since  such  formulations  are  likel>  to  provide 
improved  flow  and  stability  properties. 

2.  Examine  the  effect  of  ambient  pressure  on  the  agglomerate 
combustion  process,  particularly  emphasizing  the  higher 
pressure  levels  typically  found  in  combustors. 

.  Develope  a  pore  structure  model  to  aid  in  further 
understanding  of  the  microscopic  combustion  processes. 
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APPENDIX  A 


COMPUTATION  OF  PROPERTIES 
A. 1  Mean  Property  State 

The  average  properties  of  the  gas  surrounding  the  droplet  were 
computed  at  the  mean  state,  defined  as  follows: 

^avg  °  a9p  (1  "  a)  (A.l) 

where  d  is  a  generic  property  at  a  specific  temperature  and  mass 
fraction,  a  can  range  between  0  and  1  and  a  value  of  0.9  was  selected 
for  the  present  investigation  as  suggested  by  Faeth  and  coworkers 

[15",! 54, 155]. 


A. 2  Mixture  Thermal  Conductivity 

The  thermal  conductivity  of  the  gaseous  mixture  was  calculated 
with  the  Mason  and  Saxena  formulation  of  the  Wassiljewa  equation 
[160] . 


where 
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j=l 


(A.2) 
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ij 


1  ,  iir-i/v  ii  +  /M/2  /A ,2 

P72  <1  +  mJ>  ll~  11  +  (TT)  1  . 


(A.  3) 


The  thermal  conductivity  and  viscosity  for  gases  other  than  JP-10 
was  obtained  from  Reference  161.  The  properties  of  JP-10  are  derived 


in  Section  A. 6. 


Equation  (A.2)  was  evaluated  at  the  average 
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temperature  and  composition  of  the  gas  phase  in  Equation  (A.l).  In 
Equation  (A. 3),  -j  is  set  equal  to  1,  when  i  «  j. 


A. 3  Mixture  Viscosity 

The  viscosity  calculations  of  mixture  employed  the  method  of 
Wilke,  cited  in  Reference  160. 

For  a  mixture  of  N  components 
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i=l 
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(A. 4) 


where  j  is  given  by  Equation  (A. 3).  The  viscosity  data  used  in  the 
calculations  was  obtained  in  the  same  manner  as  described  in  Section 
A. 2. 


A. 4  Mixture  Dif fusivity 

Since  only  liquid  vapor  diffuses  through  the  stagnant  mixture, 
when  liquid  JF-10  is  present  the  diffusion  coefficient  is  given  as 


lm 


X  - 


j=2  Dlj 


-1 


D  . 


(A. 5) 


Subscript  i  represents  fuel  vapor,  and  j  indicates  all  other  gases. 
Binary  difrusivity  calculations  employed  Chapman-Enskog  theory  along 
with  Lennard-Jones  characteristic  energy  and  length  concept  (<■;'  and 
a  ),  cited  in  Reference  161.  The  data  for  C?  and  e'/k.'  are  listed  in 


Table  10. 


Table  10. 


Properties  Used  for  Binary  Diffusivity  Calculations 


a 


Substance 

a  (A°) 

e'/k?  (K) 

h2° 

2.641 

809.1 

CO 

3.690 

91.7 

co2 

3.941 

195.2 

N2 

3.798 

71.4 

H2 

2.827 

59.7 

°2 

3.467 

106.7 

JP-10b 

6.927 

531.5 

aTaken  from  Reid  et  al .  [161 J . 
^Determined  in  Section  A. 6. 
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When  all  of  the  liquid  fuel  has  evaporated  and  only  a  carbon- 
agglomerated  particle  regains  the  diffusion  coefficient  of  the  gas  is 
given  by: 


D 


1 

2 


N  X 
2  — 
L  j=2  D 


(A*6) 


A.  3  Mixture  Specific  Heat 

The  mixture  specific  heat  was  calculated  as  a  mole  fraction 
weighted  average  of  the  specific  heats  of  each  contributing  species 
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i=l 


Y  C 

Xi  Pi 
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(A.  7) 


A. 6  Yheraophysical  and  Transport  Properties  for  JP-10 

JP-10  or  exo-tetrahydrodi  (cyclopentadiene)  has  the  chemical 

formula  and  a  molecular  weight  of  136.2,  the  chemical  structure 

is  shown  in  Figure  74. 

The  critical  properties  of  JP-10  were  calculated  from  the 
chemical  structure  using  Lyderson's  method  discussed  on  page  12  of 
Reid,  Prausnitz  and  Sherwood  [161)  to  yield' 

Critical  Temperature,  Tc  =  672.4  X 

Critical  Pressure,  Pc  =  27.84  atm 

Critical  Volume,  VQ  =*  3.605  x  10~^  n^/kg. 


(A.d) 
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The  correlation  for  vapor  pressure  liquid  density  and  liquid 
specific  heat  was  obtained  by  curve-fitting  experimental  data  obtained 
from  McCoy  [163J.  The  vapor  pressure  equation  determined  where  the 
pressure  is  in  N/m  was 

P  -  3.0690xl09  exp(~— (a.9) 

i  • 

Liquid  density  in  kg/m3  was  correlated  as 

P  =  1166.4  -  0.792(T).  (A. 10) 

And,  the  liquid  specific  heat  in  units  of  cal/gm°C  was  determined  as 

Cp  =  0.0615  +  1.08  x  10"3(T)  (A. 11) 

where  the  temperature  in  Equation  (A-9)-(A.ll)  is  in  degrees  Kelvin. 

The  heat  of  vaporization  was  obtained  using  the  Clausius- 
Clapayron  equation  assuming  ideal  gas  behavior,  and  that  the  vapor 
specific  volume  is  much  greater  than  the  liquid  specific  volume,  which 
leads  to 

ifg  -  286.7  kJ/kg.  (A. 12) 

The  gas-phase  specific  heat  was  calculated  using  the  method  of 
Rehani  and  Doraiswany  discussed  on  page  234  of  Reid,  Prausnitz  and 
Sherwood  [I61J  to  give: 
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C  =  0.35669  +  0.23542xlQ-2  (T)  -  1.5827xlO-6  (T)2 

-10  3  (A*13) 

+  3.9699x10  10  (T)3 

where  temperature  Is  again  in  degrees  Kelvin. 

The  gas-phase  viscosity  for  JP-10  was  calculated  using  the 
Reichenberg  method  discussed  in  page  404  of  Reid,  Prausnitz  and 
Sherwood  [161).  The  gas-phase  viscosity  equation  determined  was 

U  *  170.32(Tr)/[l  +  0.36<Tr)(Tr-l))1/6  (A. 14) 

where 

Tr  -  T/Tcr  (A. 15) 

and  T,  Tcr  are  in  degrees  Kelvin  and  has  the  units  of  micro-poise. 

An  equation  for  the  gas-phase  thermal  conductivity  was  obtained 
by  the  Miscie  and  Thodos  method  discussed  on  page  480  of  Reid, 
Prausnitz  and  Sherwood  [161).  The  resulting  equation  with  thermal 
conductivity  in  units  of  cal/cm  s  K  was 

X  -  26.59  x  i0-8  (14.52(Tr)  -  5.14]2/3  Cp  (A. 16) 

where  Cp  and  Tr  are  given  in  Equations  (A. 13)  and  (A- 15)  respectively. 

The  method  discussed  by  Shehla  [ 162 J  was  used  to  obtain  the 
Lennard-Jones  characteristic  energy  and  length  employed  in  the 
Chapmau-Enskog  binary  diffusivity  calculations.  The  resulting 
Lennard-Jones  parameters  were: 

e'/k'  »  504.6  K 


o 


6.927  A°. 


(A. 17) 


o  o  o 
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A. 7  Equations  Used  to  Calculate  Transport  and  Thermophysical 
Properties  of  Each  Species 


The  equations 

developed 

to 

determine 

the  transport 

and 

thermophysical  properties  of  the 

gas- 

-phase  for 

the  species  02, 

co2, 

co,  h2o,  h2,  OH, 

and  N2  were 

obtained 

by  curve-fitting 

hign 

temperature  data  tabulated  by  Shehla  [161].  The  equations  used  to 
determine  the  gas-phase  properties  for  JP-10  were  discussed  in  the 
previous  section.  The  computer  subroutine  to  calculate  the  gas-phase 
properties  that  was  used  in  the  carbon— agglomerate  analysis  is 
presented  in  the  remainder  of  this  section.  The  nomenclature  used 
within  the  program  is  explained  in  the  program  comments. 


SUBROUTINE  PROPG  (YI,YS,TI,TS,P,VM,CDM,DFM,RHM,CPM,AVMS, 

$X02 , XCO , XC02 , XH2 , XH20 , XOH) 

C - SPECIFY  PROPERTIES  OF  AMBIENT  GAS - - - 

C - SPECIES  WHICH  EXIST  IN  FLOW  FIELD  : 

C  (1) - HYDROXYL  (OH) 

C  (2) - WATER  VAPOR  (H20) 

C  (3) - CARBON  MONOXIDE  (CO) 

C  (4) - CARBON  DIOXIDE  (C02) 

C  (5) - NITROGEN  (N2) 

C  (6) - HYDROGEN  (H2) 

C  (7) - OXYGEN  (02) 

- - 

IMPLICIT  REAL*8(A-H,0-Z) 

DIMENSION  AM(7),XS(7),XI(7),X(7),Y(7)tYS(7),YI(7) 

DIMENSION  DF(7 ,7) ,  D(7),  CP(7),  V(7),  CD(7) 

DIMENSION  EDKC(7,7),  SIGMA(7,7),  0M(7,7),  PHI(7),  EDK(7) ,  SIGM(7) 
DO  20  1=1,7 

TO  CHECK  MASS  FRACTION  VALUES  REMOVE  C'S 

IF  ((YS(I).GT. 1 .0D0) .0R.(YS(I) .LT.O.ODO))  WRITE(6,lO)  I,YS(I) 

10  FORMAT ( '  ',5X,'YS(',I1,')  =',1PD13.6) 

IF  (YI(I). LT.O.ODO)  YI(I)=O.ODO 
IF  (YS(I). LT.O.ODO)  YS(I)-O.ODO 
IF  (YI(I).GT.l.ODO)  YI(I)-1.0D0 
IF  (YS(I) .GT. 1.0D0)  YS(l)=l,0DQ 
20  CONTINUE 

C - LIST  MOLECULAR  WEIGHT  OF  EACH  SPECIES 

AM( I )=l 7 .0074DQ 
AM(2)=13.016D0 
AM( 3 )=23 .OllDO 
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AM(4)=44.01lD0 

AM(5)-28.016D0 

AM(6)=2.016D0 

AM(7)=32.00000 

C - COMPUTE  MASS  FRACTION  OF  FUEL  AT  SURFACE 

AVMI=O.ODOO 
DO  30  1*1,7 
AVMI=AV.MI+YI (  I )/AM( I ) 

30  CONTINUE 

AVMI=1 .ODOO/AVMI 
PRESS=P*101325.D0 
DO  40  1*1,7 

XI < I )= AVMI*YI( I )/ AM( I ) 

40  CONTINUE 
AVMI=O.ODOO 
DO  30  1=1,7 
AVMI=AVMI+AM( I ) *XI C I ) 

50  CONTINUE 

RHI=PRESS*AVMI/TI/3314. 3D00 
AVMSO.ODOO 
DO  60  1=1,7 
AVMS=AVMS+YS(I)/AM(I) 

60  CONTINUE 

A  VMS'*  1 .  DO/AVMS 
DO  70  1=1,7 

XS(I)=AVMS*YS(I)/AM(I) 

70  CONTINUE 
X02=XS(7) 

XC0=XS(3) 

XC02=XS(4) 

XH2=XS(6) 

XH20=XS<  2) 

XOH-XSCi} 

C - A  IS  A  MULTIPLICATION  FACTOR  APPEARS  IN  TB=A*TS+<  1-A) *TI  EQUATION 

C - IT  WAS  OBTAINED  BY  MATCHING  WITH  CALIBRATION  TEST. - 

A=0.9D00 

TB=A*TS+( 1 ,ODOO-A)*TI 
DO  80  1=1,7 

X(I)=A*XS(I)K1.0D00-A)*XI(I) 

80  CONTINUE 
AVM=O.QDGO 
DO  90  1=1,7 
AVM=AVM+AM(  I  )*X(  I ) 

90  CONTINUE 
DO  10Q  1=1,7 
Y(I)=AM(I)*X(I)/AVM 
100  CONTINUE 

C - SPECIFIC  HEAT  ANU  TRANSPORT  PROPERTIES  (UNIT) - 

C  C?  SPECIFIC  HEAT  AT  CONST.  PRESSURE  (CAL/G-M0LE*DEG) 

C  — * - (KJ/KG*DEG)  1/HEN  RETURNING - 

C  V  VISCOSITY  (MICRO  POISE) 

C  - (KC/M*S  OR  PA*S  )  WHEN  RETURNING- - 
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C  CD  THERMAL  CONDUCTIVITY 

C  - (KJ/M*S*DEG) 

C  Dr  MASS  DIFFUSIVITY 

C  - (M**2/S) 

C  RH  DENSITY 

C  - (KG/M**3  ) 

C 


C - MEANING  OF  SUFFIX- 


( *1 .0D-06  CAL/CM*S*DEG) 

WHEN  RETURNING - 

(CM**2/S) 

WHEN  RETURNING - 

(KG/H**3) 

WHEN  RETURNING - 


C  I - INFINITE  M - MEAN  VALUE  S - SURFACE 

C  F - FUEL  GAS  A - AMBIENT  GAS  B - BULK  GAS 

IF  (TB.GE.500.DO)  GO  TO  110 

CP( 1 )=7 . 67484 2DO-2 . 38931D-3*TB+2 . 2992D-6*TB**2 
CP( 2 )=7 . 972473D0-8 .71300D-4*TB+3 . 5127D-6*TB**2 
CP(  3 >=*7 .  18284  1D0-L  -  6 2840D-3*TB+3 . 0095D-6*TB**2 
CP(  4  )=4 . 7987  56D0+1 .  65402D-2*TB-9 . 61 1  5D-6*TB**2 
CP(  5 )**7 . 174363D0-1 .46168D-3*TB+2 . 50l9D-6*TB**2 
CP(o)=6.274392D0+3.00860D-3*TB-3. 1429D-6*TB**2 
CP(  7  )=*6 .87508500-4 . 35900D-4*TB+3 . 0954D-6*TB**2 
GO  TO  120 

no  tt=tb/iooo.do 

CP(1)=5.24639000DO+2.42980000DO*TT-5.66500000D-1*TT**2 
$  +4 . 90800000D-2*TT**3+1 . 87  38000 D- 1 /TT**2 

CP<  2 )=6 . 4 107445 1 DO+4 . 40994628D0*TT-8 . 875  2 1069D-1 *TT**2 
$  +6 . 28994081 D-2*TT**3-9 . 2660052D-2/TT**2 

CP( 3)=0. 72635540D1+0. 1 t035250Dl*TT-0 . 2358 I1400D0*TT**2 
$  +0. 17634440D-l*TT**3-O.15939236DO/TT**2 

CP(4 )=0. 12226062D2+0. 178935 10Dl*TT-0. 380 1322Q0D0*TT**2 
$  -K) . 29 1 7 1 2 1 OD- 1 *TT**3-0 . 59083890D0/TT**2 

CP(5)=0.69731460Di+0.12742370Di*TT-0.274262700DO*TT**2 
$  +0 . 2Q276900D- 1 *TT**3-0 .1185  5850D0/TT**2 

CP( 6)-Q. 52903300D1+0. 2024 1900Dl*TT-0 . 352047000D0*TT**2 
$  +0.24859530D-l*TT**3+0. 19386990D0/TT**2 

CP( 7 )*0. 8021 1900D 1+0.  54604200DO*TT-K) .  77746000D-2*TT**2 
$  -0.61 l56400D-2*TT**3-0 .2l60974ODO/TT**2 

120  IF(TB.GT.i. D3.AND.TB.LT. 2. 5D3)  GO  TO  i30 
IF(TB.GE.2.5D3)  GO  TO  140 
GO  TO  150 
130  TT-TB-IOOO.DO 

V ( 1 )- 44  3 . 2449D0+2 . 88  3 1 5  D- 1 *TT-4 . 1 5  84D-5  *TT*  *2+6 . 743  4D- 9*TT**3 
V(2)=358. 5945D0+3 . 5 1 007D-1 *TT-5 . 3575D~5*TT**2+7 . 7Q15D-9*TT**3 
V(3)=405.715  3D0+2 . 5  89  7 1 D- 1 *TT-3 . 1 3  70D-5  *TT*  *2+  3 . 8 1 3  8D-9  *TT**3 
V(4 )=393 . 1504D0+2 . 69698D- l*TT-4 . 5523D-5*TT**2+7 .71 99D-9*TT**3 
V(5)=397. 187 7DO+2. 59811 D-l*TT-4 .0295D-5*TT**2+7 . 3 i 83D-9*TT**3 
V(6)  =  197. 6894DO+1 . 27990D-  1*TT-1 .9142D-5*TT**2i-3.4007D-9*TT**3 
V( 7 )=480.9907DO+3 .02530D-1 *TT-3 . 3070D-5*TT**2+3. 29  33D-9*TT**3 
GO  TO  160 
140  TT»TB-2500.D0 

V(  t)=>304. 3804 DO+2. 079240- l*TT-l  .4765D-5*TT**2+2.033iD-9*TT**3 
V(2)=790. 5G20DO+2.42314D- l*TT-2 . 4886D-5*TT**2+5 . IQ63D-9*TT**3 
V( 3)*7  36 .  10 1  lDO-r  1 . 900 16D- i *TT- 1 . 2502D-5  *TT**2+ 1 .  16  32D-9*TT**3 
V(4  )=*72  1 . 363900+1 .83 1  92D- 1  *TT+0 .0200D-5*TT**2-5 . 7908D-9*TT**3 
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V(5)=720.8079D0+1.85955D-1*TT-1 .2555D-5*TT**2+i.3l52D-9 *TT**3 

V(6)=358. 195600+9. 17810D-2*TT-Q.4861D-5*TT**2-2. 1608D-10*TT**3 
V(7)=371 . 5132D0+2. 26 188D-1*TT-1 . 6785D-5  *TT**2+2.4879D-9*TT**3 
GO  TO  160 

150  V( 1 )=34.9350D0+6.2460D-l*TB-3. 1150D-4*TB**2+9 . 5453D-8*TB**3 
'  V( 2 )-  21 . 2350D0+2 . 5099D-1  *TB+1 . 6890D-4*TB* *2-8 . 257  1D-8*TS**3 

V  (  3  )= 2  7 . 0  290D0+5 . 82  7 OD- 1 *T3-2 . 9  5  70D-4  *TB*  *2+9 . 19  2 1 D-8  *T3*  *3 

V  (  4  )  =-S .  5  780D0+6 . 2 1 44D- 1 *TB- 3 . 1 44  OD-4  *TB*  *2+9 . 4  699  D-3+TB*  *3 

V  ( 5 ) = 32 . 7  427  DO+5 . 6251 D-l*TB-2. 9181 D-4  *TB* *2+9 . 3  694D- 8*T  B*  *3 
V ( 6 ) =1 6 . 39  53DO+2 .867  5D-1 *TB-1  . 6439D-4*TB**2+5 . 9090D-8*TB**3 
V ( 7 ) = 2  2 . 5 1 8000+7 .  1 7  90D- 1 *TB-3 . 8  5 1 OD-4  *TB*  *2+ 1 . 2  550D- 7*T3* *3 

160  CONTINUE 
C 

IF(TB.G+.l. D3.AND.TB.LT. 2. 5D3)  GO  TO  170 
IF(TB.GE.2.5D3)  GO  TO  180 
GO  TO  190 
170  n>TB-lQOO.DO 

CD( l)-274. 4993DO+2.21650D-1*TT-0.8940D-5*TT*  *2-2. 7137D-9*TT**3 
CD( 2 )=2  76 . 3345DO+3 . 60970D- 1 *TT+0 . 1080D-5*TT**2-9 . 896 1 D-9*TT**3 
CD( 3)-164. 5385DO+1 .3 1915D- 1*TT-1 .8360D-5*TT**2+2. 2402D-9*TT**3 

CD(4)-I60.5447D0+1.44716D-1*TT-2.7369D-5*TT**2+4.6602D-9*TT**3 

CD( 5)-158.8701D0+l . 30969D-l*TT-2 .0384D-5*TT**2+3. 1328D-9*TT**3 
CD( 6 )= 1021 .868D0+7 .86030D-1 *TT+0 - 6080D-5*TT**2- 1 . 7096D-8*TT**3 
CD( 7)-l 78.8670DO+1 .32902D-1*TT-1 . 2880D-5*TT**2+2. 6768D-9*TT**3 
GO  TO  200 
130  TT-TB-2  500 . DO 

CD(  1  )-57 7.9192DO+1 .79665D-1*TT-1 . 3345D-5*TT**2+1 . 88l9D-9*TT**3 
CD ( 2 )-737 . 2 1 1 9DO+3 . 0429  7D-1 *TT-2 .767  4D-5  *TT**2+4 . 3856D-9*TT**3 

CD(3)»  328. 7.07  5D0+9.24I40D-2*TT-0.7259D-5*TT*  *2+1. 4081D-9*TT**3 
CD ( 4 ) -33 1 . 7 3 1 2 DO+9 . 3 1 960D-2 *TT-0 . 1 3 1 l D-5 *TT* *2- 2 . 4 11 5D- 9*TT* *3 
CD( 5)-319 .941 400+9 . 19400D-2*TT-1 .0570D-5*TT**2+2. 8529D-9*TT**3 
CD( 6)-2 1 57 . 398DO+6 . 98236D-1 *TT-4 . 4239D-5  *TT**2+5. 26 1 7D-9*TT**3 
CD(7)-358. 1880D0+1 . 1 1655D-1*TT-Q .4735D-5*TT**2-1 .91 53D-lt*TT**3 
GO  TO  200 

190  C0( 1)*20. 6530D0+3.7420D-1 *TB-2 .054 OD-4 *TB* *2+8. 560 3D-8*TB**3 

CD(2)-19.4052DO+1. 1330 D-1*TB+1.9945D-4*TB* *2-5. 5305D-8*TB**3 
CD(3)=17. 1240D0+1 .5949D-l*TB-0. 1020D-4*TB**2- l .76840-9 *TB**3 
CD(4)— 14. 414D0+1.9866D-1*TB-1. 71 90D-5*TB* *2-6. 56680-9 *TB**3 
CD(5)-18.2390D0+1.6053D-l*TB-2.99  30D-5  *TB*  *2+ 1 . 0 10 1 D-3  *TB*  *3 
CD( 6)=99.760OD0+l .3812 D+0*TB- 7 . 9  830D-4  *TB*  *2+3 . 398  OD-7  *TB*  *3 
CD (7)- 10. 5004  DO+1 . 86440- 1 *TB-0 . 5  740D-5*TB*  *2- 1 . 2626D~8*TB*  *3 
200  CONTINUE 
C 

SIGM( 5)=3 . 798D0 

SIGM(7)=3.467D0 

EDK(5)-71.4DO 

EDK(7)»106.7D0 

1-5 

J-7 

EDKC( I , J)-DSQRT(£OK( I)*EDK(J )) 

SIGMA ( I , J )=(S ICM( I )+SIGM( J ) )/2 . DO 
VBCO.47535D0*(TB/5DKC(  1 FJ)) 


o  r> 
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VBB-l .52996D0*(TB/EDKC(I ,J)) 

VBA-3 .8941  IDO* (TB/EDKC( I , J  )  ) 

IF  (VBC.GT.174.6D0)  OM(I , J)-l .06036D0/( TB/EDKC( I, J))**0. 1561 
IF  ( VBC.GT. 174 . 6D0)  GO  TO  210 

IF  (VBB.GT.174.6DG)  OM( I , J)»l .06036D0/(TB/EDKC( I, J) )**0. 1561+ 

$  O.l93DO/DEXP(0.47635DQ*(TB/EDKC(I»J))) 

IF  (VBB.GT.174.6D0)  GO  TO  210 

IF  (VBA.GT.174.6D0)  0M(I,J)»1 .06036DO/(TB/EDKC(  I, J))**0.  1561  + 

$  0.193DO/DEXP(0.47635DO*(TB/EDKC(I,J)))+1.3O537DO/ 

$  DEXP( 1 . 52996DQ*(TB/EDKC(I, J ) ) ) 

IF  (VBA.GT.174.6D0)  GO  TO  210 

0M( 1 , J)* 1 .06036D0/ (TB/EDKC( I ,  J) )**0. 1561+0. 19 3D0/DEXP(Q. 47635DO 
§  *(TB/EDKC( I, J)))+l .03587D0/DEXP( 1 . 52996D0* (TB/SDKC( I,  J  ) )) 

$  +1 .76474D0/D£XP(3.8941 1DQ*(TB/EDKC( I, J) )) 

210  DF(I,J)-1.85SD-3*TB**1.5*DSQRT((AM(I)+AM(J))/(AM(I)*AM(J))) 

$  /P/SIGMA( I , J)**2/0M( I , J) 

- - - CALCULATION  OF  MEAN  VALUES - 

- VISCOSITY - - - 

VM=0 . DO 
DO  230  1-1,7 
XPHI-Q.DO 
DO  220  J-1,7 

PHI( J)«(l.DO+(V(I)/V(J))**.5*(AM(J)/AM(I))**.25)**2/ 

$  (8.D0*(1.D0+(AM(I)/AM(J))))**.5 

IF(J.EQ.I)  PHI(J)«l.DO 
XPHI-XPHI+X(J)*PHI(J) 

220  CONTINUE 

VK- VM+-X  ( I )  *V<  I )  /  XPHI 
230  CONTINUE 

VM=VM* 1 .0D-07 

C - CONDUCTIVITY - 

CDM-O.DO 
DO  250  K=i ,7 
XP11I-0.D0 
DO  240  L-l,7 

PHI(L)-(l.DO+(V(K)/V(L))**.5*(AM(L)/AM(K))**.25)**2/ 

$  (8.DO*(l.DO+(AM(X)/AM(L))))**.5 

IF(L.EQ.K)  PUI( L)- 1 . DO 
XrtiI-X£-nI+Xa)*PHI(L) 

240  CONTINUE 

CDM=CDM+X( X)*CD( K)/XPHI 
250  CONTINUE 

CDM-CDM* 1 . OD-Q  7 
CDM=4. 1879D0*CDM 

C - DI7FUSIVITY - 

DFM=DF( 5 , 7)*1 .D-4 

C - DENSITY - 

RHM=PR£SS*AVM/TS/33  14. 3D00 

C - SPECIFIC  HEAT - - - 

CPM-O.ODOO 
DO  260  1-1,7 
CPM-CPM+X( I)*CP( I ) 
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260  CONTINUE 
CPM=CPM/AVM 
CPM-4. 1868D0*CPM 
IF  (VM.LT.O.DO)  WRITE(6,270)  VM 
IF  (CDM.LT.O.DO)  WRITEC6.280)  CDM 
IF  (DFM.LT.O.DO)  WRITE(6,290)  DFM 
IF  (RHM. LT.O.DO)  WRITE(6,300)  RHM 
IF  (CPM. LT.O.DO)  WRITE( 6,310)  CPM 
IF  (AVMS. LT.O.DO)  WRITE(6,320)  AVMS 
270  FORMAT( '  ' , 10X , ' VM  HAS  THE  VALUE' , 1PD15.6) 
280  FORMATC '  '.lOX.'CDM  HAS  THE  VALUE' , 1PD15.6) 
290  FORMATC'  ',10X,'DFM  HAS  THE  VALUE' , IPDI5.6) 
300  FORMATC'  ' , 10X, 'RHM  HAS  THE  VALUE' , 1PD15.6) 
310  FORMATC'  '.iOX.'CPM  HAS  THE  VALUE' , 1PD15.6) 
320  FORMATC'  '.lOX.'AVMS  HAS  THE  VALUE' , 1PDI5.6) 
RETURN 
END 
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EXPERIMENTAL  DATA 

B.l  Data  From  the  Turbulent  Diffusion  Flame  Test 


3.1.1  Droo  Environment  Measurements 


Temperature  Measurements 


x/d  Temperature  (K) 
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42.5 
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74.5 

1008 

85.0 
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1543 

170.0 
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1541 
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1271 
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Time 

(s) 


0.00 

0.01 

0.05 

0.08 

0.10 

0.13 

0.16 

0.19 

0.22 

0.26 

0.29 

0.45 

0.52 

0.63 

0.78 

0.84 

1.14 

1.37 

1.56 

1.85 
2.07 
2.33 
2.52 
2.69 

2.86 
3.11 
3.24 
3.40 
3.61 
3.88 
4.26 
4.46 
4.73 
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,  x/d  =  170.0,  $  =  1.350 


d 

T 

P 

P 

{ym) 

(K) 

1230 

294 

1229 

311 

1224 

321 

1216 

369 

1190 

384 

1162 

402 

1112 

410 

1070 

413 

1056 

414 

1020 

415 

1002 

416 

996 

417 

987 

449 

985 

506 

981 

628 

973 

974 

951 

1329 

926 

1369 

914 

1458 

883 

1484 

855 

1506 

817 

1502 

783 

1543 

761 

1587 

746 

1558 

707 

1584 

667 

1587 

641 

1595 

625 

1607 

589 

1618 

524 

1632 

481 

1661 

432 

1705 

417 

1743 

Time 

(s) 


0.00 

0.03 

0.07 

0.10 

0.14 

0.18 

0.22 

0.26 

0.29 

0.41 

0.52 

0.59 

0.74 

0.81 

0.85 

0.91 

1.13 

1.22 

1.45 
1.66 
1.89 
2.05 
2.30 
2.50 
2.75 

2.97 
3.18 

3.45 
3.72 

3.98 
4.28 
4.64 
5.06 
5.37 
5.69 
5.92 

6.14 
6.27 
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-Life  History,  x/d  =  212.5,  $  -  0.939 


d 

p 

(ym) 

T 

P 

<K) 

1261 

299 

1258 

344 

1248 

377 

1226 

410 

1189 

423 

1154 

431 

1113 

433 

1081 

434 

1054 

436 

1023 

480 

1008 

554 

1006 

685 

1005 

656 

1004 

784 

1002 

1030 

993 

1169 

988 

1272 

960 

1334 

964 

1395 

927 

1456 

940 

1444 

909 

1476 

899 

1484 

864 

1497 

861 

1479 

822 

1481 

788 

1526 

779 

1534 

742 

1537 

‘05 

1517 

698 

1559 

655 

1539 

589 

1541 

582 

1575 

534 

1601 

502 

1605 

458 

1682 

447 

1724 

Particle-Life  His: 


Tisae  After 
All  Liquid 
has  Evaporated 
(s) 

0.00 

0.15 

0.30 

0.50 

0.90 

1.21 

1.62 

2.01 

2.41 
2.83 
3.24 
4.01 

4.42 
4.79 
5.18 
5.62 
6.00 
6.30 


,  x/d  =  255.0,  $  =  0.S70 


(uin) 

GO 

947 

455 

942 

508 

932 

634 

924 

634 

894 

788 

870 

990 

835 

1229 

826 

1389 

801 

1471 

771 

1525 

717 

1594 

687 

1617 

665 

1635 

629 

1653 

582 

1677 

552 

1699 

515 

1737 

471 

1796 
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Particle  Life-History,  x/d  =  297.5,  tp  m  0.761 


E^23es*ssssnsrw 


Time  After 

All  Liquid 

Has  Evaporated 
(s) 

d 

P 

(Um) 

T 

P 

00 

0.00 

921 

461 

0.25 

916 

517 

0.41 

913 

588 

0.59 

897 

807 

0.78 

881 

939 

1.01 

867 

1003 

1.23 

864 

1094 

1.38 

847 

1198 

1-64 

844 

1244 

1.81 

821 

1317 

2.03 

810 

1364 

2.37 

794 

1417 

2.80 

767 

1473 

3.21 

754 

1495 

3.63 

729 

1509 

4.03 

691 

1517 

4.42 

664 

1539 

4.77 

652 

1546 

5.20 

624 

1559 

5.81 

583 

1580 

6.43 

530 

1607 

6.84 

483 

1626 

7.15 

451 

1642 

7.41 

420 

1651 

7.78 

401 

1657 

Particle-Life  History  x/d  =  340,  <{>  ■  0.357 


Time 

(s) 

d 

P 

(ym) 

T 

P 

(K) 

0.00 

1081 

303 

0.04 

1080 

331 

J.10 

1077 

353 

3.13 

1067 

381 

0.18 

1048 

390 

0.22 

1006 

405 

0.31 

978 

414 

0.33 

949 

415 

0.38 

915 

417 

0.41 

896 

419 

0.46 

875 

423 

0.59 

874 

451 

0.70 

873 

530 

0.82 

864 

621 

1.08 

843 

754 

1.31 

834 

887 

1.53 

817 

989 

1.64 

805 

1074 

1.96 

796 

1231 

2.17 

779 

1264 

2.34 

763 

1339 

2.66 

750 

1384 

2.98 

723 

1421 

3.27 

702 

1487 

3.69 

665 

1517 

4.11 

642 

1545 

4.61 

601 

1562 

5.13 

557 

1597 

5.78 

498 

1645 

6.16 

450 

166* 
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Particle-Life  History,  x/d  =  382.5,  <p  =  484 


Time 

(s) 

d 

P 

(pm) 

T 

P 

00 

0.00 

955 

302 

0.04 

955 

328 

0.10 

952 

366 

0.14 

945 

385 

0.17 

936 

397 

0.21 

918 

406 

0.25 

893 

410 

'0.29 

865 

412 

0.33 

836 

413 

0.36 

813 

413 

0.40 

781 

413 

0.45 

779 

415 

0.50 

778 

453 

0.73 

774 

509 

0.93 

769 

592 

1.09 

758 

659 

1.31 

755 

773 

1.63 

742 

881 

1.86 

731 

921 

2.15 

725 

969 

2.37 

712 

1008 

2.54 

703 

1046 

2.92 

693 

1089 

3.33 

685 

1123 

3.79 

670 

1154 

4.31 

661 

1174 

4.69 

647 

1183 

5.17 

638 

1196 

5.53 

630 

1202 

5.97 

604 

1203 

6.42 

581 

1208 

6.87 

575 

1210 

7.35 

551 

1211 

7.74 

533 

1213 

8.38 

504 

1215 

9.71 

477 

1217 

10.38 

442 

1220 

Particle-Life  History,  x/d  =  467.5,  4>  =  0.348 


Time 

(s) 

d 

P 

(lim) 

T 

P 

<K) 

0.00 

885 

298 

0.06 

884 

334 

0.11 

882 

358 

0.15 

877 

375 

0.19 

869 

387 

0.25 

848 

400 

0.31 

818 

405 

0.36 

787 

407 

0.41 

757 

408 

0.47 

726 

409 

0.52 

725 

408 

0.56 

725 

448 

0.90 

725 

482 

1.61 

724 

532 

2.69 

724 

615 

3.97 

723 

654 

5.77 

721 

719 

8.73 

719 

726 

10.88 

719 

729 

13.07 

717 

731 

16.70 

714 

731 

20.95 

711 

732 

24.67 

707 

733 

29.05 

704 

733 

31.73 

700 

734 

36.77 

698 

735 

40.53 

694 

734 
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Particle-Life  History,  x/d  *  510.0,  $  “  0.272 


Time 

(s) 

d 

P 

CUm) 

T 

P 

(K) 

0.00 

845.0 

295 

0.05 

845.1 

323 

0.11 

843.2 

347 

0-20 

835.4 

376 

0.29 

815.2 

394 

0.31 

809.7 

397 

0.40 

773.1 

402 

0.51 

720.3 

404 

0.59 

695.5 

405 

0.68 

695.4 

439 

1.09 

695.4 

471 

1.83 

695.4 

506 

3-88 

695.3 

586 

5.71 

695.2 

615 

8.91 

695.1 

642 

11.46 

694.9 

645 

13.64 

694.8 

648 

17.22 

694.6 

649 

21.11 

694.5 

650 

24,95 

694.3 

650 

29.00 

694.2 

652 

33.17 

693.9 

651 

36.78 

693.7 

650 

40.85 

693.6 

650 
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B. 2  Data  From  the  Lamlnar-Premixed  Flame  Tests 

A 

B.  2 . 1  Particle  Environment  Measurements 


Gas 

Temperature 

and 

Composition 

,  <J>  “  0.2 

X 

T 

Mole  Fraction 

fa=0 

(K) 

°2 

oo2 

CO 

H2° 

N2  H2 

OH 

10 

1667 

0.5314 

0.0664 

— 

0.1328 

0.2691 

0.0004 

20 

1668 

0.5314 

0.0664 

— 

0.1328 

0.2691 

0.0004 

30 

1667 

0.5314 

0.0664 

— 

0.1328 

0.2691 

0.0004 

40 

1667 

0.5314 

0.0664 

— 

0.1328 

0.2691 

0.0004 

50 

1667 

0.5314 

0.0664 

— 

0.1328 

0.2691  — 

0.0004 

Gas 

Temperature 

and 

Composition 

,  ip  =  0.6 

X 

(mm) 

T 

(K) 

Mole  Fraction 

°2 

c°2 

CO 

H2° 

N2 

H2 

OH 

5 

1662 

0.0860 

0.0637 

0.1293 

0.7208 

— 

0.0002 

10 

1663 

0.0858 

0.0640 

— 

0.1291 

0.7208 

— 

0.0003 

20 

1663 

0.0858 

0.0644 

— 

0.1291 

0.7208 

— 

0.0003 

30 

1664 

0.0858 

0.0644 

— 

0.1291 

0.7208 

— 

0.0003 

40 

1663 

0.0858 

0.0644 

— 

0.1291 

0.7208 

— 

0.0003 

50 

1664 

0.0858 

0.0644 

— 

0.1291 

0.7208 

— 

0.0003 

OH  mole  fraction  was  obtained  from  CEC-72  equilibrium  calculation. 
H^O  mole  fraction  was  obtained  from  H  and  0  element  balances. 
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Gas  Temperature  and  Composition,  4>  *  0.6 


X 

(nun) 

T 

<K) 

Mole  Fraction 

°2 

c°2 

CO 

h2° 

*2 

H2 

OH 

10 

1810 

0.0973 

0.0730 

0.0001 

0.1458 

0.6832 

0.0007 

20 

1814 

0.0973 

0.0730 

0.0001 

0.1458 

0.6832 

— 

0.0007 

30 

1815 

0.0973 

0.0730 

0.0001 

0.1458 

0.6832 

— 

0.0007 

40 

1814 

0.0973 

0.0730 

0.0001 

0.1458 

0.6832 

— 

0.0007 

50 

1814 

0.0973 

0.0730 

0.0001 

0.1458 

0.6832 

— 

0.0007 

Gas 

Temperature  and  Composition 

>  4>  -  0.6 

X 

(nm) 

T 

(K) 

Mole  Fraction 

°2 

C°2 

CO 

h2° 

N2 

H2 

OK 

10 

1952 

0.1096 

0.0822 

0.0002 

0.1639 

0.6425 

0.0001 

0.0017 

20 

1954 

0.1095 

0.0821 

0.0002 

0.1641 

0.6425 

0.0001 

0.0017 

30 

1953 

0.1096 

0.0822 

0.0002 

0.1639 

0.6425 

0.0001 

0.0017 

40 

1953 

0.1096 

0.0822 

0.0002 

0.1639 

0.6425 

0.0001 

0.0017 

50 

1953 

0.1095 

0.0822 

0.0002 

0.1640 

0.6425 

0.0001 

0.0017 
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Gas  Temperature  and  Composition,  <p  =  1.0 


x  T  _ _ Mole  Fraction 


(mm) 

(K) 

°2 

C°2 

CO 

V 

N2 

H2 

OH 

5 

1683 

0.0002 

0.0645 

0.0003 

0.1295 

0.8054 

0.0002 

0.0001 

10 

1683 

0.0002 

0.0645 

0.0003 

0.1294 

0.8054 

0.0002 

0.0001 

20 

1684 

0.0002 

0.0645 

0.0003 

0.1294 

0.8054 

0. 0002 

0.0001 

30 

1683 

0.0002 

0.0645 

0.0003 

0.1295 

0.8054 

0.0002 

0.0001 

40 

1683 

0.0002 

0.0645 

0.0003 

0.1295 

0.805^ 

0.0002 

0.0001 

50 

1683 

0.0002 

0.0645 

0.0003 

0.1295 

0.8054 

0.0002 

0.0001 

60 

1683 

0.0002 

0.0645 

0.0003 

0.1295 

0.8054 

0.0002 

0.0001 

70 

1681 

0.0002 

0.0645 

0.0003 

0.1295 

0.8054 

0.0002 

0.0001 

80 

1678 

0.0002 

0.0646 

0.0003 

0.1295 

0.8054 

0.0002 

0.0001 

90 

1675 

0.0002 

0.0646 

0.0002 

0.1295 

0.8054 

0.0002 

0.0001 

100 

1668 

0.0002 

0.0646 

0.0002 

0.1295 

0.8054 

0.0002 

0.0001 

110 

1665 

0.0002 

0.0646 

0.0002 

0.1295 

0.8054 

0.0002 

0.0001 

120 

1661 

0.0002 

0.0646 

0.0002 

0.1295 

0.8054 

0.0001 

0.0001 

130 

1657 

0.0002 

0.0646 

0.0002 

0.1295 

0.8054 

0.0C01 

0.0001 

140 

1648 

0.0002 

0.0646 

0.0002 

0.1295 

0.8054 

0.0001 

0.0001 

150 

1640 

0.0001 

0.0646 

0.0002 

0.1295 

0.8055 

0.0001 

0.0001 

160 

1625 

0.0001 

0.0646 

0.0002 

0.1295 

0.8055 

0.0001 

0.0 

170 

1606 

0.0001 

0.0647 

0.0001 

0.1295 

0.8055 

0.0001 

0.0 

180 

1597 

0.0001 

0.0647 

0.0001 

0.1295 

0.8055 

0.0001 

0.0 

190 

1578 

0.0001 

0.0647 

0.0001 

0.1296 

0.8055 

0.0001 

0.0 

200 

1564 

0.0001 

0.0647 

0.0001 

0.1296 

0.8055 

0.0001 

0.0 

Gas  Temperature  and  Composition,  $  *  1.0 


x  T  _ Mole  Fraction 


(run) 

(K) 

°2 

co2 

CO 

V 

N2 

H2 

OH 

10 

1816 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

0.0002 

20 

1819 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

0.0002 

30 

1819 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

0.0002 

40 

1820 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

0.0002 

50 

1819 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

0.0002 

60 

1819 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

0.0002 

70 

1817 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

0.0002 

80 

1817 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

0.0002 

90 

1815 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.006'? 

0.0002 

100 

1812 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

0.0002 

110 

1809 

0.0005 

0.0722 

0.0007 

0.1456 

0.7806 

0.0004 

0.0002 

120 

1805 

0.0005 

0.0723 

0.0007 

0.1455 

0.7806 

0.0004 

0.0002 

130 

1801 

0.0005 

0.0723 

0.0007 

0.1455 

0.7806 

0.0004 

0,0002 

140 

1793 

0.0005 

0.0723 

0.0006 

0.1456 

0.7806 

0.0003 

0.0002 

150 

1785 

0.0009 

0.0724 

0.0006 

0.1456 

0.7806 

0.0003 

0.0002 

160 

1773 

0.0004 

0.0724 

0.0006 

0.1457 

0.7806 

0.0003 

0.0001 

170 

1754 

0.0004 

0.0725 

0.0005 

0.1456 

0.7807 

0.0003 

0.0001 

180 

1742 

0.0003 

0.0725 

0.0004 

0.1458 

0.7807 

0.0003 

0.0001 

190 

1730 

0.0003 

0.0726 

0.0004 

0.1457 

0.7808 

0.0002 

0.0001 

200 

1718 

0.0003 

0.0726 

0.0004 

0.1457 

0.7808 

0.0002 

0.0001 
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Gas  Temperature  and  Composition,  <j>  =  1.0 


x  T  _ .  Mole  Fraction 


(nm) 

(K) 

°2 

C°2 

CO 

h2° 

*2 

H2 

oil 

10 

1949 

0.0013 

0.0806 

0.0019 

0.1641 

0.7509 

0.0009 

0.0006 

20 

1952 

0.0013 

0.0806 

0.0020 

0.1641 

0.7509 

0.0009 

0.0006 

30 

1953 

0.0013 

0.0806 

0.0020 

0.1641 

0.7509 

0.0009 

0.0006 

40 

1952 

0.0013 

0.0806 

0.0020 

0.1641 

0.7509 

G.0Q09 

0.0006 

50 

1952 

0.0013 

0.0806 

0.0020 

0.1641 

0.7509 

0.0009 

0.0006 

60 

1951 

0.0013 

0.0806 

0.0019 

0.1642 

0.7509 

0.0009 

0.0006 

70 

1950 

0.0013 

0.0806 

0.0019 

0.1642 

0.7509 

0.0009 

0.0006 

80 

1948 

0.0013 

0.0806 

0.0019 

0.1642 

0.7509 

0.0009 

0.0006 

90 

1946 

0.0012 

0.0806 

0.0019 

0.1644 

0.7509 

0.0009 

0.0005 

100 

1944 

0.0012 

0.0807 

0.0019 

0.1643 

0.7509 

0.0009 

0.0005 

no 

1940 

0.0012 

0.0807 

0.0018 

0.1643 

0.7510 

0.0009 

0.0005 

120 

1935 

0.0012 

0.0808 

0.0018 

0.1643 

0.7510 

0.0008 

0.0005 

130 

1930 

0.0011 

0.0808 

0.0017 

0.1644 

0.7510 

0.0009 

0.0005 

140 

1924 

0.0011 

0.0809 

0.0016 

0.1644 

0.7511 

0.0008 

0.0005 

150 

1913 

0.0C10 

0,0810 

0.0015 

0.1745 

0.7511 

0.0007 

0.0004 

160 

1901 

0.0010 

0.0811 

0.0014 

0.1645 

0.7512 

0.0007 

0.0004 

170 

1891 

0.0009 

0.0812 

0.0013 

0.1645 

0.7512 

0.0007 

0.0004 

ISO 

1879 

0.0008 

0.0813 

0.0012 

0.1647 

0.0006 

0.0006 

0.0003 

190 

1860 

0.0008 

0.0814 

0.0011 

0.1647 

0.7514 

0.0006 

0.0003 

200 

1851 

0.0007 

0.815 

0.0010 

0.1648 

0.7614 

0.0006 

0.0003 

Gas  Temperature  and  Composition,  $  *  1.4 


x  T  _  _ Mole  Fraction 


(tsa) 

(K) 

°2 

C°2 

CO 

H2° 

N2 

«2 

OH 

5 

1622 

— 

0.0421 

0.0534 

0.1432 

0.7092 

0.0542 

— 

10 

1624 

— 

0.0427 

0.0540 

0.1412 

0.7092 

0.0563 

— 

20 

1625 

— 

0.0428 

0.0542 

0.1409 

0.7092 

0.0566 

— 

30 

1625 

— 

0.0428 

0.0543 

0.1408 

0.7092 

0.0566 

— 

40 

1624 

— 

0.0427 

0.0542 

0.1409 

0.7092 

0.0566 

— 

30 

1624 

— 

0.0427 

0.0341 

0.1410 

0.7092 

0.0566 

— 

60 

1624 

— 

0.0428 

0.0541 

0.1410 

0.7092 

0.0567 

— 

70 

1621 

— 

0.0428 

0.0541 

0.1410 

0.7092 

0.0566 

— 

80 

1620 

— 

0.0429 

0.0541 

0.1408 

0.7092 

0.0566 

— 

90 

1618 

— 

0.0430 

0.0541 

0.1407 

0.7092 

0.0567 

— 

100 

1612 

— 

0.0430 

0.0540 

0.1407 

0.7092 

0.0568 

— 

110 

1609 

— 

0.0430 

0.0540 

0.1406 

0.7092 

0.0569 

— 

120 

1604 

— 

0.0431 

0.0538 

0.1406 

0.7092 

0.0569 

— 

130 

1501 

— 

0.0432 

0.0537 

0.1405 

0.7092 

0.0571 

— 

140 

1594 

— 

0.0434 

0.0536 

0.1403 

0.7092 

0.0572 

— 

150 

1589 

— 

0.0435 

0.0535 

0.1402 

0.7092 

0.0573 

— 

160 

1570 

— 

0.0439 

0.0532 

0.1397 

0.7092 

0.0577 

— 

170 

1557 

— 

0.0442 

0.0529 

0.1394 

0.7092 

0.0580 

— 

180 

1543 

— 

0.0444 

0.0525 

0.1393 

0.7092 

0.0583 

— 

190 

1522 

— 

0.0449 

0.0520 

0.13S8 

0.7092 

0.05S8 

— 

200 

1507 

— 

0.0450 

0.0519 

0.1387 

0.7092 

0.0539 

— 
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